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The IRE “Affiliate” Plan—A New Venture in 
Engineering Society Structure and Service 
W. R. G. BAKER, Chairman IRE Professional Groups Committee 


On January 4, 1957 the IRE Board of Directors 
arrived at a decision which may in time prove to 
be one of the most far-reaching in its 45-year 
history. On that date the Board adopted a plan 
which will enable non-I[RE members whose main 
professional interests lie outside the sphere of [RE 
activities to become affiliated with certain of the 
IRE Professional Groups without first having to 
join the IRE itself. 

This plan is aimed at those specialists in other 
fields of science and technology whose work 
touches upon our own electronics and communica- 
tions field only in specialized areas. In effect, the 
IRE is extending the specialized services of its 
Professional Groups to every field of science and 
engineering. 

An outstanding example of where these services 
are needed may be found in the case of the medi- 
cal and biological sciences. At the present time 
some 1400 IRE members enjoy the privileges of 
membership in the Professional Group on Medical 
Electronics. And yet there are hundreds, perhaps 
thousands, of medical doctors, biologists, and 
others to whom the activities of this Group would 
be of interest and value. Both they and the Group 
would benefit from their participation. To require 
these persons, who have no interest in radio engi- 
neering, to join the IRE in order to join the Group 
is unreasonable, and probably futile as well. In 
fact, it was largely to provide an answer to this 
particular problem that the “Affiliate” Plan was 
first conceived, although it pertains to other fields 
as well, such as Computers, etc. 

The “Affiliate” Plan is admittedly an experiment. 
So far as is known, no other society has ever tried 
a similar scheme. The Board of Directors feels 
strongly that the benefits afforded by the plan 
justify the risk that some persons who should join 
the IRE will instead become Affiliates. To mini- 
mize this risk, the plan has been carefully worked 
out along the following lines: 

1) Participation in the Plan is at the option of 
each Professional Group. It is not expected that 
all Groups will adopt it; only those which feel it 
serves a need in their particular field. 


2) Each Group interested in initiating the “Affil- 
iate” Plan must submit to the Chairman of the 
Professional Groups Committee a list of accredited 
organizations which has been selected and ap- 
proved by its Administrative Committee, for 
official approval by the IRE Executive Commit- 
tee; 

3) To be an Affiliate of a Professional Group, a 
person must belong to an accredited organization 
approved by that Group and the IRE Executive 
Committee. Moreover, he shall not have been an 
IRE member during the five years prior to his ap- 
plication. He may affiliate with more than one 
Group, provided the accredited organization to 
which he belongs is recognized by the Groups con- 
cerned. 

4) The fee for Affiliates shall be the assessment 
fee of the Group, plus $4.50. The latter covers 
IRE subsidies to the Group, Professional Group 
overhead expenses borne by IRE Headquarters, 
and 50 cents which is to be rebated to IRE Sec- 
tions for mailing and meeting costs. 

5) An Affiliate will be entitled to receive the 
TRANSACTIONS of his Group and that part of the 
IRE NatTIonaL CONVENTION RECORD pertaining 
to his Group. He will be eligible for a Group award, 
and may attend local or national meetings of 
the Group by payment of charges assessed Group 
members. 

6) An Affiliate cannot serve in an elective office 
in the Group or Group Chapter, nor vote for candi- 
dates for these offices. 

7) An Affiliate may hold an appointive office in 
the Group or Group Chapter. 

8) An Affiliate may not receive any IRE bene- 
fits that are derived through IRE membership, ex- 
cept as authorized by the Executive Committee. 

The “Affiliate” Plan is a bold and farsighted 
venture; one that recognizes and provides for the 
rapidly spreading influence of electronics in every 
walk of scientific and technological life, and one 
that enables the IRE to further its aims as a 
professional engineering society—the advancement 
of radio engineering and related fields of engi- 
neering and science. 
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A Time-Sequential Tabular Analysis of 
Flip-Flop Logical Operation* 


GENE W. 


Summary—In examining flip-flop response the principal concern 
of the logical designer is to find what input signals must be applied 
to the flip-flop in order to produce the output conditions that are 
desired. Equation methods of analysis and a time-sequential tabu- 
lar method of analysis are described, and some advantages of the 
tabular method are pointed out. 


INTRODUCTION 


RATHER common problem is the logical design 
of a digital computer network, consisting of flip- 


flops and gates, which is to be operated in a se- 
quential cycle including a number of separate successive 
time intervals.! The flip-flops drive the gates which in 
turn supply input signals to the flip-flops, and although 
other input and output signals are involved in the net- 
work their analysis presents no unique problem. The 
network must exhibit the correct logical response to in- 
formation conditions which may exist during any time 
interval; hence, response requirements for the various 
time intervals may be lumped together for design pur- 
poses as constituting a single set of requirements. Either 
Boolean equations or tabular methods of analysis may 
be used. 

One important phase of the problem may be described 
with reference to a single flip-flop. There is provided as 
the starting information a statement of the response 
which the particular flip-flop is required to exhibit to 
each of various possible information conditions. It then 
becomes the task of the logical designer to design gates 
which will supply correct input signals to the flip-flop. 

Although the usual description? of flip-flop logical 
operation is a table or equation specifying the observed 
responses to various input conditions, it is apparent that 
what the logical designer needs is essentially an inverse 
representation of the flip-flop performance character- 
istic. That is, the logical designer needs a representation 
ot the flip-flop indicating the input signals which must 
be supplied in order to produce a desired effect upon the 
stable state of the flip-flop. 


* Manuscript received by the PGEC, October 26, 1955; revised 
manuscript received, December 22, 1956. The paper describes theoret- 
ical developments with which the author became familiar while em- 
ployed as a patent attorney at Hughes Aircraft Co., Culver City, 
Calif., and The Ramo-Wooldridge Corp., Los Angeles, Calif. 

i Patent lawyer, 19 Pine Ave., Long Beach, Calif. 

1J.R. Harris, “A transistor shift register and serial adder,” 
IRE, vol. 40, pp. 1597-1602; November, 1952. 

2'E. C. Nelson, “Algebraic theory for use in digital computer de- 
sign,” IRE TRANS., vol. EC-3, pp. 12-21; September, 1954. 

3 J. A. Postley, “A method for the evaluation of a system of Boo- 
lean algebraic equations,” Mathematical Tables and Other Aids to 
Computation, vol. 1X, 49, pp. 5-8; January, 1955. 
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EQUATION ANALYSIS 


The functional requirements placed upon a flip-flop in 
any time interval may be represented in a generalized 
form by the following: 


OF — Ort eed ae On aVn—1 (e 


Another way of saying it is that (1) indicates the infor- 
mation required to be represented by the stable state of 
the flip-flop. In (1) Qn_1 represents a binary digit, or bit, 
of information which is presently stored in the particular 
flip-flop and which has possible values of 1 and 0. O,-1 
represents the binary complement of Qn. X,1 and 
Y,1 are functions of other information then available 
in the network, and Q, is the new bit of information 
which must be formed by the gates and stored in the 
flip-flop in the ensuing time interval. 

Of particular interest is the “triggering” flip-flop 
which has two separately triggerable input circuits, and 
which reliably responds to the simultaneous pulsing of 
both inputs to change from one to the other of its stable 
states. Solving (1) and the flip-flop response equation‘ 
simultaneously, several solutions for the input functions 
may be obtained as follows, where J and K respectively 
represent signals applied to the two input circuits: 


—1 = On = On Xt ain OLIV 
iene a On aed Ona + Onn ears 


Eq. (2) implies the generation of the new bit of informa- 
tion entirely within the gates and its subsequent trans- 
fer from the gates to the flip-flop. In the physical realiza- 
tion of (2) the function for J may be generated by one 
gating circuit while the function for K is simultaneously 
generated by a second and separate gating circuit; or 
only a single function may be developed in a logical gate 
and then supplied to a complementer device for produc- 
ing a pair of signals respectively representing the func- 
tion and its complement. Note that insofar as the flip- 
flop itself is concerned, (2) also corresponds to in- 
formation shifting as normally accomplished within a 
shifting register. | 


(2) 


1 >= Optatin 


He 3 
Keay = On vA \ ' 


Eq. (3) is a simplification of (2), however, the functions 
for J and K are no longer mutually complementary. The 


* Nelson, op. cit., Postley, op. cit. 
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physical realization of (3) therefore requires two sepa- 
rate gating circuits. 


ae ~~ NESs 


Pon 4 
jen = b. Grare \ ) 


Eq. (4) is a further simplification® in which use is made 
of the fact that the flip-flop responds with knowledge of 
its existing state. The physical realization of (4) also re- 
quires two separate gating circuits except in special cases 
where either one or both of the functions X¥ and Y isa 
function of a single binary variable and hence directly 
available without the need for being synthesized in a 
special gate. 


Pes = Tope bag Ort 4 aia Oa aes (5) 


Eq. (5) represents the special case where the two input 
circuits are physically connected together and con- 
trolled by a single input function. 

Also of interest is the “set-reset” flip-flop which is 
like the “triggering” flip-flop except that it does not 
reliably respond to the simultaneous pulsing of both 
its inputs and must, therefore, be subject to the restric- 
tion that only one of its inputs may be pulsed in any 
time interval. Input functions according to (5) can- 
not, therefore, be used, and input functions according 
to (4) can be used only if Y and X are never simultane- 
ously 1 (as for example, in the carry flip-flop for a 
binary adder or subtracter). The set-reset flip-flop may 
be converted to a triggering flip-flop by an appropriate 
feedback gate connected from each of its outputs to 
the opposite input. Input functions may then be sup- 
plied in accordance with (3), which corresponds to the 
simplified input functions of (4) plus the feedback 
terms. 


TABULAR ANALYSIS 


If the functional requirements represented by (1) are 
stated in tabular form it becomes convenient to repre- 
sent the flip-flop characteristic in tabular form also. 
Thus Table I specifies for the “triggering” flip-flop the 


/2ZcABEBSI 
NECESSARY AND SUFFICIENT INPUT SIGNALS FOR THE TRIGGERING 
FLip-FLop 
Ore Qn SS Te 
0 0) 0 
0 1 1 
1 0 1 
il j 0 


input signals required to produce any desired effect up- 
on its stable state. Symbol Q representing the stable 
state of the flip-flop has two possible values designated 
as binary 1 and 0, while J and K respectively represent 


5 First proposed at Hughes Aircraft Co. by Dr. Eldred C. Nelson 
in 1951. 
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signals applied to the two input circuits and are as- 
signed binary values representing pulse or no pulse. 
The independent variable is the desired change (or no 
change) in stable state. Blank spaces in the J and K 
columns indicate that either a binary 1 or 0 may be 
used. Table I has three other equivalent forms obtain- 
able by changing the binary designation of the stable 
state, the input signals, or both. It will be noted that 
Table I is not merely the inverse of the usual description 
of flip-flop response, but is somewhat more simplified. 

The use of the time-sequential tabular method of 
analysis in which the tabular statement of functional 
requirements and the tabular description of flip-flop 
response are effectively combined® is illustrated by 
Table II, an example of a counter which continuously 
cycles through the 5 states—000—110—010—001—100. 


TABLE. It 


Time Intervals THgeorne 
ro) 


Flip-Flop Input Signals 


n—1 n 
QA |} QB|QC|Q4|Q0B|0C|JA|KA|JB|KB| JC| KC 
0 0 0) 1 1 0 1 1 0 
1 1 0 0 1 0 1 0 0 
0 1 0) 0 0 1 0 1 1 
0 0 1 1 (0) 0 1 0 1 
1 0 0 0 0 0 1 0 0 


The zeros and ones indicate the states of the three 
triggering flip-flops QA, QB, and QC in each count. The 
following equations are easily derived: 


VAL=< 

KA =" Tord or ©) 

J Be= AC 

KB =A 

JO= AB 

KC: =\“1" (ord or B orC) 


where the time subscripts have been omitted since only 
simultaneous gating’ is involved. 


COMPARISON OF METHODS 


In general, the choice of a method of analysis is 
influenced by whether or not all of the possible com- 
binations of values of the information variables from 
which the input functions are to be formed may exist 
and need to be considered. In other words, the question 
is whether the computer must be designed to respond to 


6 This method was first used at Hughes Aircraft Co. by John V. 
Blankenbaker in 1953. Mr. Blankenbaker has also supplied the par- 
ticular illustration shown in Table II. 

7™W. V. Quine, “The problem of simplifying truth-functions,” 
Amer. Math. Monthly, vol. 59, pp. 521-531; October, 1952. 

8 R. K. Richards, “Arithmetic Operations in Digital Computers,” 
D. Van Nostrand Co., Inc., New York, N. Y., pp. 51-80; 1955. 
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all possible information conditions, or to less than all. 

If all possible information conditions must be con- 
sidered in the analysis, the simplest possible solutions 
for the flip-flop input functions may be obtained either 
by the described tabular method or by the described 
equation method. The tabular method has the advan- 
tage of simplicity and may for that reason be preferred. 

If the analysis includes less than all of the possible 
information conditions the tabular method will again 
yield the simplest possible solutions. Described equation 
method may not do so, however, as is shown by deriving 
the input functions for flip-flop QB of the example. A 
general expression |like (1)| for the required state of 
flip-flop QB may be written from Table IT: 


OR p=) Aya Bn Gt Asi Bg Cnt 

St Bd eet: B= pA em ey 

and in accordance with (1): 
a 

Y= 


AC 
AC 
where the time subscripts have again been dropped. 


Hence, the simplified input functions according to (4) 
become: 
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re 
AR 6, 


JB Y 


KBeavx = 
It is apparent that the expression for KB derived by 
the tabular method is simpler as it includes only the 
term A. 

Where less than all of the possible information condi- 
tions need to be considered, the simplest possible solu- 
tions for the flip-flop input functions may nevertheless 
be obtained by means of a more elaborate equation 
analysis. The procedure involves writing one equation 
for all of the information conditions to which a re- 
sponse is required and an additional equation for each 
of the other information conditions, the equations then 
being solved simultaneously. On a comparative basis, 
however, the described tabular analysis appears to be 
a far more efficient working tool. 
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Dynamic Accuracy as a Design Criterion of Linear 


Electronic-Analog Differential Analyzers* 
AMOS NATHANT 


Summary—A frequency error analysis of computing elements is 
presented which leads to a definition of their dynamic accuracy. 

The concept of a computing transfer function is introduced for this 
purpose, permitting the evaluation of an effective bandwidth, the lat- 
ter being connected with the variance of the output for wideband in- 
puts. Limited bandwidth is considered as equivalent to finite resolu- 
tion and thus to an additional effective error. Single frequency errors 
are dealt with separately and are shown to be of minor importance. 
Suitable optimalization of dynamic accuracy yields parameters of 
design and performance such as optimum computing time and re- 
quired base amplifier gain. 

The theory is applied to integrators and adders with base ampli- 
fiers of direct and of capacitive coupling. 


INTRODUCTION 


HE BASIC specifications of differential analyzers 
consist of the required accuracy and range of 
computing times. The problem of the dynamic 
accuracy of a computing setup for the solution of linear 
differential equations has been solved by Macnee}? 


* Manuscript received by the PGEC, December 10, 1955; re- 
vised manuscript received, March 1, 1957. 

+ Israel Inst. Tech., Haifa, Israel. 

1 A.B. Macnee, “An electronic differential analyzer,” Proc. IRE, 
vol. 37, pp. 1315-1324; November, 1949. 

2 A. B. Macnee, “Some limitations on the accuracy of electronic 
differential analyzers,” Proc. IRE, vol. 40, pp. 303-308; March, 
1952. 


who considers the perturbations of the roots of the 
differential equation which are caused by the physical 
characteristics of the computing elements. His theory 
is invaluable in the direct prediction of errors for given 
setups. For all but the simplest problems, however, it 


yields qualitative results only. Representative setups; 
1.€., representative differential equations, determine 


criteria of design. It is somewhat difficult to say which 
problems are typical and an alternative approach may 


be desirable. The nonlinear case has also been dis- 


cussed.?4 


The concept of dynamic accuracy as commonly used | 


is rather vague. Here the performance of isolated com- 


puting elements operating upon representative inputs is | 


considered and the error associated with the effective 
bandwidth, as recently introduced,’ is used to define dy- 


3F. J. Murray, “Mathematical Error Analysis for Continuous 
Computers,” Project Cyclone Symposium II, Part 2, Reeves Instru- 
ment Corp., New York, N. Y., pp. 139-146; 1952. 

4K. S. Miller and F. V. Murray, “A mathematical basis for an 
error analysis of differential analyzers,” J. Math. Phys., vol. 32, pp. 
136-163; July—October, 1953. 

° A. Nathan, “A note on bandwidth,” Proc. IRE, vol. 44, pp. 788- 
790; June, 1956. (See also “Correction,” Proc. IRE, vol. 45, p. 65; 
January, 1957.) 
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namic accuracy. It is shown in principle and explicitly 
in examples how to determine this quantity. 

Once the applicability of the definition of dynamic 
accuracy has been accepted the importance of its 
determination is clear; parameters of design required 
to meet given specifications can be based on it and its 
optimalization yields information about the desirable 
value of fundamental quantities such as the amplifica- 
tion and the gain-bandwidth product of base amplifiers. 


REPRESENTATIVE INPUTS 


A computing element has to operate upon a variety 
of input signals. The errors introduced by its imperfec- 
tions cannot be determined unless these signals be 
specified or the complete setup of the computing scheme 
be known. Some essential arbitrariness is seen to be 
involved at the very outset of any general theory of the 
accuracy of computing elements. At the basis of such 
a theory either some setups or else some inputs must 
be regarded as representative, and its success. will 
largely depend upon the judicious choice of these. It 
must be borne in mind that cases may arise which devi- 
ate significantly from these assumptions; results of the 
theory must therefore not be taken too literally. 

We intend to develop a theory based on consideration 
of typical inputs. Now such inputs may be represented 
in the time or in the frequency domain. The time de- 
pendent input represents directly the variable on which 
a mathematical operation is performed by the element. 
Indeed, time is the only independent variable in elec- 
tronic analog computers; frequency representation is 
more convenient in the description of computing ele- 
ments. Thus, a transfer function is physically more read- 
ily interpreted than equivalent impulse response. 

In transforming the input from time to frequency 
dependence, note that knowledge of the march of the 
complete time function is required, from t=— © to 
t=-+ 0. The computing time per solution of a differ- 
ential analyzer is limited to t=0 - - - T, and only those 
values of the input that lie between ¢= —e and t=T, 
where e€ is a small positive quantity, should affect the 
output during 0--- 7. Although later we shall adopt 
a slightly different approach, it will now be convenient 
to assume zero input for <0 and for t>T. 

Let us, therefore, examine several truncated functions 
which are to be considered as examples of inputs to 
computing elements. Table I (next page) contains sev- 
eral such functions and their power spectra. Here we 
have set 7 =1. The spectra are seen to fall roughly into 
one of two groups: 1) wide-band spectra which resemble 
in character that of the sampling function, no. 4) of the 
fable, [1), 2), 3),-and 4)|, and 2) those that have one 
pronounced peak [5) and 6)]; the latter resemble the 
spectrum of a continuous sine wave). It will be shown 
that the requirements imposed by the first group of sig- 
nals are usually much more severe than those of the sec- 
ond. Criteria of design and performance will, therefore, 
be based primarily on inputs of the first group. 
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WIbE-BAND INPUTS 


The performance of linear computing elements with 
wide-band inputs can be described by their effective 
bandwidth, a concept which was introduced in a recent 
paper,® and which serves to compare the operation of 
a given network NV with that of an “ideal” network J, 
with this type of input. The general idea is as follows. 

Take the band-limited input signal® 


Nee (clea | w | SS 
é;= 8, = (1a) 
Q, | w | Gs 
which corresponds to 
Wm sin W®m(t — to) 
e(t — to) = ex(t — to) = 
T Wm (t — to) 
A Om ., 
= SINC wm(t — to); (1b) 


1 


and use it as input to both N and J. (In the reference 
only fo =0 was considered.) In comparing the outputs 
we permit the adjustment of the gain of an element, 
in order to obtain the best agreement; two elements 
would, for example, be considered as equivalent if the 
ratio of their transfer functions is a real constant not 
necessarily equal to one. Thus, we consider the transfer 
function of J as given, except for a multiplicative real 
constant a. The variance of the outputs provides a 
a measure of their agreement and its minimum with 
respect to a is defined as the effective—or rms error 
squared, «7. It turns out to be 


€12(wm 


9 


[| Theae ih | Hy [*do — | fl Re i y*)de | 
é ff | Tieltae [| Fy Paw 


where the integrations extend from 0 to @m. Hz, Hw 
are the transfer functions of J and JN, respectively, and 
w» is defined as the effective (angular) bandwidth of N 
with respect to J associated with error €. 

In the limit of infinite w,, the input is proportional to 
a delta function and (2a) can then be expressed in terms 
of the impulse responses W7(t) and Wy(t) of J and N, 
yielding 


(2a) 


D) 


{ wewa f werenat—| f TOC 
2 0 = : 0 ~ (2b) 
{ weoa f Ww?(t)dt 


6 These symbols will be used throughout this paper: 
= identically equal % proportional 
4 equal by definition 6(t) unit impulse 
~ approximately equal 1(¢) unit step, 

~ very roughly equal 
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TABLE '! 


SPECTRA OF TRUNCATED FUNCTIONS 


lel 


1(7)1(1—r) 


1(r—7) 1(r2—7) the yaya, 
— sin —— 
ee eee e 
at 
%~, w 
1 
w 
Se nor) —{eé = 1(+0m)1 (om) 
&m(T—TO9) i 
(not truncated) 
Dey, wW 
sin 27N71(7)1(1—r) ; 
2w sin — 
(N =2 sketched) sy 
(277N)?—w? 
cos 2xN7r1(7)1(1—7) 4aN sin — 
(N =2 sketched) Pee Sere. see 
(27N)2—w? 
2w &W oF 8 Ww 
SINGLE-FREQUENCY INPUTS frequency domain. The complex outputs of J and N 


The other extreme consists of a steady-state sine are now Hr(jwo) and Hwy(joo). Writing aH; for Hy as 
input, €’#°', corresponding to an impulse 6(w—wo) in the before, the output error is E=aH,;— Hy (Fig. 1). This 
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E 


sirgeg JN 


expression attains its minimum value for a| H,| a4 Hy| 
provided the fractional minimum error is small, which 
is then given by 


(3) 


z.e., the fractional single-frequency error is approxi- 
mately equal to the phase error of N with respect to 
I. As such it is usually specified for amplifiers and for 
adders only; we now recognize its general importance. 
Note @ is a function of frequency. What now, if two 
frequencies, say #, and w:, must be handled simultane- 
ously by the computing elements, or if we do not allow 
a readjustment of gain for error compensation when 
changing the input frequency from a; to w.? If we ad- 
just a for minimum error at wo =a@,, which requires 


a(w:) | Ar(w1)| = | Hy(o:) | 


(3) will apply at this frequency. At w2, however, the 
error will no longer be given solely by the phase differ- 
ence, since, in general, a(w:)| Hr(w2)| | H(w2)|. Thus, 
the error for w. will be larger than that given by (3). 
These remarks show the limited usefulness of the pres- 
ent results and the advantage of the simultaneous con- 
sideration of a whole band of frequencies as in the pre- 
ceding section, where equal weight was given to all 
frequencies in the band 0 to wm. 

It will be convenient to normalize the time scale at 
this point before proceeding to apply these concepts 
to differential analyzers. 


CHANGE OF SCALE 


Let us change the time scale from real time f to a 
normalized time 7 defined by 


p Bee 


The computing interval in 7 is thus always 0--- 1. 
We have, for example, e;(r) =e;,(77), where the suffix 
t denotes quantities on the old scale. 

Letting w denote the angular frequency in T, w=wT7. 
The scale of the impulse response may be changed as 
well, writing 


Wr) 2 TWAT,) 
in order to ensure that for the outputs 


€o(T) = €o,(77). 
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These relations affect the transfer function as follows, 
H(p) = H.(P/T). 


The last two relations are proved in Appendix I. 
From now on we shall mostly work with normalized 
time. 


THE COMPUTING TRANSFER FUNCTION 


The bandwidth relation (2) is not directly applicable 
to differential analyzers where we are interested in the 
errors which occur during the computing time only; 
only that part of signals which lies within 7=0--: - 1 
must be considered. The assumed input can be made 
to meet this requirement well enough; writing w,f=7T 
it is centered on To and goes to zero fairly rapidly both 
before and after; it will, therefore, fall almost wholly 
within or near the allowed interval, provided 07) 1. 
This need not be so for the output. An ideal integrator, 
for example, will hold the integral of a delta function 
input for all eternity, yet it is only the truncated; 
4.e., T=0--+-+1, part of the output of element JN that 
must be compared with that of J in forming the vari- 
ance. Transfer function H, however, may provide non- 
zero output at 7>1, through Hé;= 2. 

Let us replace H(p) by the computing transfer func- 
tion C(p) which will now be defined. First, note that an 
element can be described either by its transfer function 
H or by its memory function, or impulse response, 
Wir). In fact W=H. We now define the computing 
memory c(r) as the truncated function 


cr) 2 WHIP) (4) 


and the computing transfer function €(p) as its Fourier or 
Laplace transform. W(r) =0 for r<0; c(r) =0 for r<0 
and fort >1. The computing memory extends over unit 
duration. Provided the input is confined within 7 
=0---1, C will generate the same output as 7 during 
this interval and may generate some additional output 
during r=1.-- - 2; the rest is silence. Formally, using 
the convolution integral, with fi(r) =e,;(7)A(7)A(1 —7) 
as input and (4) as memory, the output becomes 


é)(7r) = {noe — t)dt 


feria NW Gray Mee oe gs 


( 0, cara OF 
f{ «@ W(r — f)dt, NS okesed: 
= i 
ih OW, —f)dt, 15752: 
0, r2=2. 


Returning now to the bandwidth relations (2), writing 
e;(t —To) for e;(t—to), T for Wmt, To for Wmeo, and w for w 
in (1) and (2), and noting that practically all of input 
(1) occurs at or near T=70, we may replace the condi- 
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tion of confinement of e; within r=0 +--+ 1by0S7)S81. 
We shall now use (2a) and (2b) for the calculation of 
errors in time limited computations by substituting C 
for H and ¢ for W. It is easy to see that this procedure 
is accurate, if tT) =0 in test signal (1). Actually there is 
little justification for this additional requirement since 
the test signal could be centered on any instant within 
the computing interval with equal likelihood, but it is 
not considered worthwhile to rectify this deficiency, 
although this could be done by averaging over-all To 
from 0 to 1. In this connection note that requirement 
Ty =0 is superfluous in the consideration of computing 
elements of short-range memory such as adders or 
differentiators. 

The relation between computing transfer function C 
and ordinary transfer function H is (Appendix II) 


C(p) = H(p) — ? [Wr + 1)1(7)]. (5) 


As examples, consider ideal adders and integrators. 


For an adder W;(t) =6(¢), thus 


W(r) = To(T7) = 8(7); Gr) =. (7): 
and 
C= 1. (6) 
For an integrator 
Ve) PT ce 
and 
W(r) =f1(Tr) = T1(e)5 Lam hy pa 
Thus 
Gr) ="L1@)1a —'7) 
and 
aay aa se | C| ae (7) 
p Ww 2 
Note that C(0)| =T, whereas | H(0)| is infinite. Fig. 2 


shows the absolute value of H and C for ideal integra- 
tors. 

We are now in a position to calculate the effective 
bandwidth B,, as a function of the error «. The presen- 
tation of this relationship is of importance in order to 
determine the best computing time of a setup. It is 
sometimes possible to estimate the salient frequencies 
expected in the solution of a problem; alternatively, once 
an approximate solution of a problem has been ob- 
tained, the spectra of the outputs which appear at all 
the computing elements can be estimated. This informa- 
tion allows an adjustment, or readjustment, of com- 
puting speed so that these spectra will fall into regions 
of small error; roughly speaking, computing speed 
should be such that signal frequencies fall within the 
pass band of the computing elements. It is suggested 
that curves of bandwidth vs error be supplied for all 
elements. 


June 


6) 


Fig. 2—Ordinary and computing transfer functions of an ideal 
integrator. 


Dynamic Accuracy (DA) 


Returning to our band-limited input test signal, & 
would provide a direct measure of DA; yet actual spec- 
tra may extend to infinity and we must now turn to a 
consideration of the effects of high-spectrum frequen- 
cies. Alternatively, we pose the question whether there 
is a criterion prescribing the required bandwidth. 

The situation may be regarded as the cascade con- 
nection of two elements: VW, a filter of bandwidth B,, 
and WN, the actual computing element. M introduces 
errors by suppressing frequencies; NV, through inaccurate 
operation within B,,. The fractional rms error intro- 
duced by M will be denoted by e&; it depends on the 
nature of the input and on B,,. &, as given by (2a), de- 
notes the error of V. We define the fotal error as 


ne 
Etot = Ver? + €22. (8) 


Although €, could be determined by a consideration 
of an appropriate input test signal, such as a step func- 
tion, a different approach will be preferred. If M were 
a filter with sharp cutoff, its (numerical) bandwidth B,, 
would allow it to pass without error 2B, sampling 
pulses per unit time; 7.e., per computation. In this sense 
we may speak of a resolution 1/2B,, of WM. We assume 
that € is proportional to this resolution and write 


€2 = B/2Bm. (9) 


With a resolution of 1/2B,, one can associate a maxi- 
mum error squared of (+1/4B,,)? which corresponds to 
an average of 3(1/4B,,)?. Equating this to €? one obtains 
from (9) 8=0.29. Let us use 8 =0.20 instead, thus allow- 
ing for the fact that M is a filter which cuts off gradually 
rather than sharply. Therefore 


2? = 10-2/Bn? = 0.40/wm?. (9a) 


These arguments for fixing 8 seem rather vague; 
fortunately it turns out that the results obtained in 
ordinary applications are quite insensitive to its exact 
value. 

Errors of input components up to w, are now in- 
cluded in «, errors of higher components in €; €, even- 
tually increases while € decreases with increasing Wm, 
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yet a shift in w,, corresponds to nothing but a changed 
point of view and the computing error must not be 
affected. As €to4 is to provide a measure of this error, 
only that value where it is insensitive to smal] changes 
in W, can be of physical interest; this happens at 
détor/dWm=0; 1.e., at the minimum of €or Finally, 
this minimum will be defined as the dynamic accuracy 
(DA) of the computing element, 


DA B Etot min- (10) 


These considerations also settle the question as to 
when to use (2b) instead of (2a) in the calculation of DA. 
This will be permissible whenever 6?<e? at the fre- Fig. 3—Integrator error functions and dynamic accuracy. 
quency which minimizes €.,. An example is furnished 
by an integrator having the error curves sketched in 
Fig. 3; here the minimum of €.4 corresponds to w= ~, 
where € =0, and DA follows from (2b), (8), (9), and (10). 
The case of an adder is sketched in Fig. 4; calculation 
of DA now requires the use of (2a). 


APPLICATIONS 
The Integrator 


Consider an integrator consisting of a dc amplifier 
of gain —A; and time constant 7, input resistance R 
and feedback capacitance C (Fig. 7). Taking into ac- 
count two time-constants, 7;=7,/A; and T2=A;RC, 
we have (Appendix III) 


Fig. 4—Adder error functions and dynamic accuracy. 


oi 1 1 1 
Hy « = | — i (11) 
Ga Op 0) <0 er pata po 2 e- DT 


where 1 


De ITA EAR Be Dal, OSS bane (la) 


Note that 7) is related to the gain-bandwidth product 
(GB); of the base amplifier through 


(GB); = A;f2rT aS 1/277). (11b) 


Because constant factors in the transfer functions -at 1 v 
~are of no consequence in the sequel they will be freely € 
omitted. The inverse transform is the impulse response Fig. aaa au ECHO a ales imperfect 
_Wwy(r); we use the truncated, r=0---1, part of it Rig Ee ESCH NAIA CMR OUTS NEUE sh 
only and obtain (Fig. 5), ; 
The Integrator with Fixed Input and Feedback Elements 
cn(r) = (€°% — €-*)1(7) 101 — 7). (12) 


Proper external gain is here obtained by adjustment 

Next we calculate DA;. The situation of Fig. 3 holds of an input potentiometer, Fig. 6. To obtain the opti- 

here; 7.e., DA;=«”. For small errors, 1«<7T<«TZ>2; as- mum computing time T=T», we solve 0DA,2/dT =0; 
suming also T>>T,/k;, one obtains from (2b), (Appendix je, Ty)3=37\T.? or ab?=3. The solution is 


IV) \2 
(2) i (=) Ty = 1.4477 TY = 37) (=) (15) 
e600? T A; % 

DAS > 12a * 12 (13) ko is defined by (14) with 7 =7 > and does not include 
the gain of the potentiometer. With this value of T 

| k; is the integrator gain as defined by we obtain an optimum DA i 

A Ons 
k; cay RG. (14) Rio 3 be 0.39 : 
DA io = 0.5 “0.724/— = — gee 11D) 
Two cases must now be distinguished. A; T, ~VJ/(GB);T2 
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Fig. 6—Integrator with fixed input and feedback elements and 
adjustable gain. 


O. 


From (13) and (16) follows 


iE 3 
ne 
(—) peg To 
WA ee 


Yi 
Gy) 
To 
This expression falls off fairly slowly to its minimum 
at 7, but rises rather steeply for 7>T7 . For example, 
at T=T,/20, DA;~3DAi0, whereas at T=207», 
DA;~10DA,. Computing times much larger than 7’) 
should therefore be avoided, but considerably smaller 
times may be used without serious detriment. 


(16a) 


The Integrator of Fixed Gain 


We now demand constant external gain and achieve 
this condition by appropriate changes in the RC-prod- 
uct, Fig. 7. Thus 


Fig. 7—Integrator with fixed gain and adjustable input resistance. 


It follows from (13) that DA; decreases with increasing 
T/T,; it does so, however, more and more slowly. It is 
hardly worthwhile to increase T beyond 


AN? 0.95 /A;\? TUTE 
Tn OL; ) = ( ) = (/ ’ 
k; (GB); \ k; (GB); 


which corresponds to a DA; of 140 per cent of its limit 
for infinite 7. The associated DA; is 


(17) 


oT 0.39 


k; 
DA im 0.40 — = 0.74 = Se (18) 
A; T> ~/(GB;) T» 
Moreover 
TENS SA Soe cand as 
Le ase ay (18a) 
DAim Ue 


The value of RC which is required for operation with 
=T 718 


June 


Pesca gri 
Bree ies 


PRO == (19) 

In practice, the available RC-product is limited to 
some value (RC)max. The relations derived in this sec- 
tion hold, therefore, for k;=k’ only, where 


of see ee 
CG rome (RC) max 
T,,’ denotes the value of 7, which corresponds to R’. 
Smaller values of external gain can only be realized 
by operation according to the scheme of the previous 
section. It will then be advantageous to choose ki=k’ 
and to obtain the balance of gain reduction by adjust- 
ment of the potentiometer. In order to obtain the ap- 
propriate dynamic accuracy DA,’’, we revert to (13) 
and obtain at T=k;(RC) max 


(20) 


3 k; 2 
OT; + (=) RARG) tee 


a 


DA; = 
1 » RIE rae 


: re Py 
~ pe OR bly 
2k:( RC) max 7 


provided k;<~3k’. 

For fixed GB, the gain of the operational amplifier of 
an integrator should be made large, provided one operates 
with optimum computing time. If operation over a wide 
range of computing times J is required, the gain should 
be sufficient for good performance at the longest times; 
for shorter times, one may have to work with T<T» 
and, from (13), DA becomes a function of T and GB 
only; the gain can then be reduced without impairment 
of performance until T is about equal to the new value 
Oredee 

We shall regard the integrator with fixed gain as 
typical throughout the rest of this paper. 

To review: 


(21) 


1) An integrator has an optimum computing time; 
this optimum is not critical. 

2) Amplifier gain determines the value of optimum 
DA; the amplifier gain bandwidth product deter- 
mines the available range of computing times. 


A Numerical Example 
Given: 


A =82 << 104: Lg-= 5) 21072 sec; (RO age = el Seer 


It follows that 


R= 6X2 1025 10 = 18. 


1) For k:=>k', DAn=20/2.5X2X10!=0.04 per 
cent at T=T,=6X5 X10-?X2 X104/202=15 sec, and 
RC=15/20=0.75 sec <1 sec. 
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2) 8 For ep 5'< Be 
DA F=DAS = V/5X10-47/2 5X2 10+ 


=0.05 per cent at T=5 sec. 


The Isolated Adder 


It turns out that an adder cannot function with a 
finite error unless the bandwidth of the input signal is 
limited. In this section we shall consider the perform- 
ance of an isolated adder and in the next that of an 
adder in cascade with an integrator. 

We approximate the transfer function of the adder 
by (Appendix IIT) 


ES 2 


Het fa 
27(GB) a 


1 i Ge apes 
we ’ 3 ae. Om 
p Tr, 


T., is again the time constant of the base amplifer, 
—A, its gain, and —k,= —R,/R; denotes the external 
gain of the adder. Thus 


where c¢ = 


: (22) 
T: 


Moreover C;=1. 
If the input band be limited to 0 - - 
for €? (2a) 


Wm, we obtain 


wm wm 2 
i | Cy [2dw — | Re if Craw 
0 0 
[ "les 
0 


The integrations are straightforward and result in 


Wm 
arc tan? {| —— 
Cc 


Wm Wm 
= NRE ea, | b= 
C C 


Ie Ace if 
1.€., Wm & == : 
ToT - Rha La 


1 (=) 
a s ee Wm a a 
3 if 


Equivalently, we may say that the computing band- 
width B,=Wm/27 is given by 


€ = 


2dw 


ie 


and, provided 


Wm/o K 1; 


(23) 


1.76 
(GB) aT. 


23 
Ae he ee 


fe 
Bm = 0.28 «= 
T3 1 
The computing bandwidth is thus proportional to 
the GB product of the amplifier and to the computing 
time. Note that the frequencies and the bandwidth B,, 
are referred to the numerical time-scale r. 
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We shall now determine the dynamic accuracy of 
the adder. For the total error (8) we have, using (9a) 
for €:, 


0.40 


2 
Wm 


TENE 
Etot? = €17 + €2? = 0.58 =) Wm + 


Now the minimum of f(x) =ax+0/x is 2\/abat x =~+/a/b- 
The dynamic accuracy is, therefore, 


T3 ae T3 
DAa« = €totmman ae 2-0.58 — /0.40 = (Ofek a 
Ife IE 
1+ ke 
221, 3a Ay (24) 
(GB) aT 
at a bandwidth of 
fe O40. (1 Wes a 
By = Wy,/ 24 ar ( ) = 0.166 
Qa 0.582 T3 Ts; 
(GB) .T 
Mssiy vay iabtenaas y (24a) 
1+ k, 


We note that DA, decreases with increasing comput- 
ing time 7. The optimum computing time of a differen- 
tial analyzer is therefore given by T,,, (17), the optimum 
computing time of the integrator. Replacing T by T,, 
in (24) we have 


GB); ki/1 + ka 
DAS — 0.364/ © ) = Ne (25) 
(GB)a A; 
Comparison with DA;,, (18), yields 
DAD 1+ ka)(GB); 
$ Ce 
Darn (GB). 


For best results, (GB), should, therefore, be larger than 
(GB);. For Ramax =3, for example, (GB), =3(GB),; equal- 
izes the accuracies of integrators and adders at optimum 
integrator computing times. 

Note that it is the GB product alone that matters 
for adders. For integrators on the other hand we have 
found that the gain itself should be high. A minimum 
gain exists only if the adder is to be directly calibrated: 
a fractional devation ¢€ of its output from —ke, requires 


tk k 


€ 
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We regard an adder as isolated if it connects two 
branch points of a computing setup. A branch point 
is here to be understood as a point at which three or 
more computing elements connect. Frequently, how- 
ever, adders are cascaded with integrators. Let us in- 
vestigate this important case in some detail. 
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Integrator and Adder in Cascade 


The combined transfer function is now 


1 ee 


r~ - 


(taal tn) 
US rae f D2 


where powers of ~ higher than the second have been 
neglected in the denominator. We have 


1 = T./A;; Shy S AAG = TA,/k; with k; = T/RE: 
Ty (1 er Ai 


The subscripts 7 and a refer to integrator and adder, 
respectively, and T;, —A; and 7,, —A, are the ampli- 
filer time-constants and gains. These results show that 
T, and 73 play similar roles. The optimum computing 
time of the combination is greater than that of the 
isolated integrator, 7, (17); this follows at once from 
(13). The best time for a complete setup will, however, 
be taken as Tn, since there is likely to be at least one 
isolated integrator present and we have seen already 
that 7 must not be larger than 7,,, whereas somewhat 
smaller times are of no consequence. Choosing T= T7\, 
and replacing 7; in (13) by 7:+73 we obtain from 
(13), (17), and (18) after some algebra, 


DA? (=) [1 pees | 
6xtd: De OMGBY: 


2 (GB)a 


Since k=; is fixed, best results are achieved if (GB), 
is several times as large as (GB),; this allows us to work 
with high & and corresponding low k;. Note also the 
possibility of obtaining high accuracy in integrators of 
comparatively low amplifier gain. Cascading an inte- 
grator with an adder allows us to reduce k;; high ac- 
curacy results, provided (GB), is sufficiently large. 


I 


(26) 
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ELEMENTS WITH RC-COUPLED AMPLIFIERS 


As a more complicated application consider a com- 
puting element based on an amplifier with capacitive 
interstage coupling such as is used occasionally in 
repetitive computers. We shall neglect the effect of 
grid current and thus partially invalidate our results. 
In an actual example grid current would have to be 
considered separately. 

Let 7. denote the time-constant introduced into 
the amplifier response by the coupling networks. It is 
shown in Appendix III that integrator performance 
is not substantially affected by T<T,.. For an adder 
(Appendix ITI) 
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= p 1/T; 
ins eee =: — x 
ya ec 
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mee (27) 
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Therefore, (Fig. 8) 
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Fig. 8—lIllustrating the memory function of the RC-coupled 
adder, 


In addition to the ordinary error of adders, caused by 
the nonvanishing rise-time T, of the amplifier, we 
notice an approximate integration due to the finite 
time-constant 7, of its interstage coupling networks. 
Under suitable restrictive assumptions, (Appendix V) 


1 TieN2 Tin / aN 
eateQ+2Q 
3 it Din NIE 
The first term will be recognized as the squared error of 
the direct-coupled adder, (23). The second term ex- 
presses the blocking of direct current. We now add «? 
and €”, (9a), and obtain the square of the total error 
€*tot- 

How does the additional term in &? affect the mini- 
mum of €,o:? with respect to w,,? It increases it and shifts 
it to the right. In order to obtain approximate results, 
we can, therefore, neglect &”?. The minimum of e,? with 


respect to Wm then gives 
es 4/ ise 
A ilies 


As the results for ordinary integrators remain valid we 
now have all the fundamental relations required for the 


(29) 


DA = €imin = 1.6 


(30) 
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design of RC-coupled computing elements. Further 
consideration of this problem would lead us too far 
astray from the main objects of our paper. 


SINGLE-FREQUENCY PERFORMANCE 


We shall now investigate the performance of inte- 
grators and of adders with steady state sine inputs. 
It will be shown that the corresponding errors are usual- 
ly comparatively small so that the use of wide-band 
considerations alone is then justified for our purposes. 
The fractional error was given in (3) as the phase differ- 
ence between AH; and Hy, the ideal and actual transfer 
functions of the computing element. 

For an integrator with transfer function Ay of (11) 


= w? — ab 
< Hy & arc cot ———— a> 0b; 
aw 
tH, = — 1/2 
and, therefore 
T? 
we £2 ee ee 
w? = ab IL SPS 
es Ag & = (31) 
aw wf 
—w 
T, 


Subscript s refers to single frequency performance. 
This error is zero for 
wo = T/V/TiT». 


w= 


(32) 


At the optimum computing time of an integrator T= T°", 


(17), and 


Wo = Wom = 2.45>/A;/ Ri. (32a) 


Let us find the frequency range for which DA, <DAn, 
maintaining the condition 7 =T,,. For high frequencies, 
W>>Woms 


For low frequencies, wKWom, 


-a() 
Ce Ww. 
A; 


Comparison with DA,, of (18) yields the result 


e; > DAm, provided 2.5 < w2.5(A/h). (33) 


The lower limit cannot worry us; operation at fre- 
quencies below w=2.5 is better evaluated by our wide- 
band methods. The upper limit is high enough to be 
quite safe. 


Horsadders ate —7),,, using: (3), (17);-and (22) 
Go 2 w<T,/T. (34) 


Comparison with DA,, (25), results in 
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a (AS Wiccie (oo Nr ees) 
=> ‘a WwW a —— Cs B = me 
DAa (GB)s A; 
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For all reasonable values of w, ¢,< DA, and it follows 
again that wide-band considerations furnish better 
criteria of performance. 

There are, however, exceptions. For integrators, for 
example, it can be desirable to work at very high values 
of external gain k; when inputs with a strongly dominat- 
ing single frequency are to be integrated. This is due to 
the fact that integration of a steady state sine input 
of angular frequency w (in tT) and amplitude gain k’ 
requires k;=wk’. The upper bound in the provision of 
(33) then becomes w<1.6./4;/k’; if this restriction be 
overstepped, e; should be considered. 


DRIFT 


The consideration of drift is outside the scope of 
this paper. Drift will of course influence the choice of 
design parameters. We shall confine ourselves to one 
remark. Let e;g be the input offset voltage required to 
compensate for output drift of an integrator. The out- 
put voltage due to drift, at the end of the computing 
interval T, is then Te,;2, where e;q¢ is assumed to be 
essentially constant. It is well to be aware of the fact 
that in integrators it is this quantity that matters 
rather than e;q itself. 


CONCLUSION 


We have developed a formalism sufficient for a system- 
atic evaluation of the performance of linear computing 
elements using the effective error and some new concepts 
such as the computing bandwidth and the computing 
transfer function and introducing a quantitative defini- 
tion of dynamic accuracy. Practical applications always 
deal with small errors and actual calculations can, 
therefore, be considerably simplified through suitable 
approximations. 

Although dynamic accuracies have been calculated, 
only a few conclusions were drawn in this paper about 
the choice of parameters of design. The important ques- 
tion of the requirements for the best elements within 
the limitations of existing components were not touched; 
it is hoped to cover this topic in a further paper. An- 
other question of interest would be the establishment 
of experimental techniques for the determination of DA 
and of computing bandwidth as here defined. 
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THE EFFECT OF CHANGE OF SCALE ON OUTPUT 
AND ON TRANSFER FUNCTION 


On the new time scale 7, output is given in terms of 
input and impulse response by 


ear) = f HOW oe ade 
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and on the old scale ¢ 


Tr 
[ eulx)W (Tr — x)dx 
Lt Nit x 
= HE es( ) — w(: — =) dx 
0 Sad id Te 


= [eds = y)dy = els) 
0 


=Tr by 


€oi(t) — Corl Tr) = 


which proves the first relation. For the transfer func- 


tion, 
Hi(p) = if W (t)e—?'dt 
0 
and 


H(p) =k Woerdr = f TW (Tr)e—?*dr 
0 0 
2 { W(t Ola, = Fl, (=) 
0 th 
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as stated. 


PROOF OF THE RELATION BETWEEN ORDINARY AND 
COMPUTING TRANSFER FUNCTIONS 


Transforming (4) we obtain 


[wera is i a iii 


= H — W (r)e~?"dr 


1 


C 


= Hi — of W(r + 1)e-?"dr 
0 


=H --(WGF DIF), 
which is (5). 
APPENDIX III 
DERIVATION OF TRANSFER FUNCTIONS 


Consider the computing element of Fig. 9. Z; and 
Z; are the input and feedback impedances; — A H,(p) is 
the transfer function of the amplifier, where H/,(0) =1 
and A> 1. 

The transfer function of the element is 


Ze 1 


Hay 2 ee 
AH,(p) Z; 


HEC fo (36) 


The Integrator 
We take for the amplifier transfer function 
1 
19 0 Ve 
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Fig. 9—Block diagram of a computing element. 


Moreover 
wee 1 IE 
2; = R;: Ly = —; pe 
pc RC 
and thus 
4 k; 1 
Hi(p) = — 7 ; (37) 
(p+ = Jere + 1) 
T» 
where the time constants are 
LT A; 
a > Pa Re Se (37a) 
The Adder 
With H, as in above and 
wed Ea R 
Ze Ree) 2 ee 
R 
we obtain 
H i (38) 
tient ’ 
1 cy 
with a time constant 
1+ &, 
ilies = y hath (38a) 


Now 


but 


H,= 1 5 (hod be 
Tigers a) = 
(pT. + (e+ =) 


T. is the large time constant of the capacitive inter- 
stage coupling networks of the amplifier. It follows 


1 
H,=-— : P ) 
RC * 5H oe 1 Tey 1 
e : Den Ang RC Wig 


provided A;RC > T.. 


(39) 
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For high frequencies; 7.e., small ft, we neglect the last 
term in the denomimator. The result is then identical 
with that for the direct coupled integrator, (37). For 
low frequencies, or large t, we neglect the first term in 
the denominator and 


Bre A ie PL p Mite Wes 
Des Noes UNAS eer 
RE ( ee 1 ) A;RCT, 
PY SARC) + wp? 
> 0) 
Thus 
Wi)2=- ine eVARC cos wt l(t) 


1 t 2 8 

2 ( = = )10, t << /A;RCT; 
RC OARO s 2UROT: 

£>>0. 


Combining these results we obtain without restriction 


1 t ie 
ie zl — H/T, — a | 1(t); 
RC QA:RC . QA-ROCT, 
t << /SA;RCT.. (39a) 


It follows that the effect of the coupling term is 
negligible as long as the computing time T is smaller 
than T,, the coupling time constant. 


The RC-Coupled Adder 
Z, = Ry; 


as for the direct-coupled adder and dH, as in the last 
Section. 


Hip) « — 
2) aa [ Aon 1 le 1 
P ? 1+k, T, is Ihe ead hy 
p 
SY 40 
GUS) 
p T3 p IP) 
where 
1+ k, Ala 
te a) 1Bh = lines 
As 1+ k, 
APPENDIX IV 
DETERMINATION OF INTEGRATOR DA 
For an integrator we can write (2b), using 


Cr=1(r)1(1—7), 
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ee) 
w= pian (cw Sey 
cn? 


where the wriggles indicate averages. Now a>>b and 
we shall also assume ab>>1; 7.e., 


A,T k; Z 
—>1 or TSO—- 
bp! A; k; 
Therefore 
~ exe are 1 e? ae 1 
CN = (e-* — €&")dr = — 
0 a b 
a b 
~i{-—-—-—; 
a 2 
cy? & 1—-——621; 
a 
——~ 1 ; 
(cy — cy)? = (a —-e™—1+4+—+4 3) dr 
0 a 
1 1 1 
= —+ — b+ = ~>—+ 6°, 
Del 2 a a? 2a 12 
and 
6+ ab? 
DA? = 
12a 


which is (13). 
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CALCULATION OF ACCURACY OF THE RC-COUPLED 
ADDER 


We shall use a short cut in order to avoid the lengthy 
direct calculation. Eq. (28) consists of two terms. The 
first will be recognized as the impulse response of the 
direct coupled adder. The second term is approxi- 
mately equal to (1/74)(r)1; it represents, therefore, the 
impulse response of integration. Neglecting the small 
interaction between both terms we take over (23) 
for the error caused by the first term; 7.e., 


The second term causes nothing but an error output; 
it corresponds to a computing transfer function 


or 


icp = (- ya ; 
= —(— — cos w). 
| w? T4 i 
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Our standard input is given by 
é;| = l(w)l(w — wr). 


The normalized output is the second error term and 
we have 


1 is a 2 {1 \* fm 1 — cos w 
19° = -{ WB ART Se ( ) [ — dw 
Wmv 0 . Wm ili J 90 w? 
rs (- ) 
os ’ 

WeorNed 4 


Adding both terms, 
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The integral can be calculated by contour integration. 
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provided w,, > 1. 
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The salient steps follow. 
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Trigonometric Resolution in Analog Computers 
by Means of Multiplier Elements* 


R. M. HOWE7 anv E. G. GILBERT 


Summary—A method of generating sine and cosine functions 
in analog computers by means of multiplier elements and integrators 
is discussed. Static accuracy of the method is analyzed and found to 
be essentially equal to the accuracy of the multipliers employed. 
The system accepts d6/dt as the input and generates output voltages 
of sin 6 and cos 6. Amplitude-stabilizing loops are employed to main- 
tain sin? ¢+-cos? 6=1. Advantages of the method include representa- 
tion of unlimited range in angles, dynamic capabilities far beyond that 
of the multipliers alone, and possibility of employing electronic multi- 
pliers. The method has been successfully used to compute Euler 
angles in analog computer solutions of the three-dimensional flight 
equations. 


INTRODUCTION 


ANY of the problems solved by analog com- 
NV puters require the generation of trigonometric 

functions. For example, in solving the three- 
dimensional flight equations, the Euler-angle computa- 
tion requires generation of sine and cosine functions of 
the Euler angles. In present general-purpose electronic 
differential analyzers this is usually accomplished by 
servo-driven sine-cosine potentiometers or by diode 
function generators. 


* Manuscript received by the PGEC, June 25, 1956; revised 
manuscript received, February 14, 1957. This research was sup- 
ported by the Dept. of the Air Force, Wright Air Development Cen- 
ter, Dayton, Ohio, under Contract No. AF 33(616)-2839. 

+ Dept. of Aero. Eng., University of Michigan, Ann Arbor, Mich. 


The purpose of this paper is to describe an alterna- 
tive method of generating sine and cosine functions 
which involves only multipliers as nonlinear elements. 
Furthermore, this new method allows continuous angu- 
lar travel through an indefinite number of revolutions, 
whether servo or all-electronic multipliers are used. 

This resolving method has been discovered independ- 
ently by other groups, but to our knowledge has not 
previously been published nor analyzed completely with 
regard to static and dynamic accuracy.! 


DESCRIPTION OF THE NEW RESOLVING TECHNIQUE 


The new resolving technique is based on the solu- 
tion of the equation 


ee (1) 


' Similar self-correcting methods are described for stabilization of 
harmonic oscillators and for direction-cosine computation in the fol- 
lowing notes issued by the University of Chicago, Advisory Board on 
Simulation: 
M. D. Clement, “A Stabilized Harmonic Oscillator,” IN-21, 
December 10, 1954. 

J. P. Corbett, “Orthogonalizing and Normalizing Corrections for 
Direction Cosine Computation by Differential Analyzers,” 
TN 54, December 2, 1954. 


This method has also been used by the Analog Computer Section, 
Aero. Res. Labs., Wright Air Dev. Ctr. 
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with initial conditions 
“yan (2 
dase as 


The solution to (1) with the conditions in (2) is simply 


“= sin 0 (3) 
<= cos 6 (4) 
= — = cos 
é dé 


which can then be used to provide coordinate resolu- 
tions involving the angle 6. 

In @) the independent variable is 6, whereas the 
electronic differential analyzer (electronic analog com- 
puter) can integrate only with respect to time ¢. But 
by multiplying the integrand by d@/dt before integrating, 
the integration with respect to @ is effectively accom- 


plished. Thus 
do 
f-G) Whe fa. (5) 
dt 


The electronic differential analyzer circuit for accom- 
plishing the solution to (1) is shown in Fig. 1. Note that 
the input to each integrator is multiplied by d0/dt 
(i.e., 6) before being integrated with respect to time. 
The multipliers shown schematically in the circuit could 
be either servo or all-electronic. 


xd@/dt 


One Unit int=RC Seconds 


Fig. 1—Simplified computer circuit for trigonometric resolution. 


One difficulty with the circuit as it is shown in Fig. 
1 is the tendency of the solution to drift so that cos? 
6+sin2 0641. For example, if 6=w=constant, the theo- 
retical solution would be x=sin wt, y=cos wt. Due to 
power loss in the feedback capacitors the solution may 
damp somewhat. The net result is that after several 
cycles x?+y? is no longer unity, but perhaps a per cent 
too low. Other computer errors such as dynamic multi- 
plier errors, integrator drifts, etc., can cause similar 
errors. 

This error in amplitude can be calculated using multi- 
pliers and a summer to obtain 1—(x?+y’). Thus define 
the error € by 


a aas) (6) 


The magnitude of « and y can be adjusted automatically 
upward or downward as needed by providing each inte- 
grator with negative or positive damping proportional 
to e. This is shown in Fig. 2 where ¢€ is computed by (6) 
above, multiplied by x and y, and fed back into inte- 
grators 2 and 1, respectively. 
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Fig. 2—Complete computer circuit for trigonometric resolution. 


When servo-driven potentiometers are used for multi- 
plication, this resolving circuit requires two servos, 
each with at least 3 multiplying potentiometers in 
addition to the reference potentiometer. The shaft 
angle of each servo is proportional to cos @ and sin 86, 
respectively. This means that additional ganged po- 
tentiometers on each shaft can be used to multiply cos 
6 or sin 6 by other voltages. The output voltages of 
integrators 1 and 2 are also proportional to cos @ and 
sin 9, respectively. 

If all-electronic multipliers of the time-division type 
are used, two masters (one for cos 6, one for sin @) are 
needed, with 3 slaves for each master. Additional slaves 
can be used to multiply cos @ or sin 6 by other voltages. 

To summarize, the resolving method described here 
accepts +6 and —@ voltages as inputs, providing output 
voltages and shaft angles (or master pulses) proportion- 
al to cos 8 and sin 6. The over-all requirement is two 
multipliers, each with three channels, two integrating 
amplifiers, and four summing amplifiers. 

In the next section the static and dynamic accuracy 
of this resolving technique is analyzed. It will be shown 
that servo multiplers are to be preferred over electronic 
multipliers if the servo bandwidth is adequate. 


STATIC AND Dynamic AccURACY ANALYSIS OF THE 
RESOLVING TECHNIQUE 


In this section a detailed study is made of the factors 
influencing resolver accuracy. It is shown that static ac- 
curacy is as good as the multiplier static accuracy, while 
dynamic accuracy is generally better than the multiplier 
dynamic accuracy. Exact solution of the resolver equa- 
tion is important and will be developed first. 


Exact Solution of the Resolver Equations 


It is convenient to write the second order resolver 
equation solved in Fig. 2 as two first order equations. 
Letting the first and second integrator outputs be y and 
x, respectively, one obtains 


alse johe 


time SH Ave x2 eens ie " 
Tale ee pa? +4 ve (7) 


~ 
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and 


d y dé < : 
—_= — a — p(x? + y? — 1)y. (8) 
dl dl 
The analytic solution for these equations is best ob- 
tained by expressing x and y in polar form. Thus define 
the polar variables R and ¢ by 


R sin d (9) 
R cos ¢. (10) 


I 


x 


y 
In terms of R and ¢ perfect resolution is obtained when 
R=1 and 6=¢. Substituting (9) and (10) into (7) and 
(8) gives 
do dR . 
R—cos ¢ + — sin ¢ 
dt dt 


d6 
= nee —u(R?—1)Rsing (11) 
t 


wage aR i 
— se SS roa) 
Bare eT ag 


dé 
a Rsin ¢ — u(R? — 1)Rcos¢. (12) 


Eqs. (11) and (12) can be written as single equations 
in R and ¢ in the following way. Multiply (11) by sin @ 
and (12) by cos ¢ and add, obtaining 
dR 
dt 
Similarly, multiply (11) by cos ¢ and (12) by sin ¢ and 
subtract, getting 


(13) 


—, or ¢=6+4+ ¢(0). (14) 
Eq. (14) implies that if (11) and (12) are solved exactly 
with ideal computer components, as indicated in Fig. 2, 
then ¢=6 as required, unless the initial value of @ is in 
error by (0). But this possibility of error is present also 
in a conventional resolver if d0/dt is integrated directly 
to obtain @. 

To study the behavior of (13) it is convenient to 
introduce the variable 7, which is the deviation of R 
from unity. Thus let 


r= R—-1. (15) 
Eq. (13) then becomes 
dr 
5a BA et 2 RAT (16) 
In the usual case v1 and (16) reduces to 
a = — 2uyr, AGA: (17) 
dt 
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For an initial error in x?+y?, 7.e., for R(O) #1, 7 has an 
initial value 7o, and the solution to (17) is simply 


(18) 


7 = roe 7H. 


Thus the error in R decays exponentially to zero with 
a time constant 1/(2u). This time constant can be re- 
duced by making the correction gain-factor pw large. 
Up to now, then, we have shown that accurate resolu- 
tion is obtained when ideal computing elements are 
used, despite any integrator initial-condition errors. 
At worst the error in angle is simply a constant (0) 
equal to the initial angular error caused by incorrect 
integrator initial conditions, while the error 7 in magni- 
tude decays rapidly to zero. In the following sections 
the errors in resolving due to both static and dynamic 
errors in the computing components will be analyzed. 


Errors Due to Static Inaccuracies in Components 


The static errors in computing components take the 
form of zero offsets, nonlinearities, and gain inaccuracies. 
To implement the following analysis, such com- 
ponent errors must be defined in mathematical terms. 
For a multiplier whose inputs are x and y one may write 


multiplier output = xy + E(x, y) (19) 


where xy is the desired product and E(x, y) an error 
term. If a servomultiplier is used and x is the shaft 
quantity, y does not contribute directly to the error 
and 


E(x, y) = yP(x). (20) 


P(x) is an error function dependent on servoshaft posi- 
tion error and potentiometer nonlinearity. The servo- 
multiplier possesses, therefore, a zero point error ad- 
vantage in that E(x, 0) =0. Integrator errors, primarily 
due to amplifier unbalance and grid current, may be 
represented by an equivalent voltage IJ which appears 
at the integrator input. Integrator gain errors, sum- 
ming amplifier gain errors, and summing amplifier off- 
sets are relatively minor and may be included in the 
error analysis of Fig. 2 as part of the multiplier errors. 

The component errors as just described may now be 
introduced into (7) and (8). Subscripts on the quantities 
E, P, and I refer to the component numbers in Fig. 2. 
After making the necessary substitutions and trans- 
formations (see Appendix), one obtains 


R—12r= — 3[E;(cos 0, cos 6) + E,(sin 0, sin@)] (21) 
; : ‘ 1 
¢=6 — cosé | Ba(cos 6, —6) + es 6,0) + r| 
. . i 1 
+ sin @ | za(sin 6,6) + ig E;(cos 6, 0) + 1| er a4) 


These equations give the minimum static errors and 
hold if » is made large by the inequality 
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‘ 1 
| Aa + ma Eo + I» 


1 
+ | BtoB th 
> - 


23 
| Es + E,| ae 


max 


and @ varies slowly. Hence, the errors H3; and £, in 
multipliers 3 and 4 affect only 7 and not ¢, while the 
errors /;, E:, Es, and Es in multipliers 1, 2, 5, and 6 
affect only @ and not r. Note that the influence of Es 
and £5 can be made smaller by making k as large as pos- 
sible without causing amplifiers 3 and 4 to saturate. If 
the integrator errors J; and J, are negligible, a good as- 
sumption, it then becomes clear that the error in angu- 
lar rate ¢ of the resolver is the order of the errors £,; and 
Fz, in multipliers 1 and 2. For many electronic multi- 
pliers this error may be the order of 0.1 per cent of full 
scale, even when 6=0. A drift rate in ¢ of 0.1 per cent 
of maximum ¢ would result, and this could be serious in 
certain applications. 

On the other hand, if servomultipliers are used, we 
see from (20) that @ in (22) becomes 


_ @ = 6[1 + Pi(cos 6) cos 6+ P2(sin 8) sin 0] 


— I, cos 6+ 7, sin @. (24) 


Multiplier errors P; and P: are simply equivalent to 
errors in integrator gains, and when 6=0, the drift error 
in ¢ is the order of integrator drifts, which are normally 
quite small. 

In summary it has been shown that static errors in the 
resolver outputs are minimized when uw is large as deter- 
mined by (23) and servomultipliers are employed. Under 
these conditions the following statements apply: 


1) When the outputs «=sin 6 and y=cos @ are ex- 
pressed with polar coordinates R and ¢, where 
ideally R=1 and ¢=8, then the error 7 in the value 
of R is given by (21) and is approximately equal to 
the average error in multipliers 3 and 4. 

Errors in the integration of 6 are the same order as 
errors in multipliers 1 and 2. 

Drifts in ¢ are at least as small as integrator drifts. 
Errors in multipliers 5 and 6 are relatively unim- 
portant. 


2) 


3) 
4) 


Errors due to Potentiometer Granularity in Servomulti- 
pliers 


Another type of static error is created by the step dis- 
continuities in wire-wound potentiometers and is char- 
acteristic of both new and conventional resolving tech- 
niques which employ servo-driven potentiometers. As- 
sume that one-turn potentiometers with a total of NV 
turns of wire are used. In the conventional sine-cosine 
potentiometer method of resolving, a full revolution rep- 
resents 27 radians, so that there are NV /27 turns of wire 
per radian. In the new resolving method a full multiplier 
potentiometer revolution represents a range in sin @ or 
cos §@ from —1 to +1, 7.e., a total range of 2. Thus in the 
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region where @=0 or z (the most critical regions for po- 
tentiometer granularity in computing sin @) sin 0&6, 
and there are N/2 turns of wire per radian. Hence the 
effect of wire-wound potentiometer granularity is re- 
duced by a factor of t when comparing new and con- 
ventional resolving techniques. This feature is im- 
portant in obtaining the maximum smoothness of re- 
solver output. 


Dynamic Errors 


In this section dynamic errors of the new resolving 
method are analyzed by considering sinusoidal response 
characteristics. Two general cases will be considered: 
1) resolver performance when d@/dt=constant, so that 
ideally the resolver outputs are sin wt and cos wt, where 
w =6; 2) resolver performance when 6=w6 cos wt, where 
61. Only servomultipliers will be considered. 

First let us analyze the case where 6=constant. In 
the steady state it is clear that the outputs y and x of 
integrators 1 and 2 will be sinusoidal with constant am- 
plitude. If both multiplier dynamics are the same, it is 
evident from symmetry that x and y will have the same 
steady-state amplitude Rk. For the sinusoidal outputs 
to remain constant it is also clear that the output of in- 
tegrator 1 must lead the output of integrator 2 by 7/2, 
so that total phase shift around the loop is 7. Thus let 


(25) 
(26) 


x = RK sin wt 


R cos wt 


where the constants R and w are to be determined. 
Next let us assume that the output voltages of the 
multipliers can be described in terms of a magnitude MZ 
and a phase-shift NV for sinusoidal inputs of magnitude 
R and frequency w. M and N can either be calculated 
from previous information or can be measured experi- 
mentally. Thus the output of multiplier 1 can be written 


yi 


output of multiplier 1 = — R6M cos (wt+ N). (27) 
In the same way 
output of multiplier 2 = ROM sin (wi+ WN). (28) 


Similarly, 
output of multiplier 3 


R°M cos wt cos (wt + N) (29) 


and 


output of multiplier 4 = R?M sin wt sin (wi + N). (30) 
Thus the output of amplifier 4 is given by 

output of amplifier 4 

= ku{ R?M [sin wtsin (wf + N) + coswtcos (wt + N)] — 1} 
= ku[R?M cos N — 1]. (31) 
This signal is fed through multiplier 5 to give 

output of multiplier 5 


= kuMR(R°M cos N — 1) cos (wt + WN) (32) 
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which, finally, is multiplied by 1/k and summed with 
the output of multiplier 2 to give the input to integrator 
1. But from (26) the input to integrator 1 1s also given 
by wR sin wt. Hence from (28) and (32) we have 

R6M sin (wi + NV) 

+ uMR(R2M cos N — 1) cos (wt + N) = oR sin ot 


or 


MJ + u2(R2M cos N — 1)? sin (wf + N + 7) 


= sin wl 


(33) 
where 
u(R?2M cos N — 1) 


tan y = ; mi 
6 


(34) 


Equating magnitudes and phase angles on both sides of 
(33), we obtain 


1 6 i 
ya ae eee (: — — tan ) (35) 
M cos N Me 


and 


M6 


yo = —— + (36) 

cos V 
For the case where the phase error NV of the multiplier is 
small and where N6/u<1, it follows that 


1 6N \ a 
Se —), N «1, —<«1 (37) 
V/M 2u Me 


and 
(38) 


The effect of the multiplier magnitude factor M and 
phase error N is now clear when 6=constant; the mag- 
nitude of the resulting sinusoidal response is shifted 
slightly from unity in accordance with (37), while the 
angular frequency is in error by the factor M—1. 

It is important to note that FX in (37) is the magnitude 
of the sinusoidal output of integrators 1 and 2. If addi- 
tional trigonometric multiplications are to be performed 
using extra pots on the shafts of servos 1 and 2, then the 
magnitude R, of the sin wt and cos wt factors in the ex- 
ternal trigonometric multiplications will be given by 


6 1/2 
(1 — — tan v) 
MU 


[Re = WOR = M}!?, (39) 
(cos V)1/2 
For small N this reduces to 
ON N6 
IRs = (: -~) M1? N «1,5 <1. (40) 
2u a 


The phase shift will of course be NV radians with respect 
to the integrator 1 and 2 outputs respectively. 
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As an example of the order of magnitude of errors in- 
volved, consider a servomultiplier with an undamped 
natural frequency of 100 radians/second (~16 cps) and 
let 6=3 radians/second (~0.5 cps). Reasonable values 
for M and N under these conditions would be MJ = 1.001 
and N=0.01 radian. Assume that »=10. From (40) the 
error in magnitude at the servoshafts will be 0.1 per cent 
while from (38) the error in frequency will be 0.1 per 
cent. 

Next let us consider the dynamic errors in the resolver 
method when servos are used and the input angular rate 
§ is a sinusoidal oscillation given by 


6 = dw cos wi (41) 

from which 
6 = 0) + 6 sin at. (42) 
We will consider the case where 6<1, 7.e., where the 


angle # consists of a small sinusoidal oscillation about a 
constant angle 6). This type of dynamic performance is 
often required in solving aircraft or missile flight equa- 
tions. 

It is difficult to solve this problem analytically for 
the general case where the multipliers have a magnitude 
M and phase shift N for each of the harmonic fre- 
quencies which will make up the steady-state resolver 
outputs x and y. However, of particular interest is the 
solution where M=0, 7.e., where the frequency w is so 
high that the servomultipliers cannot follow the outputs 
x and y at all. In this case the average output voltages 
x and y from integrators 2 and 1 in Fig. 2 will be sin 6 
and cos 6) respectively, and the shaft angles of multi- 
pliers 1 and 2 will have constant displacements given by 
sin 4) and cos 4, respectively. Thus the input to integra- 
grator 2 will be —6 cos 0) = —6w cos wt cos 80, from which 


the output x is given by 
x% = sin 4) + 6 cos 4 sin wt. (43) 


Similarly, the input to integrator 1 is 6 sin 0) =éw cos wt 
sin 4, so that the output y is given by 


y = CoS 8) — 6 Sin 4% sin wt. 


(44) 


The outputs of a perfect resolver would be given from 
(42) by 


x = sin # = sin (4) + 6 sin wf) (45) 
y = cos 6 = cos (6) + 6 sin wi). (46) 

Eq. (45) can be rewritten as 
*% = sin 8) Cos (6 sin wt) + cos 6 sin (6 sin wi). (47) 


Expanding cos (6 sin wt) and sin (6 sin wf) in power 


series, we obtain 


: (6 sin wt)? 
4 = sin 9 Lee a 
2 


(6 sin wt)? 


4+ 608 ] Asin af — a tee |. (48) 
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It is now clear that the resolver output « given by (43) 
when multiplier dynamic response is zero equals the 
zeroth and first order terms in the series representation 
for the ideal resolver, as given in (48). This means that 
the actual resolver output will be very nearly equal to 
the ideal output providing 61. Similarly, from (46) 


(6 sin wt)? 
y= cosy 1 — S28 | 


(6 sin wt)? 


— sin 4 in wt — 
3! 


co 


which shows that the actual resolver output y given in 
(44) represents the first two terms of the ideal resolver 
output.. 

To summarize, if the angle @ consists of a constant 4 
plus a small-amplitude sinusoidal oscillation, the re- 
solving method gives accurate sin 6 and cos @ outputs 
even when the frequency w is so high that the servo- 
multipliers have no response at all at that frequency. 
Unless rather severe resonant peaks are present in the 
servos, it seems reasonable to assume that at lower fre- 
quencies the results would be somewhat better. Note 
that these conclusions apply to the voltages x and y ap- 
pearing at the integrator outputs. The servoshaft mo- 
tion will not be accurate at all, and electronic multi- 
pliers would need to be used to multiply x and y by 
other voltages, as required in the particular problem. 

This suggests the possibility of employing servo- 
multipliers for the resolving method because of drift 
considerations, but still using electronic multipliers for 
further required multiplications by the resolver outputs 
e—sin 0 and y=cos 0. 

One disadvantage of this resolving technique is the 
requirement of the time derivative d0@/dt of the angle as 
an input voltage. However, in many cases d6/dt is 
readily available because of the problem mechanization. 
For example, this is true in computing the sines and 
cosines of the Euler angles, as needed in a three- 
dimensional flight simulation. 


RECTILINEAR TO POLAR RESOLUTION 


Thus far we have described how to use the new re- 
solving method to produce sin 6 and cos 0, given d6/dt. 
Actually, this is the polar to rectilinear conversion prob- 
lem, 1.e., given 6 and Rk, compute a=R cos 6 and b=R 
sin 6. An equally important problem is to calculate 0 
and R, given a and 8, 7.e., the rectilinear to polar con- 
version problem. This can be accomplished by recogniz- 
ing that 


b sin 6 ' 
or bcos @ — asin #@ = 0. 


a cos 6 


Using the new resolving method 6 cos 6—a sin @ is com- 
puted. This will in general not be zero because @ will not 
have its proper value. The resulting nonzero signal is 


Resolution in Analog Computers 


91 


amplified and fed into the d6/dt inputs to the resolver, 
so that the sin @ and cos 8 outputs are driven to the 
appropriate values representing the correct 6.2 By in- 
tegrating d@/dt at the same time, the voltage @ is ob- 
tained. R is calculated from the formula 


R = acos 0 + bsin 6. 
If R varies over wide limits in any particular problem, it 
will be necessary to use an automatic gain control in the 
d6/dt input to make the loop gain proportional to 1/R. 
This method of rectilinear to polar conversion has 
several disadvantages. The angle @ (a voltage) and the 
functions sin @ and cos @ are self-consistent only if the 
integrator which produces @ from 6 is perfect. As a re- 
sult, integrator drift and gain error limit the time of 
accurate computation. Furthermore, the angular range 
of @ is restricted by amplifier saturation. The more con- 
ventional servo resolver which has a shaft angle 6 does 
not possess these limitations. 


APPENDIX 
STATIC ERROR EQUATIONS 


Introduction of error functions in (7) and (8) gives 


z = by — ula? + 9? — 1)x — Bx(y, —6) 


—n[Es(y, y) + Ea(x, x) |x — — Bal, — hye) — Iz (50) 
y = — bx — p(x? + y? — 1)y — Ea(s, 8) 
—nlEs(y,y) + Es(x,x) ly — - Ey, — kee) dy OD) 
where kué is now the output of amplifier 3 or 
€é= 1 — (x? + y’) — By — Ez. (52) 


These equations are next converted to polar form ac- 
cording to (9) and (10) and the following expressions 
for Rand ¢ are obtained in the same manner as (13) and 
(14) were obtained. 


R = p(R2— 1) R RE eA 
1 
= sing] B+ Be +h 
1 
= Os 6| Bs + Bae 4 1 | (53) 
5 = 6 | a tg +1| 
o— R 1 a ~6 2 
sin 1, sh 
ans R | Bat b Be + 1 | (94) 


2 The amount of amplification in the dé/dt input will determine 
the time constant of the corrective action. 
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. 1 1 
itn | Es + Es| >—|Ei+— Eat Is 
a= 1 — FR? — By — Ey, (55) i 
1 1 
If the resolver errors are small, r&1 and ¢&@ which ay te E» + Fedor Ty}. (59) 
allows (53)-(55) to be simplified to 
Thus a minimum error 7 is obtained when (59) is satis- 
l 1 sin 6 l fied and is given by 
on eed Pag On ee nel B+ Bet be] Soy eigre : 
2u ’ 2 k r = — 4/E;(cos 6, cos 0) + E,(sin 8, sin 6)| (60) 
= £98 aoe Ne a J | (56) provided @ varies slowly enough for */2u to be small 
Qu ae compared to r. Note that large » makes é&0, a condi- 
1 tion which might be expected from the nature of the 
é = 6 — cos of oh ae es 1 | feedback loops. Under the same conditions 
; 2 
1 
1 62Zo- cos | #s(cos 6, —6) + — E,(sin 6,0) + h| 
+ sin o| Ba + 7 Es+ 1| (57) k 
j 1 
ay oe (58) + sin 6| za(sin 6,0) + a ee 6,0) + |. (61) 


Eq. (56) shows that increases in yp will effect decreases in 
r until pw is so large that 


Eqs. (60) and (61) are the ones presented in the section 
on static accuracy. 


Minimization of the Partially-Developed 
Transfer Tree’ 
MIECHBEEIE.P: MARCUST 


Summary—A transfer tree is a particular type of multiterminal 
network having a single input which may be connected to any one of 
a number of outputs. 

An n-relay transfer tree is partially developed if it has less than 
the 2” possible output terminals. Rearrangement of a partially-de- 
veloped tree can lead to a reduction in the total number of transfers 
required. 

This paper presents a method of obtaining a required partially- 
developed transfer tree with the minimum number of transfers. 


INTRODUCTION 


TRANSFER TREE! is a particular type of 
multiterminal relay network having a single in- 


put which may be connected to any one of a 
number of outputs. One and only one output is con- 
nected to the input at a given time, the selection being 
controlled by the combination of network relays op- 


* Manuscript received by the PGEC, July 3, 1956; revised manu- 
script received, January 17, 1957. 

{ Eng. Lab., IBM Corp., Endicott, N. Y. 

1S. H. Washburn, “Relay trees and symmetric circuits,” A7EE 
Trans. vol. 68, pp. 582-586; 1949. 

2 W. Keister, A. E. Ritchie, and S. H. Washburn, “The Design of 
Switching Circuits,” D. Van Nostrand Co., Inc., New York, N. Y.; 
1951. 


erated. Each input-to-output path passes through one 
and only one contact on each relay in the network, and 
all outputs are disjunctive. 

The number of possible relay combinations with n re- 
lays is 2”; therefore, in an n-relay transfer tree there 
are 2” possible outputs. A fully-developed tree has an 
output terminal for each of the 2” possible relay com- 
binations, and the total number of transfers in the tree 
is 2”—1. A partially-developed tree has less than 2” out- 
put terminals, and the total number of transfers in the 
tree may vary. 

By the rearrangement of a fully-developed tree, the 
contact load on all relays except the first may be made 
more uniform; however, regardless of the arrangement, 
the total number of transfers in the tree remains 2”—1. 
Rearrangement of a partially-developed tree, on the 
other hand, can lead to a reduction in the total number 
of transfers required. 

The majority of transfer trees used in switching cir- 
cuits are partially developed, and it is desirable that 
there be a method of minimizing them, 7.e., obtaining 
a required tree with the minimum number of transfers, 
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and recognizing that the solution is a minimal one. All 
possible arrangements of a tree could, of course, be ob- 
tained, and the minimal one selected, but this process 
is long and laborious, prohibitively so in many cases. 
Table I gives some indication of the progressive com- 
plexity of a trial and error approach as m increases. 


TABLE I 


Number of 
relays in tree 


Number of possible arrangements of fully- 
developed tree 


i) jie 2 
3 31592 12 
4 41.32.24 = 576 
5 51.42. 34.28 = 1,658,880 
6 61-52-44. 38.216 = 16,511,297, 126,400 


The number of possible arrangements of an n-relay tree 
is represented by 


[m?*][( — 1)" ][(m — 2)" ][(m — 3)**] - 


If P represents the number of possible arrangements of 
an n-relay tree, then the number of possible arrange- 
ments of an (7+1)-relay tree is (n+1)P?. 

This paper presents a method of minimizing partially- 
developed transfer trees. 


aan 


BASIS FOR METHOD 


In a fully-developed transfer tree, such as the one 
shown for reference in Fig. 1, a branch (minor tree) 
representing 2% combinations contains 2*—1 transfers. 
For example, in Fig. 1, the branch labeled (2) represents 
four combinations and contains three transfers; the 
branch labeled (4) represents one combination and con- 
tains no transfers. 

A particular partially-developed tree can be obtained 
by starting with a fully-developed tree and removing 
the branches representing the groups of invalid output 
combinations. The removal of a branch representing a 
group of 2° invalid combinations results in the elimina- 
tion of 2*—1 transfers. For instance, in Fig. 1, if the four 
combinations A’BCD, A’BCD’, A'BC'D and A’BC'’D’ 
were invalid, the branch labeled (2), representing this 
group of four combinations, could be removed, eliminat- 
ing three transfers. 

It follows that the total number of transfers elimi- 
nated equals the number of invalid combinations minus 
the number of groups into which these combinations are 
combined. The method of minimizing a _ partially- 
developed tree, therefore, consists of starting with a 
hypothetical fully-developed tree and eliminating the 
maximum number of transfers by arranging the order 
of the relays in the tree such that the invalid combina- 
tions can be combined into the minimum number of 
groups. 

The key to the method, then, lies not in the analysis 
of the valid relay combinations, but rather in the 
analysis of the invalid combinations. 
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(2) 


Branch (1) represents 8 combinations and contains 7 transfers. 

Branch (2) represents 4 combinations and contains 3 transfers. 

Branch (3) represents 2 combinations and contains 1 transfer. 

Branch (4) represents 1 combination and contains 0 transfers. 

Removal of a branch representing 27 combinations results in the elim- 
ination of 2*—1 transfers. 


Fig. 1—Fully-developed transfer tree. 


METHOD 


In this paper, maps (charts)*4 are used as the means 
for combining the invalid combinations and obtaining 
the optimum order of the relays in the tree. There are 
notable differences, however, between the normal use of 
the maps and their use here; in this method 1) the in- 
valid rather than the valid combinations are of prime 
consideration, 2) each combination is not considered as 
often as desired but is considered only once, and 3) the 
groups formed must be compatible with fundamental 
transfer tree configuration. 

The use of the map to obtain any desired partially- 
developed n-relay tree (not necessarily a minimal one) 
will be described first. 

An n-variable map is drawn, and an X is entered in 
each square representing a valid combination, as illus- 
trated in Fig. 2 (next page). The map is then divided in- 
to two (1—1)-variable submaps, the divided variable 
becoming an adjunct to one submap, and the negative of 
the divided variable becoming an adjunct to the other 
submap. In each submap, the adjunct to that submap is 
written in the lower left-hand corner of each square 
representing a valid combination. (See Fig. 3.) (The 
particular order of subdivision used in the running ex- 
ample leads to a minimal tree; the basis for arriving at 
this optimum order will be apparent later, and the ex- 
ample will be reviewed from this standpoint.) 


3 E. W. Veitch, “A chart method for simplifying truth functions, ” 
Proc. Association for Computing Machinery; May, 1952. 

4M. Karnaugh, “The Map Method for Synthesis of Combina- 
tional Logic Circuits,” AIEE Tech. Paper 53-217; 1953. 
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Valid Combinations 


June 


3—First subdivision. Fig. 4—Second and third 


subdivisions. 


ABCD 
ABCD’ 
AB C'D 
Nl 3 Kk Oe) Be 
Mo st Gal Die 
A‘B’'C. D 
iN AOU DM 
ABS Err 
Fig. 2—Map related to partially- Fig. 
developed tree. 
Fig. 5—Fourth and fifth 
subdivisions. 
Each (n—1)-variable submap containing a_ valid 


combination is subdivided into two (7 —2)-variable sub- 
maps, the newly-formed adjunct in each case being writ- 
ten to the right of the previously-written adjunct. (See 
Fig. 4.) A submap containing only invalid combinations 
is not further subdivided. An invalid-combination sub- 
map containing 2” squares represents the elimination of 
2*—1 transfers. 

The subdivision process is continued until each valid 
square becomes a submap; all other submaps will con- 
tain only invalid squares. (See Figs. 5-7.) The desig- 
nation of each valid combination will be completely 
written at the bottom of its representative square. The 
transfer tree can then be drawn with reference to these 
written designations, each input-to-output path from 
left to right corresponding to the equivalent order of the 
related written combination. (See Fig. 8.) 

The order of subdivision determines the order of the 
relay contacts in the tree; the writing of the adjuncts 
in the valid squares is the means of recording the order 
of subdivision so that the transfer tree can be drawn by 
reference to the map; the subdivision process insures 
that the submaps formed will be compatible with 
transfer tree configuration. 

Any order of subdivision will lead to a legitimate tree. 
However, to obtain a desired tree having the minimum 
number of transfers, the map must be subdivided such 


Fig. 6—Sixth and seventh 
subdivisions. 


Fig. 7—Eighth subdivision. 


Fig. 8—Partially-developed transfer tree obtained from Fig. 7. 


that the number of invalid-combination submaps is a 
minimum. The procedure is to combine the invalid com- 
binations into the minimum number of groups that can 
be obtained by the subdivision process® (there may be 


>In the following map, it is possible to combine the invalid 
combinations into three groups (AB, A’CD’, B’C’D); however, 
these three groups cannot be collectively obtained by the subdivision 
process (they are not compatible with transfer tree configuration). 
The minimum number of groups of invalid combinations that can 
be obtained by the subdivision process is four. 
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more than one way of doing this), and then subdivide 
the map to form the desired groups (there may be more 
than one way of doing this, also). 

This was the procedure followed in the running ex- 
ample, as a review of the example will illustrate. 
Analysis of the map in Fig. 2 shows that the minimum 
possible number of groups of invalid combinations is 
three (A’B, B’C’D’ and AB’D), and that these groups 
can be obtained by the subdivision process. Further- 
more, to form these three particular groups, the map 
must be subdivided as illustrated in Figs. 3-5. The 
tree obtained (Fig. 8) is, therefore, a minimal one. 

The foregoing method may be modified; for instance, 
the subdivision process may be continued only until all 
invalid-combination submaps have been formed, as il- 
lustrated by Fig. 5; from the map at this stage, the tree 
can be partially constructed, and the rest completed by 
inspection. Fig. 9 shows the portion of the tree obtained 
from the map in Fig. 5. Branch removal has been com- 
pleted and the rest of the tree can be drawn by inspec- 
tion. 

Following is a summary of the relationships involved: 
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Fig. 9—Portion of tree obtained from Fig. 5. 


nm =number of relays in tree. 
2”=number of output combinations in fully- 
developed tree. 
(2”—1) =number of transfers in fully-developed tree. 
m =number of invalid output combinations. 
(2”—m) =number of valid output combinations. 
b=number of invalid-combination submaps. 
(m—p) =number of transfers eliminated from fully- 
developed tree. 
(2™—1) —(m—p) =(2"—m) +p—1=number of trans- 
fers in partially-developed tree. 


A New Diode Function Generator* 
T. MIURA, H. AMEMIYAt, anp T. NUMAKURAt 


Summary—With the diode function generators that are currently 
in use, generation of functions is made by combining straight lines. 
The main drawback of these function generators is that the slope of 
the line segment cannot be changed independently. With the new 
function generator described in this paper, functions are generated 
by connecting independent line segments. Accordingly, the slope of 
each segment is given independently and also quantitatively. It is 
possible to approximate any desired function without recourse to 
an oscilloscope for inspection. These advantages are realized by 
using ganged potentiometers differentially. The operating principle 
and a practical generator experimentally built are described. 


INTRODUCTION 
7 | YHE DIODE function generators are very useful 


as components of electronic analog computers. 

However, those generators currently available are 
not convenient to use because of their inherent disad- 
vantage that the slopes of line segments constituting a 
function cannot be changed independently. This is be- 
cause a function is generated by superposing linearly- 
increasing functions starting at certain input voltages, 
such as shown in Fig. 2(a). With the new function gen- 
erator, a function is produced by superposing constitu- 


* Manuscript received by the PGEC, October 2, 1956; revised 
manuscript received, January 14, 1957. 

+ Hitachi Central Res. Lab., Tokyo, Japan. 

t Showa Denshi, Ltd., Yokohama, Japan. 


ent functions, which increase linearly for certain input 
voltage ranges and which level off beyond, [Fig. 1(a) ]. 
Therefore, the slope of each line segment of the resultant 
function is produced in the same manner but independ- 
ently of others. Any desired function can be generated 
very easily and yet quantitatively. 


OPERATING PRINCIPLE OF THE NEW DIODE 
FUNCTION GENERATOR 


Figs. 1 and 2 show how a function is generated by the 
new and by the conventional diode function generators. 
The function to be generated is shown at the top in both 
figures and is expressed by 


[cole <ei<Ey 41 Fi (Ee y EF) tan A, 
+ (E3; — Ey) tan 62 


ar (ar = Jib San) tan Oral 
= (e; vy Ee) tan bn 


n—1 


= (Ej41 — E;) tan6é;+ (e; — E,) tan @, 
7=1 

4 De C07; (1) 
j=1 
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(a) 


OUTPUT 


En>En-1 


Fig. 1—Operating principle and basic circuit of the new diode 
function generator. 


where 

€o; = (Ey. — E;) tan 6,;, 
and 

€on = (€; — En) tan Op. 

Eq. (1) may be rearranged as 
(e; — F) tan 6; 
+ (e; — E2)(tan #2 — tan 4;) 
+ (e; — E3)(tan 63 — tan 62) 


[eo Jen< eck, +1 = 


+ (e; — E,)(tan 6, — tan 6n_1) 


es (ee 4) taniOs —tan'6,2,) 


j=1 


= Sy) (e; — E;) tan 0; 


j=1 


= do e0/ (2) 
j=1 
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@,=tan' (tan Gz - tan@1) 
| | 


Eq 


OUTPUT 
) Te) 


(b) 


Fig. 2—Operating principle and basic circuit of the conventional 
diode function generator. 


where 


€o; = (e; — E;) tan@;, tan O; = tan 0; — tan 6;-1, 


and 
tan A = () 


With the conventional function generator, the desired 
function is obtained by realizing the relation expressed 
by (2). The terms of (2) are shown separately in Fig. 
2(a). By adding these voltages, the desired function is 
obtained. 

Fig. 2(b) is a basic circuit to achieve this. E’s are in- 
dividual bias voltages, V’s are diodes, and k’s are po- 
tentiometers. The diodes conduct only if the input volt- 
age e; is higher than the respective bias voltages, thus 
giving straight lines starting at the bias voltages. The 
potentiometers are used to give slopes to these lines. 
Now it is clear that the slopes of line segments of the 
resulted function are closely interrelated. For instance, 
if tan @, is changed to vary the slope tan 62, all the 
slopes of the following segments, such as tan 63, tan 64, 
etc., change simultaneously. This is a disadvantage of 


1 Although not shown in the figure to avoid confusion, a sign 
changer is necessary between the potentiometer and the adder when 
the slope tan 0, is negative. In the example given here, a sign changer 
must be used after Re, because tan @» is negative. 
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the function generator of this type, because this makes 
it troublesome to adjust the settings of potentiometers 
to obtain precisely the desired function. 

With the new diode function generator, the function is 
generated by realizing (1) literally term by term. As 
shown in Fig. 1(a), each constituting voltage eo; has a 
slope prescribed by the function for E;<e;< Ej. and 
then levels off. The function is obtained by adding these 
voltages. Fig. 1(b) shows how this is done. Here, E’s 
are bias voltages and V’s are diodes as in Fig. 2(b), and 
R’s are ganged potentiometers. The bias voltage £;, the 
diode Yi, and the potentiometer Ry give a line with a 
slope of tan 0, starting at e;=&,. Similarly, Fo, V2, and 
Ry give a line with the same slope, because Ry and Ri» 
are ganged, but starting at e;=F». The output from 
Ry is fed directly to the adder, whereas the output from 
Ry is fed first to the sign changer and then to the adder. 
Therefore, the output of the adder contributed by these 
two circuits has the slope of tan 6; for FE; <e, < EF» and is 
constant for E,Se; as shown by eq of Fig. 1(a). Voltages 
€o2, €o3, etc. are obtained in exactly the same way. Since 
the output voltage eo is the sum of these voltages, the 
slopes of the line segments are given independently by 
individual potentiometers, and there is no relation either 
explicit or implicit, between these slopes. 

The practical diode function generator built on the 
principle described in the preceding paragraphs has 
construction somewhat different from Fig. 1(b). The 
basic circuit is shown in Fig. 3, where potentiometers are 
not grounded so that the slopes of the line segments can 
have either positive or negative values according to their 
settings. In Fig. 3, the output voltage is given by 


Ri Ro 
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INPUT E Vv, 


OUTPUT 
° Co 


Fig. 3—Basic circuit of the new diode function generator 
experimentally built. 


The equation shows that the slopes, which are given by 


1 1 \= (; 1 \= 
= ) = GTO ae 
€ L1—a/R b 1—b/R 


may take either positive or negative values and are 
given independently by respective ganged potentiom- 
eters. 


PRACTICAL DIODE FUNCTION GENERATOR 


Fig. 4 (next page) is the circuit diagram of the new 
diode function generator that gives eight independent 
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Fig. 4—Circuit diagram of the new diode function generator. 


S: SLOPE 
B: BIAS 


Fig. 5—Slope and bias potentiometers are used to control 
the segments shown. 


line segments, namely, four on the plus side and four on 
the minus side, as shown in Fig. 5. 

The constant current circuit consisting of V; and of 
its associated potentiometers, By, By, B3, and Bs, and the 
cathode followers, V2 and V3, replace the individual 
sources of bias voltage, £,, £2, etc., in Fig. 3.2 A perfect 
constant current characteristic insures the same gain 
from the input to the grids of all the cathode followers, 
V, and V3, which is very important to obtain the afore- 
mentioned independency of slopes. Capacitors, Ci, Co, 
C;, and C,;, are used to maintain the same gain up to 


® The circuit was designed after Function Fitter, manufactured by 
G, A, Philbrick Researches, Inc,, Boston, Mass. 
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Fig. 6—Calibration curve of slope vs potentiometer dial setting. 


high frequencies where stray capacities tend to decrease 
it. Thirty kilo-ohm, 5-turn helical potentiometers of 
good linearity are used for slope potentiometers. The 
internal resistance of the diode circuit, that is, the re- 
sistances of the cathode follower and the diode rectifier 
in series, must be low to obtain the independency of 
slopes.* Practically, however, this requirement is some- 
what alleviated by equalizing the diode loading by at- 
taching Kz and Jy, to the first diode. 

Functions for negative values of e; are generated in 
exactly the same manner as for positive values, with the 
only difference being that a sign changer, Aj, is used at 
the input. Potentiometers, Z,; and Z», are zero (origin) 
adjustments of the bias potentiometers, and Cy and Cio 
are to compensate for the phase lag of the operational 
amplifier A>. The output dc level of A» is adjusted with 
the input of A; grounded as otherwise the output dc 
voltage of A; feeds back to the input of A» making the 
level adjustment impossible. 

Fig. 6 shows the relation between slope and poten- 
tiometer setting. The curve was experimentally ob- 
tained and applies to all the slope potentiometers with 
very little error. It is possible to set desired slopes by 
referring only to this curve, without recourse to an os- 
cilloscope for inspection. Bias voltage is 8.6 volts when 
the bias potentiometer is placed at the full-counter- 
clockwise position. 


Fig. 7 shows a function generated by the new func- 


3 Appendix I, 
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(a) 


(b) 


Fig. 7—Independency of slope and independency of bias are 
shown in 7(a) and 7(b), respectively. 


tion generator. Fig. 7(a) illustrates visually the in- 
dependency of slopes. The picture was taken by double 
exposure changing the settings of potentiometers, S) and 
S»’. It should be noted that the other slopes are not 
changed at all. Fig. 7(b) shows the independency of 
bias adjustment. Bias potentiometer B; was turned with 
a mere shift of a part of the generated function. In Fig. 
8, the absolute value characteristic and square charac- 
teristic are approximated. Fig. 8(a) and 8(c) are by 400 
cps input; 8(b) and 8(d) by 1 kc input.. It is supposed 
that the phase shift at 1 kc is about one degree or two. 


CONCLUSION 


A new diode function generator described in this 
paper is believed to be a very powerful addition to 
analog computers because of its extreme versatility and 
also of its easy but quantitative generation of any de- 
sired function. The principle on which this new function 
generator is based may also be applied to obtain a 
limiter with ideal characteristic.* 
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If the internal resistances of the diode circuits are 
taken into consideration, (3) becomes 
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(a) (b) 


(c) (d) 


Fig. 8—Absolute value characteristic and square characteristic 
are approximated. Figs. 8(a) and 8(c) are by 400 cps input, and 
8(b) and 8(d) by 1 ke input. 


Rg = a1 — aR, 
a(i — a)b(1 — 3d) 


Bo ’ 


a(1 — a) + d(1 — 3) 


and 
m(1 — m)n(1 — n) 


m1 — m) + n(1 — n) 


Ry =a 


Independency of slopes of line segments is attained if 
the condition 


Ry Rp, Ry 
Ru SF Ra Ris aF Rez Rip SF Ryn 
is satisfied. Ideally, Ry =Re= -- - =Rin=0 is the only 


solution. If the resistance of potentiometer is made high, it 
is possible to realize Ra>>Ra, Re>Ria, , Rv >Rin to 
satisfy the above condition. This is not recommended, 
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aP (5 a i coulis \ ie ) ; 
Ret Re Riz+ Re Rig + Re Rio — Re b yh Ik 
E,) Ry \( 1 1 \= 
C= n Y , 
if Rin —- Ry nN 1 SS R 


where R,’s are the internal resistances and 
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however, for a function generator of high-speed analog 
computers for which good high-frequency response is 
essential. 
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APPENDIX II 


Ordinary limiter circuits give a final slope of 


| =| R; 
ieee RL: 


as shown in Fig. 9. To idealize the characteristic and 
make the final slope zero, R must be made very large in 
comparison with R;, which limits the useful frequency 


=a 


range. 
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Fig. 9—Ordinary limiter circuit and its characteristic. 


If a circuit is arranged so that the output from a dead 
zone circuit is subtracted from the input as shown in 
Fig. 10, a limiter with ideal characteristic, or even with 
a negative final slope, is obtained. The frequency re- 
sponse of this ideal limiter can be made sufficiently wide 
to cover the necessary frequency range of repetitive 
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Fig. 10—Block diagram of the new limiter with ideal 
characteristic. 
analog computers. It is recommended that a switch be 
provided so that the equipment may be used both as a 
limiter and as a dead zone simulator. 
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An Electronic Analog Multiplier* 


DAVID’ Cy KALBRELL} 


Summary—This multiplier uses the variable pulse area principle, 
but employs phase sensitive circuitry to operate naturally in all four 
quadrants without bias voltages. The output is zero if either input 
is zero. The X and Y channels are separately linearized with inde- 
pendent feedback loops. The circuitry is simple and lends itself to 
either transistors or vacuum tubes.! 


GOOD electronic multiplier should have s2veral 
A characteristics: 1) It should operate correctly in 

all] four quadrants; that is, X and Y may be 
either positive or negative. 2) Bias voltages should not 
be added to the variables, to permit four-quadrant 
operation, since these bias voltages may drift and cause 
errors. 3) The rise time for transient changes should be 
no greater than 10 milliseconds. 4) The circuit should be 
inherently free of drift and be fairly linear without feed- 
back. 5) Negative feedback should be incorporated in 

* Manuscript received by the PGEC, October 6, 1956. 


+ Kalbfell Electronix, San Diego, Calif. 
1 Patent pending on this device. 


the X and Y channels of the multiplier separately, so 
that each is independently linearized and compensated 
for drift. 6) The whole circuit should be simple and re- 
quire a small number of vacuum tubes. 


SYMMETRICAL 
LIMITER 


VOLTAGE 
SENSITIVE 
OSCILLATOR 


MONOSTABLE 
MULTI- 
VIBRATOR 


BALANCED 
AMPLITUDE 
MODULATOR 


LOW 
PASS 
FILTER 


XY 


INDICATING 
SYNCHRONOUS 
RECTIFIER 


Fig. 1. 


Fig. 1 illustrates the basic operating principles of this 
new electronic multiplier. One variable, X, controls the 
width of pulses; while the other variable, Y, controls 
the height of the pulses. Phase sensitive circuitry is 
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employed to keep track of the signs of X and Y, so that 
bias voltages are unnecessary. The X voltage controls 
the frequency of a voltage-sensitive oscillator, such as a 
multivibrator. This might typically operate at a fre- 
quency of 100 kc. The output of this voltage-sensitive 
oscillator fires a monostable multivibrator. The duration 
of the metastable state of this multivibrator is deter- 
mined by an inert delay line. If the voltage-sensitive 
oscillator puts out a frequency of 100 kc when X =0, 
then the delay line should be 5 microseconds long. 

This unsymmetrical waveform is now passed through 
a blocking capacitor and is clamped symmetrically with 
respect to ground. This symmetrically clamped wave 
wil] now contain a dc component which is directly pro- 
portional to the original X modulating voltage except 
for such nonlinearities as might have been introduced by 
the system up to this point. The de and low-frequency 
components of this unsymmetrical wave are now fed 
back to the input as negative feedback. For ordinary 
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factor might be employed and the frequency of the 
voltage-sensitive oscillator could be deviated through 
a wider range, thus making the latter portion of the 
equipment function more accurately. 

The first section of this system functions as a unit 
including its own stabilization, and merely feeds a 
clamped unsymmetrical waveform to the amplitude 
modulator. This modulator converts the clamped am- 
plitude into a variable amplitude which is proportional 
to Y. This might be accomplished in a ring modulator 
or by several alternative means. The output of this 
modulator is used in two ways. First, a feedback signal 
must be derived from it to buck against the Y signal 
to give negative feedback around the amplitude modu- 
lator section of the system. In order to give the correct 
sign and to be independent of X, the amplitude of the 
modulator output is sensed and this is passed through a 
switching system which insures that the correct sign of 
feedback voltage is always bucked against Y, regardless 


LIl Ll au 


SET FULL SCALE 
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FILTER 


4X Y OUTPUT 


READ BAL 
© 


READ Yo 


Fig, 2. 


computer applications, X would have a full-scale value 
of 100 volts, although 10 volts would probably be suffi- 
cient to produce the desired degree of modulation of the 
oscillator. If the feedback voltage is 90 volts when the 
input voltage is 100 volts, this will act like 20 db of 
negative feedback and will compensate drifts and non- 
linearities. For a superior system, the X voltage might 
have been fed to the voltage-sensitive oscillator through 
a preamplifier. If the feedback voltage is mixed with X 
at the input of the preamplifier, then a larger feedback 


of the original sign of Y. One method of accomplishing 
this synchronous rectification will be described later. 
The second use to which the output of the amplitude 
modulator is put is to give the product of X and Y after 
passing through a low-pass filter. The output of this 
low-pass filter may be either positive or negative, de- 
pending upon the original signs of X and Y. If either 
X or Y is 0, the output will be 0. 

Fig. 2 is a more detailed circuit diagram showing one 
method for making an electronic multiplier of the type 
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which has been described. The X input is applied to the 
bottoms of the grid leaks of a multivibrator operating 
at a frequency of 100 ke. In practice, it might be desir- 
able to make a somewhat more sophisticated form of 
multivibrator in order to give better stability and 
sharper cornered waves. This is an engineering detail 
which does not concern the present design, however. 
The multivibrator output is buffered through a cathode 
follower and triggers the monostable multivibrator. TS 
is normally cut off, since its grid leak is returned to a 
negative bias voltage. A positive pulse from T3 causes 
TS5 to fire, which cuts off T6. When T5 fires, its plate 
voltage drops and a negative pulse is fed through T4 
to a delay line. The delay line, D, should be 5 micro- 
seconds long for the example being discussed. When the 
negative pulse comes out of the delay line, it cuts off 
T5, thus resetting the monostable multivibrator to its 
initial state. When the negative pulse comes out of the 
delay line, its energy is partly drained to ground through 
the diode and resistor R4. This need not cause any 
trouble, however, since R4 may be picked to be the 
optimum terminating resistor for the delay line. 

The waveform sketches which are indicated through- 
out this figure are drawn on the assumption that X is 
positive and Y is positive. 

The voltage at the cathode of T4 is fed through a 
series resistor and capacitor to a symmetrical clamp. 
This clamped voltage is fed through a low-pass filter, 
cutting off in the vicinity of 1 kc, back to the input 
terminal. Since the cathodes of T1 and T2 are tied to a 
negative potential, their grids would be negative and 
current would flow through their grid leaks from an in- 
put point at ground potential. R2 is inserted to provide 
current from a positive supply to just equalize the nor- 
mal grid current of T1 and T2, so that the signal input 
point is normally at ground potential. If X is positive, 
it causes current to flow into these grid leaks. The feed- 
back signal, however, is then negative and causes a 
negative current to flow into the grid leaks, largely 
cancelling the X input current. This system may be ad- 
justed for drift by setting the X input point at ground 
and adjusting the variable resistor labelled “Set Xo” to 
make the potential at jack “Read Xo” come to ground 
potential. This part of the system is then adjustable 
and tends to be self-compensating. 

The voltage from the cathode of T4 is now fed to the 
primary of a modulating transformer. For best stability, 
this voltage may well be taken from the point at which 
it is symmetrically clamped with respect to ground, in 
order to avoid amplitude variations due to tube aging. 
The Y input signal is applied to one side of this balanced 
ring modulator, and the Y feedback voltage is applied 
to the other side. The ring modulator is balanced by 
potentiometers at the center of each winding. This 
makes it possible to compensate for differences in the 
diodes. The polarity of winding of all transformers is 
indicated by dots at one end of each winding, and the 
waveforms are also indicated on the assumption that 
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both X and Y are positive. The secondary voltage at 
the output of this ring modulator will have the same 
general type of wave shape as the input, although it 
would be inverted if Y were negative. In any case, its 
de component has been lost in passing through the 
transformer, and so it must again be clamped sym- 
metrically with respect to ground. Since the amplitude 
now contains part of our intelligence, however, we can 
no longer use a simple clamp for this purpose. Instead, 
peak indicating rectifiers are used to develop positive 
and negative voltages corresponding to the positive and 
negative amplitudes of the wave after passing through 
an operational amplifier. By feeding these two peak 
voltages into the input of the amplifier, the output is 
made symmetrical with respect to ground. This takes 
care of the clamping requirement. 

The next portion of the circuit provides negative feed- 
back voltage to correct for any nonlinearity or drift in 
the ring modulator and the operational amplifier. The 
pulses on the output of the operational amplifier are 
amplified further, so that they may be used even though 
their amplitude be less than one tenth (0.1) volt. T7 
performs this amplification, and its plate current flows 
through the primary of another modulator transformer. 
If desired, these transformers could be shunt fed in- 
stead of carrying de current. The other transformer on 
this modulator is fed from the cathode of T4 and pro- 
vides a de return path for the current of this tube. This 
modulator is balanced so that when the output of the 
operational amplifier is shorted to ground, 0 voltage is 
produced at the jack labelled “Read Bal.” In order to 
provide appropriate feedback voltage which depends 
upon Y and not upon X, we must extract information 
which depends only upon the amplitude of the output 
of the operational amplifier, while being independent of 
its waveform. Since this output has already been 
clamped symmetrically with respect to ground, either 
peak indicating rectifier might be used for this purpose, 
but it is also essential that the correct sign be selected. 
This balanced modulator and the diode gate system 
which is fed by the modulator act as a synchronous 
switch to permit signal of only the correct sign to be ap- 
plied to the Y feedback point of the Y amplitude modu- 
lator. 

The operation of this synchronous switch will now be 
described. If the phases of all voltages are as indicated 
for X and Y, both being positive, then the balanced 
modulator of this synchronous switch system will put 
out a positive voltage. In this case, the left-hand pair 
of diodes of the synchronous switch will be clamped to 
ground, and the negative rectified voltage of the opera- 
tional amplifier output will be shorted to ground and 
will make no contribution to the feedback signal. A 
positive signal from the modulator portion of this syn- 
chronous switch system would block off the right-hand 
pair of diodes so that they conduct no current, and 
hence the positive rectified voltage from the output of 
the operational amplifier will be applied to the Y feed- 
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back point. The proper sign of voltage will always be 
available. This system has the advantage that the mag- 
nitude of the feedback voltage is not dependent upon 
the nonlinearities of the switching system, but this am- 
plitude will simply be proportional to the voltage on the 
peak indicating rectifiers. The rest of this feedback sys- 
tem merely acts as an on-off switch of the shunting type. 

The final output of the multiplier consists of the 
average value of the output of the operational amplifier 
after passing through a low-pass filter. It may be ad- 
justed to give 100 volts out when X and Y are each 100 
volts, by adjusting the feedback resistor of the opera- 
tional amplifier. 
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From this analysis, it is now apparent that this elec- 
tronic multiplier system meets all of the requirements 
which were set forth in the first paragraph. If this ring 
modulator should get out of balance, it could only lead 
to an error for very small values of Y. Even with nega- 
tive feedback, there will be some nonlinearity of the Y 
modulator at small amplitudes, and this is corrected by 
inserting some crystal diodes in series with the Y input, 
to compensate this error. The X and Y systems are 
separately stabilized with negative feedback, and pro- 
vision is made for setting their 0 values. The system 
automatically operates in all four quadrants without the 
use of any bias voltages. 


An Algorithm for Determining Minimal Representa- 


tions of a Logic Function’ 
BERNARD HARRIST 


Summary—For each logic function, or Boolean algebraic expres- 
sion, there corresponds an appropriate computer circuit. However, the 
minimization of the appearances of the Boolean variables does not 
necessarily lead to the most economical circuit. A general approach 
to the problem therefore requires the development of techniques for 
the simple and rapid generation of a variety of near-minimal forms. 

This paper describes such a method for constructing the minimal 
representations of a logic function given as a truthtable or in one of 
its canonical forms. The minimal representations achieved are either 
sums of products, or products of sums, such that no term contains 
superfluous variables and such that no term is superfluous. The utility 
of the method lies in the conciseness of notation, which permits the 
handling of a large number of variables and simplifies the process for 
machine computation. 


INTRODUCTION 


OR EACH logical truth function, or Boolean al- 
[ieee expression, there corresponds an appropri- 

ate switching circuit. However, the minimization 
of the appearances of the Boolean variables does not 
necessarily lead to the most economical circuit. A gen- 
eral approach to the problem therefore requires the de- 
velopment of techniques for the simple and rapid gen- 
eration of a variety of near-minimal forms. 

The direct application of Boolean algebra [1] is not 
completely satisfactory, since it involves the nonsys- 
tematic trial of a very large number of alternate pro- 
cedures. A minimization chart, such as developed at 
the Harvard Computation Laboratory [2], represents 


* Manuscript received by the PGEC, October 29, 1956; revised 
manuscript received, January 12, 1957. The research reported in this 
paper was performed at New York University and sponsored by 
Air Force Cambridge Res. Ctr., Air Res. and Dev. Command, under 
Contract AF 19(604)-1049. 

+ College of Engineering, New York University, New York, N. Y. 


an improvement, but the large number of entries renders 
the system cumbersome. The map method of Karnaugh 
[3] represents a further improvement because it permits 
the over-all system to be visualized. Unfortunately, 
however, these methods and their variations become 
difficult to apply when the number of independent 
variables greatly exceeds five. 

The technique of assigning a decimal number to each 
minimal polynomial of the truthtable considerably im- 
proves the facility with which a large number of vari- 
ables may be handled. The method introduced by 
Mueller [4] requires performing arithmetic operations 
with each decimal representation and checking whether 
the numbers so derived are present in the given list of 
decimal vertex numbers. By these comparisons, the 
decimal numbers are appropriately grouped into forms 
which permit expression by the fewest possible terms. 
The main disadvantage of this method lies in the neces- 
sity for the separate bookkeeping of indexes to deter- 
mine the form of the solution. 

The method described here is an extension of Mueller’s 
idea, the terms of the canonical form being interpreted 
as though they were ternary numbers rather than 
binary. After these ternary numbers have been written 
as their decimal equivalents, the resulting numbers are 
manipulated without the necessity of having to refer to 
the truthtable again. The use of ternary notation per- 
mits not only the original entries of the truthtable to be 
assigned corresponding numbers, but the various groups 
as well. The solution then consists of a set of ternary 
numbers which directly corresponds to the appropriate 
groups. 
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THEORETICAL DISCUSSION 


There exists a one-to-one correspondence between 
any possible terms of the canonical form of an n-inde- 
pendent-variable logic function and the vertices of an 
n-dimensional cube. The terms of the canonical form 
therefore correspond to a unique subset of the vertices 
of the n-dimensional cube. Let any face of this cube, all 
of the vertices of which are present in the truth function, 
represent a cell. The logic function may then in turn be 
represented by any set of cells which contains all of the 
required vertices at least once. 

An irredundant covering of the truth function is a 
set of cells such that each vertex is included in at least 
one cell, and such that no cell is contained in a larger 
cell or in the covering of the other cells. This permits 
the truth function to be expressed as a sum of products 
such that no product is superfluous or contains a super- 
fluous variable. By covering the complement of the logic 
function in a similar fashion the truth function may be 
expressed as a product of sums such that no sum is 
superfluous or contains a superfluous variable. In gen- 
eral, several such coverings exist, and the emphasis is 
on determining those combinations having the smallest 
possible number of cells. 

Description of each vertex may be given by the ap- 
propriate binary number. Description of a cell, however, 
requires another symbol to indicate which variables are 
not required to describe the cell. This requirement may 
be met by using a ternary radix. 

A simple coding scheme between the conventional 
binary notation and this ternary notation is to let the 
ternary zero and one correspond to the binary zero and 
one. The ternary two then corresponds to a condition 
in which both the binary zero and the binary one are 
included. Vertices are thus represented by all possible 
combinations of the ternary zeros and ones, while an 
r-dimensional cell is characterized by a ternary number 
containing the digit two in exactly 7 positions. Varia- 
tions in the coding scheme are possible, but this choice 
permits the vertices to be designated by the smaller 
numbers and the higher order cells by the generally 
higher numbers. 

As an example consider four vertices which may be 
combined into a second order cell. 


101000 


101001 
101202. 
101100 


101101 


Considering these numbers to be in ternary notation 
and writing the equivalent decimal numbers gives 


270 + 271 + 279 + 280 = 290. 


The rule for this “addition” is that the “sum” of the 
numbers equals twice the largest minus the smallest. 
It is important to note that this rule is only a procedure 
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for determining the proper number representing the 
(r+1) cell after it has been determined that the proper 
vertices are being combined. The remaining problem 
lies in determining if the proper vertices are present. 

Consider two vertices, numerically symbolized by 
conventional one and zeros, to be separated by a dis- 
tance defined as the number of binary places by which 
they differ. The binary representation of these vertices 
may further be considered to be actual numbers written 
with some radix P. Then the necessary condition that 
these vertices be exactly a distance K apart in u-dimen- 
sional space is that the numerical difference of the cor- 
responding vertex numbers be a number expressible in 
the form 

ES 
> asP* (1) 
t=1 
where k;#k; is the degree of the positions where the 
numbers differ and a; is either plus or minus one. 

To show that this difference is a sufficient indication 
of the distance between vertices, consider any two 
vertices which are separated by a distance H#K. The 
numerical difference of the corresponding vertex num- 
bers must then be of the form 


>, BPM (2) 


j=1 


and it is necessary to show that 
K H 
D> aP* % DY BP, (3) 
1=1 411 


Assume that an equality is possible and transfer those 
terms for which a;= —1 to the right and those terms for 
which 6;= —1 to the left. If an equality is possible, then 

PORE sein ahi Dea GALI Nia Mh 2 ie tay ei cig? a, 

iaj=l i’ ,Bj’=—1 Vag’ 5,8j=1 

However, the expressions on either side of this equation 
represent unique numbers in a P-base number system 
for P23. These numbers may be designated by numbers 
having “1” in those positions where k;%h;, “2” where 
k;=h;, and zero elsewhere. Since H#K, however, the 
resulting numbers must differ in at least one place. The 


equality cannot, therefore, be satisfied. These results are 
restated as a theorem. : 


Theorem: If the binary numbers corresponding to the 
canonical form of a logic function are interpreted as 
numbers of radix P>3, a numerical difference of 
value 


K 
Di a,P*i, Y= ae 1 
ci 
is a necessary and sufficient condition that the ver- 


tices are exactly K apart in m dimensional space. 


While this theorem is applicable for any radix P>3, 
greatest conciseness of notation is achieved by choosing 
a radix 3. This follows from the fact that the numbers 
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are to be recast into numerically equal decimal num- 
bers. The largest number employed using a radix P 


with m independent variables will be 222 + - - 2, or 
Sel b 2 pr ie ie 
2 aS Pi= ac tinaae) : 
j=0 (Peael) 


If this number is recast into decimal notation the result- 
ing number of digits is 


OF 2(P" — 1) 
x = logig ———-—— 
nude (he at) 
or, for large n, 
x 
a logio Ie, 
n 


Hence, if a ternary radix is used, the number of deci- 
mal digits required is approximately 0.48 times the 
number of independent variables, while for a quaternary 
notation this ratio is 0.60. 


METHOD FOR DETERMINING [RREDUNDANT COVERINGS 


Based on the work of Samson and Mueller [4-6], this 
ternary notation can be used to determine the irre- 
dundant coverings of a truth function. To do this it is 
necessary to determine for each vertex those cells which 
contain the vertex but are themselves not contained in 
a larger cell. Such cells will be called basic. To do this, 
first write all the vertices in ternary (binary) form, and 
then write the corresponding decimal numbers. 

For a given vertex JN, only those vertices whose vertex 
numbers are of the form N+3" can form one-dimen- 
sional cells containing V. Now ¢ one-dimensional cells 
define an r-dimensional cell provided all vertices of the 
form below are present. 


N+ >> a3", 
t=1 
Consider two vertices 
Ny = N + ashi 
N» N+ a3 


which form one-dimensional cells with the vertex JN. 
For these vertices, plus the vertex N, to form a two- 
dimensional cell, the vertex below must also be present. 


Nw= N+ ashi + aol 
Nye = Ny + Ne >; N 


For three vertices 


Ni = N 4 a3", 
Ne = N+ a;3", 
Ns = N + a3", 


which form a one-dimensional cell with the vertex NV to 
form a three-dimensional cell, the three two-dimensional 
cells requiring the additional vertices 
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Function 


Nw=Ni+tNe—N 

Ni3 = Ni+ N3 — N 

Nos = Ne+ N3 — N 

must be present, as well as the vertex 

Neat eta 37 aah ao 
or 

Niex = Nia + Ng — NV 

= Ni; + No—N 
= Nea+ Mi — KN. 


This process is readily extended to higher order cells. 

Consider a vertex which is embodied in only one 
basic celi. Such a basic cell is essential and must be 
present in any irredundant covering. Furthermore, the 
vertices embodied in this cell are thereby also covered 
and need not be separately considered. This results in a 
considerable decrease in the labor involved in deter- 
mining an irredundant covering. 

Consider those vertices which are not contained in a 
single basic cell. Tabulate these vertices together with 
the various basic cells in which they are embodied. If a 
larger cell covers the same vertices as a smaller cell, the 
smaller cell is not essential and may be deleted. If, asa 
result of this deletion, a vertex is now contained in only 
one remaining basic cell, this cell is essential. 

The possible variations in the covering of the remain- 
ing vertices form the possible variations of the irre- 
dundant form. To form an irredundant covering con- 
taining the minimum number of cells, preference is 
given to those cells covering the largest number of the 
remaining vertices. 


Example 1 


The method will be applied to determine the irre- 
dundant form of a four-variable case. Consider the 
truthtable 


0000, 0100, 0110, 0111, 1000, 1100, 1101, 1111 


interpreted as numbers in ternary notation. The corre- 
sponding numbers in decimal notation are 


O90 p12, 15, 21, 50 90F waOs 


To determine the basic cells embodying the vertices 
start with the vertex 0. For this vertex, the vertices of 
the form 0+3? and 0+3% are present, so that these 
vertices form one-dimensional cells with the vertex 0. 


hh Oa 


Since there are two one-dimensional cells present, there 
is a possibility of a two-dimensional cell being present. 
This requires that the vertex (9+27—0) or 36 be pres- 
ent. Since this vertex is present, these four vertices are 
embodied in a single two-dimensional cell. The numeri- 
cal value of this cell is given by twice the largest vertex 
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minus the smallest or 72. Symbolically the process is 


36 Ce 


The circle about the cell number indicates that it is the 


written 
Cio 274 


only basic cell containing the vertex and hence is es- 
sential. The basic cells embodying 9, 27, and 36 need 
not be considered, although these vertices may be em- 
ployed to simplify the cells about other vertices. 

The remaining vertices form groups as follows 


12% fe 9. Tg 
(3/12; ae 
37 36,4840 
AOS 3 Tee Sif). F 


Since no additional essential cells appear, the various 
possible coverings of the remaining vertices must be 
considered. To do this, list the vertices and the cells 
covering them. 


Vertices One-dimensional cells 
12 15 14 
is 14 67 
ST 38 43 
40 43 67 


To find the irredundant coverings that employ the 
fewest possible cells, preference is given to those cells 
which cover the greatest number of vertices. Thus a 
minimum irredundant covering is 


72+ 14+ 43. 


Other irredundant coverings are found by covering each 
vertex at least once, with each cell covering at least one 
vertex not otherwise covered. 


15+ 67+ ie 
72+ 38 


{4-+ 38 + 67. 


If desired, these forms may be put into truthtable 
notation by converting the decimal numbers to a 
ternary notation. Thus, the minimum irredundant form 


becomes 
2200 + 0112 + 1121 
or 
Xy/X9/ + 49’ H2%1 + X3X2X0. 


Example 2 


The following example has seven variables and 43 
terms in the canonical form. It is identical to the ex- 
ample described by Mueller [4], who generated it by 
reading modulo 128, the first 45 page numbers in an 
alphabetical index of a 2000-page volume. The entire 


calculation is shown below. 


PLO no Ooms, Ti2otdy TLS el 2O0n 1216243" 24-7, 
255 0200,280,3283; 325,336; 351 63525 3545:300,,361.-363, 
729, 730, 732, 133, 739, 756, 760, 765, 810, 811, 849, 973, 
975, 1000, 1002, 1008, 1054, 1065. 
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3 / 732 /461) 


100/253. 73008 
36 / 117 765 / 
93 / 120 336 / 363 /(633) 


109% 1175 351° / 118% 352# 360% / 361 / 14) 


109 121 / 118 / 33) 


283 / 280 256 / 253 /G13 

325 / 352 1054 / 

354 / 351 363 / 360 /G75) 

739 / 730 10 / 

729 / 730 732 756 810 / defer 
730 / 729 733 739 973 811 / defer 
733 / 732% 730° 760° / 729% / 
756 / 765 729 / 

760 / 733 /@87) 

765 / 756 1008 36 / 

810 / 811 729 / 730 /(93) 
849 / 120 / 4578) 

973 / 1000% 1054 730° / 811% / 
975 / 1002 732 / 

1000 / 973 / G02 

1002 / 975 / 4029 

1008 / 765 / (251) 

1054 / 973% 811% 325¢ / 730% / 
1065 / 336 / (794) 


The calculation about those vertices embodied in 
many one-dimensional cells is deferred until the end be- 
cause of the high probability that these vertices will 
appear in an essential cell. The calculation of the nu- 
merical value of the basic cells about a vertex that has 
more than one basic cell has also been deferred. This is 
done because these vertices may later appear in an es- 
sential cell and thus require no further consideration. 
Letter superscripts are employed, where necessary, to 
note what vertices are being considered. 

Table I is formed by collecting those vertices not 
covered by essential cells. 
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TABLE I the average value of the numbers representing the ver- 
wos arene pean tices. While this coding has some elegance, it is never- 
theless more work to average several numbers than to 
Be foe a perform the operation of “twice the largest minus the 
325 379 1783 smallest.” 
ie i ares In some instances it may be desirable to perform oper- 
1054 1783 1378 ations directly with the cells. The following theorem is 


Smaller cells which are dominated by larger cells (in 
the sense that there is a larger cell which covers at least 
all those vertices covered by the smaller) are deleted. 
In this example the larger cell 1378 does not dominate 
1783. 

Various irredundant coverings are given by all of the 
essential cells plus those combinations of cells from this 
table which cover each vertex once. The minimum ir- 
redundant covering (in the sense of the least number of 
cells) is found by giving preference to those cells in this 
table which cover the largest number of vertices in this 
table. Thus 1468 and 1783 must be included in any mini- 
mum irredundant covering. 

The irredundant coverings are given by the cells 
133 + 243 + 265 + 213+ 375 + 614+ 633 + 787 4+ 893 + 
1027+10294+ 1251414614 1578+1794 and 


(1468 + 1783 — minimum 
1468 + 1378+ 379 
496+ 748+ 1783 

496+ 748+ 1378 + 379. 


A: { eee on 
1494 783 


In conventional truthtable notation the minimum ir- 
redundant forms are 


133 

243 

265 

els 

375 

614 

633 

787 

893 

1027, 
1029 
1255 = 
1461 al) 
1468 = AU) 
1578 ==) 
1 

1 


Sorooue. 
S || 
re Oo | 


I ae ee > > SP 
— 
| 

oO rF FP OF 


RPBPrRerRerRrRocooCcoCcoCc Oo 
| 
re oO | 


1783 ss 
1794 = 


198 or 1494 00 — 1100 or — 001100 
774 or 783 1001 — 00 or 100 — 00. 


eEPeoc ol\lco © | 
POoOrRrPROrROCOCOrF, | 
ore Or = Ol One es aS || 
Ororoco9$ce | 


ALTERNATE PROCEDURES 


It is possible to develop an alternate procedure where- 
in the ternary one and two are interchanged. The 
number designating a particular cell is then simply 


then of interest. 


Theorem: If the binary numbers corresponding to 
the canonical form of a logic function are interpreted 
as numbers of radix P>4 and the corresponding cells 
are designated by the proper locations of the digit 
two for those variables not required to describe the 
cell, then the necessary and sufficient condition that 
two r-dimensional cells form a larger (7+1) cell is 
that the difference of their corresponding r-cell num- 
bers be of the form P”. The proper (r-+1)-cell num- 
ber is then given by twice the larger r-cell number 
minus the smaller. 


When a radix 3 is used, it is possible to have two 7 
cells with a difference of 3” which have an (r—1) cell 
in common. Such +7 cells effectively form a “corner” 
of the m-dimensional cell and are thus not combinable 
into an (r+1) cell unless the appropriate additional 
(r—1) cell is present. 

This theorem permits the entire cell structure of a 
truth function to be calculated and displayed in com- 
pact form. While it is possible to utilize a radix 3, the 
additional precautions required appear to make a radix 
4 more desirable. 


Example 3 


The vertices chosen for example 1 will be again util- 
ized, but now the entire cell structure instead of just 
the minimum form will be determined. 

The corresponding decimal numbers of the vertices, 
using a radix 4, are 


Or6. 20021. O40 S081 oo. 


The results of the calculation are shown in Table II. 


TABLE II 
Vertices 1-cells 2-cells 

(0) oe 128 160 
16 BY 24 144 160 
20 22 24 
21 22 149 
64 96 128 160 
80 96 82 144 160 
81 89 82 
85 89 149 


Calculation of the 2-cells requires searching the list 
of 1-cells. This constitutes one of the main disadvan- 
tages of the method because the 1-cells are not in numeri- 
cal order. However, the main advantage is that the cell 
structure of the truth function is directly displayed. 
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Computing Techniques for the Sampling 


Parametric Computer” 
C. J. HIRSCH# anv F. C. HALLDENT 


Summary—This paper describes novel calculating techniques 
suitable for an electronic analog computer using exponential dis- 
charges to simulate the logarithmic scales of the slide rule. Among 
others, the device can perform the following operations: z=xy; 
z=x/y; z=xy"; log. x; a*; and evaluate such series as y= Ax*+ Bx’ 
+Cx+ +--+ where the exponents can be fractional and z, x, and y 
can be time variables. The problems can be solved explicitly or 
implicitly; thus in the above series y (or x) can be determined with 
equal ease if x (or y) are given. The Appendix describes an actual 
computer which gave accuracies of 1 to 2 per cent of full range from 
0 to 95 per cent of full range and 4 per cent from 95 to 100 per cent of 
full range without recalibration for different groups of problems. The 
accuracy can be increased by calibrating the device for a specific group 
of problems. 


HISTORY 


HE COMPUTING techniques described in this 
| eee were originally developed for a RHO- 

THETA navigational electronic analog com- 
puter using distance measuring equipment (dme) to 
obtain distance (RHO) and a VOR receiver to obtain 
bearing (THETA). The distance information is sup- 
plied by the dme as a de voltage whose value is propor- 
tional to distance. A computer capable of handling this 
data directly was then developed. Since its original 
concept the computer has been greatly expanded and 
diversified. It has been used for navigation, equation 
solving, the generation of unusual mathematical func- 
tions, and fire control problems. 


INTRODUCTION 


This paper describes calculating techniques which are 
suitable for a novel type of electronic analog computer 
* Manuscript received by the PGEC, Feburary 20, 1957; re- 


vised manuscript received, February 28, 1957. 
+ Hazeltine Res. Corp., Little Neck, N. Y. 


known as the sampling parametric computer. These 
techniques permit computations according to all the 
operations of algebra and, in addition, provide simple 
means for certain calculations which heretofore have 
not been conveniently performed. Among the latter 
are the evaluation of power and transcendental series 
and the solution of higher-degree equations whose 
exponents need not be integers. The computer operating 
with these techniques provides the convenient genera- 
tion of many of those functions which can be expanded 
into series. 

For example, the computer uses circuits which can 
calculate y or x with equal ease in the equation, 


is I Nn eV adhd nal Ore es Ou 


(where x and y are variable functions of time and a, 8, ¢, 
etc. need not be integers). 

Since this computer can use fractional exponents, it 
can compress long series into only a few terms which 
closely approximate the function and require only 
simple circuitry to implement, such as 


Sin ey Sa — 479/644, 


It can calculate more complicated functions with more 
elaborate circuits. 

As used here, an electronic analog computer is one that 
performs mathematical operations on physical quan- 
tities whose magnitudes are expressed in electrical 
analogs, such as volts per mile, volts per kilogram, etc. 
In addition, the analog computer handles magnitudes 
expressed in terms of the complete class of real numbers, 
integers and fractions. While less accurate than the 
digital machines, its simplicity of programming and 
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rapidity of operation make it desirable for many prob- 
lems of control and computation for which large digital 
machines are cumbersome. The sampling parametric 
computer utilizes inputs expressed as voltage analogs. 
It carries out the computation by electrical means, 
and its output is a voltage which is an analog represen- 
tation of the magnitude of the physical quantity de- 
sired. 

The techniques described here are similar to the re- 
peated observation of a slide rule and the circuits in- 
volved are electrical analogs of the scales and indicators 
of the slide rule. The magnitudes engraved on the 
different scales of the slide rule can be expressed by 
simultaneous equations in terms of the common param- 
eter of distance along the scale. The computing tech- 
niques described in this paper similarly make use of 
magnitude scales which are disposed in time, instead of 
distance, and can be expressed by simultaneous equa- 
tions in terms of the common parameter of time. 

The quantities used in most control operations vary 
with time. Computations involving such variables are 
performed by rapid repetition of the computation on 
instantaneous values of the variable which are assumed 
to be constant during each computation. As soon as 
the computation is completed, the process is repeated 
on anew sample of the variable. Hence the name sam- 
pling parametric computer. 


PRINCIPLE OF OPERATION 


The principle of operation will be made clearer by 
examining the process of division, z=y/w, which is 
illustrated in Fig. 1(a). Specifically we shall compute 
48/16 =3, a very simple example which will allow check- 
ing the validity of the process. The independent varia- 
bles y and w are expressed as analog voltages E, =48v 
and E,=16v which are caused to charge respectively 
two capacitors C, and C,, having equal capacitance C, 
through switches S, and S,. At a given instant ¢=0, the 
two switches are opened simultaneously causing the 
charges on C, and C,, to discharge respectively through 
two equal resistors of value R. The voltages on these 
capacitors decay exponentially as shown in Fig. 1(b), 
and their instantaneous values e, and e,, are respec- 
tively 

n= Eye #/RE (1) 


Cy ele, URS (2) 
The first of these, e,, reaches the value E, at the in- 
stant 4; when, simultaneously, e, reaches the value 


i, so that 
E, = Eye #!8C (3) 


Ee eRe (4) 
The relations among E£,, E,,, E., and E, established by 


the rundown process can be examined by dividing (4) by 
(3), which results in 
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COMPUTATION OF Ez = 


Ey /Ew 


(a) 


E, = — E,. (5) 


If E, is made equal to unity in the scale used for all 
voltages, the answer is 


E, = E,/Ew 


or, by analogy, z= y/w. 


For the numerical example given above, Fig. 1(b) 
shows that e, reaches the value ey=#,=3 v at the in- 
stant that e,=£,=1 v. The answer is #,=3 v. The 
computer requires a device to measure £, at the instant 
that é)=1 v. 

The analogy to the slide rule is apparent. To divide 
48 by 16, magnitudes corresponding to 48 and 16 are 
lined up opposite each other on the D and C scales re- 
spectively of the slide rule. The scales vary logarithmi- 
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cally from this point and the answer 48/16=3 is read 
on scale £ at the value corresponding to unity on scale 
Ce 

The computer operates on variables by repeating the 
process of solution for successive values of the independ- 
ent variable. Successive solutions can be reinserted into 
the problem to synthesize answers to special problems. 
Thus integration is performed by storing the sum of 
successive values of the function to be integrated. Dif- 
ferentiation is performed by successively subtracting 
the value of the function at one instant from its value 
after a fixed increment of time. 

Although means of solving many problems will sug- 
gest themselves, the following examples are illustrated: 


1) s=w(x/y)”. This relation is used for multiplica- 
tion, division, and powers. 

2) z=log, x, and x=a’. 

3) s=Axt+Bx'+Cxe+ +--+. This relation can be 
used to generate, or find the roots of, functions 
which can be expressed in series form. 

4) z=fotxdt and fouxdy. 

5) s=dx/dt and dx/dy. 


The exponents in 1) and 3) need not be integers but 
should be real. Quantities z, w, x, and y can be variable 
functions of time. 


INSTRUMENTATION 


In addition to the actual generation of the auxiliary 
functions, which correspond to the scales of a slide rule 
and which may obey the appropriate mathematical 
laws desired, the computer uses three special devices 
called 1) a comparator, 2) a sampling switch, and 3) a 
solution iterative switch. The operation of these de- 
vices is symbolized in Fig. 2. A practical form is de- 
tailed in the Appendix. 


1) Comparator 


This device is shown as K in Fig. 2(a). Its purpose is 
to recognize the instant tf; when a scanning voltage E&,, 
applied at point 1 is just equal to the data voltage Ez ap- 
plied at point 2. At that instant a marking signal ap- 
pears at point 3. The marking signal usually acts to 
trigger the sampling switch. 


2) Sampling Switch 

This device is shown as SS in Fig. 2(a). Switch SS 
closes for an instant At whenever EF, = Eq. The time rela- 
tions relating K and SS are shown by the graphs of Fig. 
2 


3) Solution Iterative Switch 


This is shown as S/S in Fig. 2(b). Its purposes are a) 
by closing to place the various circuits into operating 
condition by allowing suitable capacitors to be charged, 
b) to open and thereby initiate the solution by causing 
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COMPARATOR 


SAMPLING SWITCH 


E 
IE 
Eq 
Tae TIME 
SWITCH SS 
RY CONNECTS POINTS 4 AND 5 
_~JEF-2" ror snort TIME INTERVAL At 
(a) 
Sis 
ial SAR oT lk A 
|_= LINKAGE 
CAM 
DRIVING ar 
MOTOR 
(b) 
Fig. 2. 


the discharge of capacitors, c) to remain open during 
the solution, and d) to close to place the various circuits 
in condition to repeat the solution. Although the “solu- 
tion iterative switch” is symbolized for simplicity by 
a cam-operated mechanical switch, an all-electronic 
equivalent is usually preferred. 


SOLUTION OF z=w(x/y)” 


As a first illustration of the operation of this com- 
puter, suppose that it is desired to solve z=w(x/y)”, 
where m is a known constant, usually not an integer, 
and z, w, x, and y may be variables. This particular 
problem is chosen because it combines multiplication, 
division, and raising to a power, and is conveniently 
handled by the computer. 


z= w(x/y)". (1) 


The first step in the solution is to assign analog volt- 
ages E,,, E,, and E, which are proportional to the inde- 
pendent variables w, x, and y. The answer z will appear 
as an analog voltage E, proportional to g. 

These voltages are applied as shown to the circuit 
of Fig. 3(a) in which capacitors C; and C2 and resistors 
R, and R, are related by 


CiR1/C2R a | (2) 
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100 Ey 


TO MULTIPLY 
ES x Ey = ES x E 
40 x25 *10 x 100 VOLTS 


n=l 
-at 
VOLTAGE ACROSS C, =e, * Eye 


25 VOLTAGE ACROSS Coc 9 


SCALE FOR Ew AND Ey 
SCALE FOR €, AND Ey 


ss TY 


(Ey) 


SOUUTION OR Ez = 


(d) 


Fig. 3. 


Switches S; and S» are closed to charge C; and C, to 
E, and Ew, respectively. These switches are then opened 
simultaneously, causing the charges on C, and C, to 
leak through Rk; and R». The instantaneous voltages e 
and é, on C, and C, decay exponentially according to 
the time constants 


CiRi = 1/a (3) 
CoR2 = 1/na (4) 
as shown in Fig. 3(b), so that 
é; = Eye *! (5) 
and 
és = Hye? = Ey (e-%*)*. (6) 


Voltage e, reaches the value E, at time 4, so that 


= Iie = (eae ol (7) 


ee = E,/E,. (8) 


The comparator K recognizes the instant 4 (when 
¢:;=F,) and causes the sampling switch SS to close for 


a short time and immediately reopen thus trapping 
the value of e) at time 4; on capacitor Co. This value of ee 
is called £,. 

At this instant 


é2 = H,(e*4)" = E,, (9) 


Substituting (8),.(e- <4 =, 2.) an (9), 


Ig WE 
Oi i ae By ( ) 
Ey 


since E,, E,, Ey, and E, are proportional to w, x, y, and 
z, respectively, 


(10) 


w(x/y)”. (1) 


After each solution, switches S,; and S», are closed, 
by motion of the cam M, to repeat the cycle of solution 
for new values of w, x, and y. 

To illustrate the process of solution, numerical values 
are assigned in Fig. 3(b) to H, =10 v, E,=16v, E,=9 v, 
and 2=2.41. In the time required for e; to decay from 
16 v to 9 v, e. has decayed from 10 v to 2.5 v, which is 


ZZ >= 


10(9/16)?41 = 2.5 vy. 
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MULTIPLICATION 


For this computation the general equation z= w(x/y)” 
is changed to 


WX 


(11) 


by making n=C,R,/C,R,=1. w is the multiplicand, x is 
the multiplier, and z is the answer which must be 
multiplied by a constant y to give the absolute value. 
The value of y can be made equal to unity if desired. 
Figs. 3(a) and 3(c) illustrate a numerical example for 
which w=40 and x=25. 

Capacitor C; is charged through S; to a value of E, 
which is larger than #,=25 v, because the voltage e, on 
C,, must decay to 25 v. E, should be chosen also for its 
convenient use as a factor. The value /, = 100 v satisfies 
these conditions. Capacitor C, is charged to F,=40 v 
through S». Discharge of C; and C, is simultaneously 
initiated at time ¢=0 by opening S; and Sy. 

When [Fig. 3(c)] e:=E,e~**=100e7* v_ reaches 
equality with H,=25 v, at the time 4, the comparator 
K operates the sampling switch SS to trap the instan- 
taneous value of é (at t,) which is found to be equal 
to 10 v as shown below. 

im, Cpa yen | = 14,— 25790 that 


Z2y = 


Binet 25 
eva = = = 0,25, 


E, 100 


Cb 2} 


E,E, 40X25 
Pa e100 


E,= = 10Vv. 


Eye ot = 


This value of #,=10 v is then multiplied by the scale 
factor E,=100 v to give the correct absolute answer of 
1000. 


DIVISION 


Division has already been described in the introduc- 
tion. 


POWERS AND Roots OF VARIABLES 


For this problem the general equation z=w/(x/y)” is 
rewritten as = (w/y")x”, where 


(13) 


Quantity A is a scale factor which is usually made 
equal to unity or some other convenient value. In the 
following example A =1. 

As a numerical example compute z=5?-!. The struc- 
ture is shown in Fig. 3(a) for which E,=5-v and n=2.1. 
The value of E, is chosen to be greater than E,=5 v 
but not so great as to make E,?-:! excessive; a suitable 
Walie is Ho 7 vimaking £,*t= 1, =77;:=060: y.-The 
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time constants are adjusted in the ratio of C,R)/CrRy 
=n=2.1. 

Referring to Fig. 3(d) in time required for e; = Lye~*! 
to decay from FE, =7 v to E,=5 v, ¢:= Ey (e—*)" has de- 
cayed from 60 v to 29.3 v which is #, =5?=29.3 v. 

Note that it is just as practical to find the 2.1 root 
of 29.3, 


29.3121 = 5 
by sampling e, when e2=29.3 v. 


LoGs AND ANTILOGS OF VARIABLES 


For these relations we write 


z= loge x 


and 
(14a) 


x= a’. 


To evaluate z=log, x, a decaying exponential e 


; 


(14) 


=E,e—*' is started simultaneously with a decreasing 


sawtooth e. =, —kt as shown in Fig. 4. Voltages e; and 


10 - Ey 


= ON434 (fOr a =a10 


7 
= he O jfor acs « 


Ez=loggE, for a-e where K, = logge { 
TO FIND Ez = logg&y, SAMPLE @5 when e, = Ey 


TO FIND E, = a&z2, SAMPLE e, when en = E, 


Fig. 4. 


é. simultaneously reach E, and E, respectively at time 
t,. This is expressed as 


= FE, = Eye *" (15) 
and 
ai hh eae (16) 
From (15) 
€-“ti=sh./ Fi, 
so that 
a ane ae . (17) 


a y 


| 
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Substituting this value in e, (16) 
k Ee 


éo = EF, = Ey» + —log, (18) 
a y 
k k 
| = FE, + — log. Ez — — log, Ey. (19) 
| a Q 
‘Let 
| k 
Ey = log. Hy = 0: (20) 
a 
Then 
k 
éo = be = — loge Ex. (21) 
(64 
From the standard mathematical formula, 
log. « = log, € log. x (22) 


it is seen that FE, can be made equal to log, x by setting 


— = logae. (23) 
a 
For logarithms to the base 10, k/a=0.434. 
Referring to (20),. 
k 
Ew. — — log. Ey = 0 (20) 
a 
can be rewritten as 
Ew — loga€ log. Ey = Ew — logs Ey = 0 
so that 
En = log,-E,. (24) 


Therefore in Fig. 4, the initial amplitude of the saw- 
tooth (16) is made equal to the logarithm of the maxi- 
mum value of the independent variable E,. The slope k 
of the sawtooth is adjusted so that it is equal to the time 
constant (a) of the exponential multiplied by loga € (23). 
This last condition can be expressed also by saying that 
the sawtooth is equal to zero (e2=0) at the instant that 
é=1. 

It should be noted that the base of the logarithms 
can be changed by altering £, and &. It should also be 
noted that 

x= a’ 


(14a) 


can be measured equally well by sampling e at the 
instant that e.=£,. 

The numerical example in Fig. 4 illustrates the prin- 
ciple for obtaining common logarithms (@='10) for 
values of x <10. The initial value of the sawtooth at 
t=0 is equal to E,=logio 10=1. It is equal to zero 
when e,=1. For example, the value of logio 2 is seen 
to be 0.3. Other values may be checked. 
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GENERAL EVALUATION AND SOLUTION OF POWER 
SERIES 


The process previously described of raising a variable 
to a given power (z=Ax”) can be extended to generate 
many functions which can be expressed as power series 
of an independent variable, such as 


2=Ax*+ Bob+ Cae+---: (25) 


The method also affords the solution of power equa- 
tions. Thus (25) can be solved for x as well as z. The 
exponents (a, b, c, etc.) to which the variable x is raised 
need not be integers. As stated in the introduction, this 
often allows the use of more rapidly converging series, 
which require fewer terms for the required accuracy, 
than would be the case if only integers could be used. 
For example, if only the first two terms of the series 
oe eC aes 
Pin GE: Bu aes |: Wee 


toe (26) 


are taken, the series will be in error by as much as 7.7 
per cent in the range 0<x<7/2 radians. 
If, however, it is modified to 


42-9 


6.44 


sinx~ x — (27) 
the maximum error is only 0.52 per cent over the same 
range. The instrumentation for (27) will be discussed 
after the general case of (25) is demonstrated. 

Another advantage of the ability to use nonintegral 
exponents is that it greatly facilitates the experimental 
adjustment of the computer to fit a given curve. Once 
the curve is fitted, the computer can generate instan- 
taneous values of the function for instantaneous values 
of the independent variable or, vice versa, often solve 
for the independent variable for instantaneous values of 
the function. The ability to solve such series is believed 
to be much easier than with earlier types of computers. 


EVALUATION AND SOLUTION OF 
2=Axt+Bx'+Cxue+ --- 


Fig. 5 shows a block diagram of the arrangement to 
evaluate z in this case. 

When the switch S; is closed, C, is charged to the 
value FE, and C,, G, C., etc. are charged respectively 
to 1Ra, tk, 1R., etc. When the switch S, is opened C, 
discharges through R, so that the voltage e; across C, is 


€j = Eye"! Belz, (28) 
Similarly the voltages across Cz, Cy, C., ete. 

Cg = 1Rge— tl Bala 

Cy = tRye /RtCe 


Since RaCa, R,C;, R.C,, etc. are in series the voltage es, 
across the chain, is 


€o = iRgem"/RaCa + GRye-t/R0Cd + iR,et/RCe 4... (29) 


Let 


TRAC = ond Or — bRiC,y = Ghee 


(30) 


Substituting these values in (29) 
eo = iRge ttlRaCe 4 GRyePt/RaCe + GRye-ct/ReCe +... 
= iR,(ert/RxCz)@ 4. {Ry (er "/ReCz) 


+ 4ReetlPo0e)¢ + (31) 


The voltage e:, across C, reaches the value /, at the 
time f, or, 
= Vs — Eien BxGe 
which gives 


eT h/RzCz = pp 0 One (32) 


At this instant (é,) the comparator K closes the sam- 
pling switch SS momentarily thus trapping on capacitor 
C the value of e at time 4. 

At time ¢,;, @2 has reached the value E, so that sub- 
stituting ¢=¢; and (32) in the value for eg given in (31), 


(EONS Ta NG 
GS ONE WIR =) + iRs/ ) 
Ey Ey 


NS 
. Ey, 
which can be rewritten as 
. iR, pat iR, omy ik, pene (34) 
63 — 265— z" as, x ary x tans 3 


This equation is seen to be the electrical analog 
to 


g= Axt+ Boot Caoe+-::- (25) 
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where 

iR, 
A=—) B=—) 
Ee 
UR» 
gee 


(35) 


The structure described above evaluates zg in (25) 
when x is known. It is obvious that, by interchanging 
the comparator K and the sampling switch SS, it is 
possible to solve for x when 2 is given. 

Negative coefficients can be introduced by charging 
the proper capacitor with a battery having the opposite 
polarity, by using an inductor-resistor combination and 
placing a secondary winding on the inductor, which 
secondary can be poled to either polarity, or by special 
artifice such as is shown in Fig. 6 and described in the 
next section. 


EXAMPLE 


To COMPUTE 2=SIN x, AND * =ARCSIN 2 


As discussed above, 


2-9/6,44, (27) 


sin «~~ x — 


The circuit for this computation is shown in Fig. 6(a). 
When the switch S at the upper left is closed, battery 
V charges C; and C; to voltages 7k; and 7k, respectively. 
Opening the switch will cause the voltages e, across 
R, to decay according to 


Ge = tR ye / R101 — dees, 


Likewise on opening switch S, the voltage across 
C2, which is equal to 7R3, divides between R, and R; in 
such a manner that the portion across R: is equal to 
—iR2R3/(Ro+R3). Note that the current flows from 
R, into Cs, hence the negative sign. 

As shown in Fig. 6(b), the voltage e, across Ro, will 
decay according to 


RoR; 
.o— 
Ro + Rs 


& = — et] et R3)C2 


so that the sum e;+e, is 
RsR3 


es = ey + ey = iRye H/RCr — § —— 
Ro + R; 


e—t/ (RatR3)C2 


(36) 


At time t=, e, decays to the value E, which is the 
analog of x so that it must start from a value £,=iR, 
which is higher than the maximum value of x. Since this 
was taken aS. %max=7/2, H,=2. Veils a. conveniens 
value. 

Writing (27) in analog form with 


sin Z, = E, — E,?-9/6.44 
= 2-78 — (Qe—%"1)2-9/6, 44 


os Dee th ay 1.16672-9441, 


(37) 
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sin Ey 


IR; = 100 VOLTS MAX 
IRz= 58 (Ro+ R3)/Re VOLTS MAX 
R, C) A(R2+R3) C2 = 2.9 


100 V 
OR 
2 RADIANS 


\ 


=5 8 V 
INSTRUMENTATION TO COMPUTE 
50 sin x ~ 50x — 50x?-%/6-44 FOR 0 <x < 1/2 
~ 100e-%% — 58e-2.9%4, 
(b) 
Fig. 6. 


Because the comparator described in the Appendix 
operates best with higher voltages than those in the 
range 0-1 v, the scale of (37) is changed by multiply- 
ing the whole equation by some convenient factor such 
as 50. This equation then becomes 


50 sin Ey = 50(2e-**!) — 50(1.1667-?-9#4) 
1O0QE= 8 = 5 8e=2 24 


l| 


(38) 


To implement this equation, the constants in Fig. 
6(a) are adjusted as follows: 


1) iR; and i1R2R3/(R2+R3) are made respectively 
equal to 100 v and 58 v. 

2) The time constant 1/(R2+R3)C2 is made equal to 
2.97 Ri Cy, 


Then when the switch is opened, the voltage e; on 
decays according to e¢;=100e~*’ and reaches the value 
EF, at time t, so that 


E, = 100e>2". 


At the same instant (t=), e3 has reached the value 
given in (38). The comparator A causes the sampling 
switch SS to close momentarily causing the value 50 
sin EL, to be trapped on C>. 

To compute x=arcsin z, the comparator K and the 
sampling switch are interchanged. 
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INTEGRATION 


Integration with respect to time is performed by add- 
ing successive values of the independent variable which 
are spaced by equal increments of time. An extension 
of the method which allows integration with respect 
to another variable is discussed in a later section. 

Function e(t) to be integrated with respect to time 
from ¢, to fy, shown in Fig. 7(b), next page, is impressed 
across terminals 1, 2 of the circuit of Fig. 7(a). A source 
of periodic voltage preferably linear but expressed for 
convenience here as /o sin wt is produced by trans- 
former T across terminals 4, 5. Voltage Eo sin wi, in Fig. 
7 (b), is made greater in amplitude than the volt- 
age across terminals 2, 3 of Figs. 7(a) or 7(c), 


f ” di 


1 


and w/2m is made higher than the highest frequency in 
e(t). Assume for simplicity that the initial charge on C 
is zero. 

As Ey sin wt increases, it reaches a value at time 
when it is equal to e(t;) plus the voltage on C (which is 
assumed for convenience to be initially zero). At that 
instant Eo sin wt,=e(t). The comparator K closes the 
sampling switch SS for a short time which is long 
enough to charge C to the value e(t,;). The voltage e(t) 
is trapped on C by the opening of the sampling switch. 
During the next cycle, Ho sin wt reaches a value, at 
time f, when it is just equal to the voltage e(t,) stored 
on the capacitor plus the new value of e(¢) at the time 
t so that 


Eo sin wl, = e(t1) + e(te). 


The comparator again closes the sampling switch long 
enough to charge C to this value. The switch then opens, 
trapping this value on the capacitor. It is obvious that 
after cycles, the voltage on C will be 


Eo sin wt, = e(t1) + e(te) + +++ + e(tn). (39) 


Since the intervals At between samples can be made 
very nearly equal we can write 


n 


E,At = >. e(i)At 


1 


and if the samples are made very frequently 


n t 


E,At = >. e(t)At~ e(4)dl. 


1 t=0 


(40) 


While the intervals At are not strictly equal to each 
other they can be made very nearly so by making Eo sin 
wt rise very steeply. 

The capacitor C can even have appreciable leakage, 
as the charge on it is restored after each cycle and the 
leakage from C can be compensated by charging C each 
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le ee 


Eg SINwt 


e(t) 


Jecrrer 


Bae Rb ~o¥pc 7 | 
Rp hs 
WHERE At = 27w 
(c) 
Fig. 7. 


time from a voltage just great enough to compensate 
for the leakage. Thus if C is shunted by R, the voltage 
on it would be reduced to 


E.e Atle 


after the interval At. Since Af is the interval between 
cycles its value is well known. The voltage applied to 
the capacitor is therefore increased by a factor M such 
that 


Me-Stlke = 1 
as shown in Fig. 7(c). 


DIFFERENTIATION 


The circuit for differentiation is shown in Fig. 8(a) 
and its operation is shown in the waveform diagram 
of Fig. 8(b). The circuit samples the function to be 
differentiated e(t) at regular intervals. The value of the 
function at these times is written e(t,), e(tr_1), etc. The 
value of the function at time (f,_1) is stored so as to 
subtract it from its value at the time ¢t,. The difference 
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e(t)-e(th) 
O 
4 AVERAGE OF e(t)— e(tn) 
' 
ce 
(b) 
Fig. 8 


e(t,) —e(tn-1) between intervals is proportional to the 
derivative de(t) /dt. 

Referring now to Fig. 8(a), whenever Eo sin wt =e(t) 
at times hy, fo, tn—1, tn, the comparator K closes the sam- 
pling switch SS for a short interval of time during which 
capacitor C is charged to the value of e(t,) as shown by 
curve FE, of Fig. 8(b). At this instant the voltage be- 
tween terminals 1—3 is equal to zero and assumes a saw- 
tooth form between instants ft, fo, etc., the instantaneous 
amplitude of the sawtooth being equal to the difference 
between the instantaneous value of the function and its 
value at the time of last sampling. Since the intervals 
of time between samples are very short, the function 
can be assumed to be linear between these samples 
and the average value of the sawtooth is therefore 
4[E(t,) —E(tn1) |, which is nearly equal to 4(de/dt). 
The filter F in Fig. 8(a) does the averaging. 


INTEGRATION AND DIFFERENTIATION WITH RESPECT 
TO A VARIABLE OTHER THAN TIME 


Fig. 9(a) and 9(b) shows block diagrams whereby 
y can be integrated, or differentiated, with respect to 
x. The computer makes use of components which have 
been described above: 


1) Integration z= fydx [Fig. 9(a) | 


Step 1—~x is differentiated to dx/dt. 

Step 2—y and dx/dt are multiplied to obtain 
ydx / dt. 

Step 3—y(dx/dt) is integrated with respect to time 
to obtain the desired result fydx. 
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day, 


——________________ x 


(b) 
Fig. 9, 


2) Differentiation z=dy/dx [Fig. 9(b) | 
Step 1—-y is differentiated to dy/dt. 
Step 2—-~x is differentiated to dx/dt. 
Step 3—dy/dt is divided by dx/dt to give dy/dx, 
the desired derivative. 


OTHER APPLICATIONS 


Many other applications of this sampling parametric 
method of calculation will suggest themselves. Thus it is 
possible 


1) To determine sin (nx+0@) when sin x is known by 
generating the two waves and sampling one at the 
instant that the other reaches the known value. 

2) To set up special series such as z=Aa*+Bb** 
+Cc* and to solve these series explicitly and 
implicitly. 

3) To add voltages (E,, E, for example) relative to 
a common point and get an answer relative to 
this common point. 

4) To multiply by a constant n. 


APPENDIX | 


A computer was constructed to test the theory out- 
lined above and determine possible accuracy. This com- 
puter solves the equation 


Ey Eye 
E3|—| — £4} — 
1 Eo 
where any of the voltages, or combination of voltages, 
may be treated as independent variables and the ex- 
ponents a and b can be any real numbers. 
Fig. 10 shows the manner in which the equation is 
related to the exponential discharge curves of the com- 


puter. Fig. 11 shows the block diagram of the computer 
which repeats its operation 60 times per second. 


Es = (41) 
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t 
ep) = Eoemne 


t 
e pee Ee 
pore = me RC | 
e3 =E3€ RC 2 t 


_bt 
C4 =Eqge RC 


Q@o = E 
2 f TIME Se 


Es =e3 -e4 E4 
AT MOMENT OF SAMPLING b | 
E| THESE TWO 


Z ial 
Ee = &3:(—-) Seas AT MOMENT 
Str Es ee) ear OF SAMPLING 


Fig. 10. 


e4 
GENERATOR 


Fig. 11. 


Three identical circuits are used to generate the ex- 
ponential waveforms. The de voltages /2, E3, and Fy 
determine the starting voltage of the eé, e3, and e4 ex- 
ponential waveforms and each decay rate is determined 
by the RC time constant of its associated circuit. The 
exponents a and 6 are related to the time constants by: 


RoC, ‘time constant of é» 

a= = — (42) 
R3C3  ~+time constant of es; 
RC. ~+time constant of e2 

6= =— ‘ (43) 
RzC, time constant of es; 


The several decaying voltages é2, e3, e4, and é; are 
related to each other by the following equations 


€y = Eoe~/®C; whence e2/E, = et/RC (44) 
€3 = Eze tt/RC = E3[e—#/R¢ |a me E3|e2/ E> |* (45) 
ey = Ege PtlRC = Ey [e-RC]> = Eylen/Es]® (46) 
és = €3 — €, = E3[e2/Eo|* — E,leo/ Ea? (47) 


At the time ¢,, that e: voltage falls to the value equal 
to /,, (Fig. 10) the comparator (Fig. 11) operates the 
diode switch which then charges the capacitor A to a 
value (e3—e,) = 5 at that instant. 
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Fig. 


The repeater reproduces the dc voltage on capacitor 
A at the output terminals designated E;. This voltage is 


; Ey, a By b 
Es — (es — é€4 at time ty) = E3 =) co Ey, =), (41) 


2 2 


Referring to Fig. 12, the clamp waveform, which is 
shown at the bottom of each generator, is supplied 
from a large 60-cycle sine wave which is amplified and 
limited to give a square wave. This sine wave determines 
the rate at which the computation is repeated. 

During the positive half of the clamp wave, the cir- 
cuit for the e.RC time is clamped to the voltage Fo. 
Because of the clamp circuits chosen, the voltages Fi, 
Fy, Es, Ey, and &; are all negative with respect to ground. 
When the negative half of the clamp wave occurs, the 
clamp tube stops conducting and RC circuit begins to 
rise toward zero at a rate determined by the RC time 
constant of the values chosen for the resistor and capaci- 
tor. The cathode follower provides a low-impedance 
output for the generator, but introduces a dc differ- 
ence between its grid and cathode. In order that it be 
independently accurate in dc level, the stationary end 
of the RC circuit is returned to a point below ground 


so that when the RC circuit is fully discharged, the 
cathode follower gives zero de volts output. 

The output of the e2 generator goes to the comparator, 
where it is compared with the dc voltage /;. As the 
voltage e; becomes less negative, it finally reaches the 
value of £; volts, and causes plate current to flow in 
the 1-meg plate resistor of the comparator. The sudden 
negative swing of the plate is amplified through the 
12AU7 amplifier and triggers the blocking oscillator. 

The blocking oscillator causes conduction through 
the diode switch, and connects e3 and e, to the capacitor 
A for a period of 10 microseconds. The dc voltage across 
A represents the answer, but both sides of the capacitor 
have a common potential which has both ac and de 
components relative to ground. 

The method used for repeating the de voltage differ- 
ence with one side grounded, is to reverse with respect 
to ground, and take the average of the voltage e3; and 
— és which is $(e3—e4) or 35. Both voltages e3 and es are 
repeated through cathode followers to avoid loading 
the charge on capacitor A. The voltage e, is reversed 
by using feedback across the two 1-megohm resistors at 
the left in Fig. 3. The midpoint of the two resistors is 
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kept at ground potential by applying to one end of the 
resistors a feedback signal which is equal and opposite 
to es. The two 1-megohm resistors to the right are used 
to obtain the average between e3; and —e, which gives 
the desired output. 

In operating the equipment, the value of A, must 
always be less than £2, and the minimum value of F; 
that can be used is determined by the minimum value 
that e: reaches at the end of its exponential decay cycle. 
If the RC product is chosen equal to 0.2 of the available 
decay time, then values down to 0.006 of full scale are 
possible. 

The accuracy obtained with the present equipment 
is between 1 and 2 divisions over most of the meter 
scale of 100 divisions, except at the end of scale above 95 
where some problems produce errors of 4 divisions. 

By reconnecting the elements of the computer as 
shown in Fig. 13, the output obtained is the time integral 
of the input voltage. The potentiometer is normally set 
near the top of its range. If E, is zero, the comparator 
will close the switch when the supply voltage reaches 
the voltage e,, thus maintaining whatever charge is in C. 

If the voltage F, is greater than zero, the switch 
will close at the time that the supply voltage reaches 
EF, volts above e.. Thus e.+; volts are placed across C. 
Because the change that occurs in e, is always equal to 
F,, it follows that e, is the integral of /, for the number 
of samples taken, or if the repetition frequency is con- 
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stant, for the amount of time elasped. 

Because some leakage, R, exists, the supply to the 
switch is made greater than the supply to the sample, 
and thus keeps the charge from draining off from C. 
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A Gray Code Counter* 


A binary counter may be programmed to 
count according to the reflected (Gray) 
code.! A suitable block diagram for four 
stages is shown in Fig. 1. The counts are fed 
into flip-flop Fy whose two differentiated 
outputs are utilized to trigger the counter 
proper consisting of Fi, Fo, F3, and F,. 

Each flip-flop consists of two tubes, one 
of them conducting at a time, and the fol- 


pulses in 


“one.” In Fig. 1, the next trigger from Fo will 
flip Fi, consequently changing the state of 
the counter to 1110. This corresponds to 5S. 
The next trigger from Fp is transmitted to 
all gates. However, only gate 1 is open, and, 
consequently, F2 is triggered. The count is 
now 1010, corresponding to 6. Fig. 2 shows a 
practical diagram of a reflected code binary 
counter containing six stages. 

Throughout a regular counting sequence, 


Fig. 1. 


where it appeared last. For instance, starting 
with 1010, corresponding to 6, the last trig- 
ger appeared at the bus leading to all gates. 
Another trigger applied to the same point 
will find gate 1 still open, and return F, to its 
previous state 1. The count is now 1110, 
corresponding to 5. The next regular count 
flips F, and produces 0110, corresponding 
to 4. 

This reversing property of the counter 
enables one to build a forward backward 
counter, and Fig. 3 shows the logical arrange- 
ment for the input stage of such a counter. 
A separate input is provided for additive and 
subtractive counts. In the state shown, an 
input pulse at the “add” terminal first flops 
Fy. Then it produces, after the delay neces- 
sary for Fo to reach the steady state, an out- 
put at the gate bus. A subsequent pulse at 
the “subtract” terminal returns Fo to its 
previous state and again produces an output 
at the gate bus, as required. 


TRIGGER IN 


1050 


lowing conventions will be followed through- 
out this discussion in describing their states. 
A tube at cutoff represents a zero, and a con- 
ducting tube a one; in flip-flop Fo, a trigger 
will be generated by the tube which is chang- 
ing from zero to one; finally, it will be as- 
sumed that the lower half of F,, Fs, F3, and 
F, represents the binary number contained 
in the counter. In our case, this number is 
0110 which corresponds to 4, in the reflected 
code. 

The operation of the counter will now be 
described. One trigger produced by F» is con- 
nected directly to Fi, and the other one toa 
bus bar feeding all the other counter stages 
through “and” gates. The proper gating ac- 
tion is provided by the flip-flops F,, Fs, and 
F;. The “and” gates will transmit the trigger. 
from Fo only if all their other inputs are 
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the trigger will appear alternatively at the 
input of Fi, or at the bus leading to all gates. 
The counting direction will be reversed by a 
single repetition of the trigger at that point 
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A Method for Obtaining Complete 
Digital Coding Chains* 

Given K different symbols (e.g., the 
digits 0, 1, 2,- +--+, K—1), one can always 
write a closed loop in which each symbol ap- 
pears K” times, arranging the symbols in 
such order that every subsequence of 
contiguous symbols is unique. It is then pos- 
sible to read K” different code combinations 
or numbers, corresponding to the K” differ- 
ent reading stations for the set of 2 symbols. 
This assertion can be proved by showing 
that the permutation scheme described be- 
low is exhaustive. A forthcoming SCEL re- 
port! will contain such proof. 

Johnson? has described sequences of this 
kind for binary digits. He calls them “binary 
chains” and attributes their discovery to 
Baudot. Johnson shows methods which as- 
sist in the discovery of “complete” binary 
chains, (¢.e., sequences which contain all 
2” n-bit combinations), but his rules do not 
work predictably in all cases. By contrast, a 
method in use at SCEL gives complete 
chains for any radix and any . This method 
is as follows: 


1) Assign to K symbols an arbitrary 
“order of preference.” 

2) Imagine a “register” consisting of 1 
boxes. Fill this register with the last preferred 
symbol written n times. 

3) Shift the contents of the register one 
place in either direction. One symbol will 
overflow out of the register and a box will 
become vacant at the opposite end. Write 
the first preferred symbol in the vacant box. 

4) Continue to shift the contents of the 
register in the same direction. Continue to 
write the first preferred symbol in the vacant 
box after every one-place shift; then com- 
pare the new combination of symbols in the 
register with all previous combinations. 
Eventually, a combination will repeat. When 
this occurs, remove the first preferred sym- 
bol just written and replace it with the sec- 
ond preferred symbol. If the second preferred 
symbol also gives a repeat, use the third, etc., 
but do not violate the order of preference. 

The above process gives all possible per- 
mutations of the symbols before it is ex- 
hausted. It follows that the sequence of 
symbols occupying any one box forms a 
complete coding chain when closed upon it- 
self. 

For example, let K be three, let 0, 1, and 
2 be the three symbols in order of preference, 
and let be two. Shifting the register con- 
tents to the left, we get in succession: 


2 
20 
00 
01 
10 
02 
zal 
11 
12. 
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In the endless chain made by closing 
220010211 on itself, all nine two-place 
ternary numbers can be obtained by reading 
two adjacent digits. 

For binary chains, one may start with all 
0’s and try 1’s first. For large n, it is con- 
venient to consider the successive n-digit 
combinations as conventional binary num- 
bers in order to use their decimal equivalents 
as a shorthand for generating the sequences. 

The rule for left-hand shift is, then, to 
start with zero and to perform repeatedly 
the following process: 

Double modulo 2”, and add unity to 
make an odd number. If the resulting odd 
number is a repeat, subtract unity from it, 
to make it even. Either the odd number or 
the even will be new, until both unity and 
zero repeat after 2” different numbers. Thus, 
giving 2 the value four, we get 0, 1, 3, 7, 15, 
14, 13, 11, 6, 12, 9, 2, 5, 10, 4, 8. By writing 
an 0 for every even number and a 1 for every 
odd, we get 0111101100101000 which is a 
complete binary chain of 16 places. 

In general, the process described defines 
2K! different chains for all K and all n=2, 
corresponding to all possible orders of prefer- 
ence for K symbols and to both directions of 
shift. (Known exceptions are » =2 or 3 when 
K=2, for which unique chains occur.) 
Nevertheless, there are other complete 
chains not obtainable in this way. Thus, 
neither 1111010000101100 nor 221001120 is 
obtainable by these rules. 

In general, the binary chains described 
herein do not coincide with those of Johnson. 

Neither Johnson nor the authors have 
yet found a simple means of translating to 
conventional number codes. 

One possible use for chain codes is in in- 
ventory or catalog systems employing 
cyclicly interrogated digital memories such 
as magnetic drums. With a complete binary 
chain, only 2” bits are needed to record 2” 
catalog numbers or addresses, while still re- 
cording the entire address for every location 
—a possible saving of (7 —1)2” bits. 

Here, the difficulty of translating to nat- 
ural numbers is no handicap. Other applica- 
tions can be visualized in data transmission, 
in certain special analog-to-digital conver- 
sions, and in some digital data processing. 
Where “incomplete” chains can be devised 
with superior capability for translation, etc., 
they may be preferable in all the latter cases. 
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A Survey of the Characteristics of 
Currently Used Bistable Multi- 
vibrators* 


A survey of the characteristics of cur- 
rently used bistable multivibrators has re- 
cently been carried out in the College of 
Engineering at the University of Wisconsin, 
with the aim of collecting a body of data and 
comparisons which might assist those faced 
with the selection or design of such circuits. 

The analysis of a bistable multivibrator 
or even the synthesis of this circuit in ac- 
cordance with dc or steady-state considera- 
tions is a relatively simple matter. The 
transient analysis, however, which en- 
deavors to determine the resolving time or 
maximum repetition rate, is much more in- 
volved. For ease of discussion, the resolving 
time may be separated into: 1) the settling 
time, which is dependent on various RC 
time constants and is a measure of the time 
which must elapse from the completion of 
flip-flop transition until the circuit can be 
successfully triggered again, and 2) the 
transition time which involves all the circuit 
elements including the tube characteristics 
and their associated nonlinear effects. The 
transition time is particularly difficult to 
calculate accurately, and in most cases is a 
byproduct of the de circuit design. 

In an effort to overcome these difficulties, 
an attempt has been made to determine ex- 
perimentally various circuit characteristics 
by compiling data on a number of flip-flop 
circuits which have been used commercially 
in sizeable quantities by various manufac- 
turers and universities. A standard ques- 
tionaire was sent to 143 organizations in the 
United States, Canada, and England. In re- 
sponse to this, 88 replies were received, 58 of 
which comprised useful data. From over 75 
circuits, 31 circuits were picked as repre- 
sentative, or chosen to be tested for purposes 
of comparison. 

Whereas it is probably true that no two 
flip-flop applications can be satisfied by the 
same optimum circuit, it is hoped that the 
circuits included in this report will be di- 
verse enough to satisfy specifications ap- 
proximating those desired. Modifications 
made to these circuits in accordance with 
well-known theory should produce a satis- 
factory circuit with a minimum of effort. 

A report describing the results of the 
survey has been published by the Engineer- 
ing Experiment Station of the College of 
Engineering, University of Wisconsin. Re- 
quests for single copies of this report may be 
sent to the undersigned. 

CuHarLes H. Davipson 
Dept. of Elec. Eng. 
University of Wisconsin 
Madison, Wis. 
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Unit-Distance Binary-Decimal 
Code Translators* 


In a recent discussion of unit-distance 
binary-decimal codes,! Tompkins remarked 
that codes with simpler rules for translation 
to conventional constant weight codes might 
be found. Gilbert has enumerated all classes 
of unit-distance 4-bit binary-decimal codes.? 
From these it can be shown that there is 
only one distinct class of unit-distance 4-bit 
binary-decimal codes for two-track com- 
mutation in analog-to-digital converters. 
No unit-distance binary decimal code using 
four brushes in one track may be devised.® 
At least one unit-distance code using five 
brushes in one track is possible.*4 

First, consider a unit-distance code used 
ona particular angular position encoder and 
its associated translator which produces a 
desired output code. Permutation of this 
code corresponds to relabeling the brushes 
of the encoder in a different order or shifting 
the positions of tracks of each decimal digit 
radially, among themselves, without chang- 
ing the labels on individual brushes. Ap- 
propriate relabeling of translator leads 
leaves its elements unchanged. Comple- 
mentation of a code corresponds to inter- 
changing conducting and nonconducting 
segments of selected tracks. This does not 
change the required total number of ele- 
ments in the translator. All codes which may 
be generated by permuting and/or comple- 
menting a given code have translators of 
identical complexity and will be considered 
to bea single class of codes. Only the number 
of different classes of codes need be deter- 
mined. One member of each class of code is 
then investigated in detail, to determine the 
number of translator elements for all codes 
of that class. 

It should be noted that horizontal or 
vertical translation, rotation, reflection 
about a horizontal or vertical centerline or a 
left or right diagonal produces a map of a 
new code which in the same class as the origi- 
nal code and does not have an essentially 
different translator. Also, reflections about 
minor axes between the first and second 
square and between the third and fourth 
squares in either the vertical or horizontal 
direction do not produce a code of a new 
class. Some simple algebraic transformations 
of a code, such as permuting a and d, pro- 
duce complex geometric changes in the oc- 
cupied squares of a map. 

In Tompkins’ Fig. 2, the first and third 
maps in the left column are not related to 
each other by a single translation or a single 
rotation as are the other maps of the figure. 
Using McCluskey’s method for detecting 
symmetry,° it can be shown that these codes 
are related by: 


* Received by the PGEC, February 20, 1957; re- 
vised manuscript received, March 14, 1957. 
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ao = C3 
bo = bs 
Co = 43 
do = ds’. 


Using either of these maps, but not both, 
decimal zero may be placed in any of the 10 
squares containing an x and the sequence 
may proceed in either direction. Evaluating 
the 20 different translators completes the 
analysis. Placing decimal zero in the upper 
left hand corner of the map and decimal 1 
to the right of it permits the code words for 
0, 1, 2, 3, 8, and 9 to differ only slightly from 
their usual binary values. That is, the en- 
coder map fits the map of the translator out- 
put code well in a number of places. 

Fig. 1, presented below, shows maps of 


Fig, 1. 


two of the twenty maps using two different 
sets of coordinate labels. 

The figure also shows values to illustrate 
the relationship of these maps to Tompkins’ 
maps. However, this is not a rigorous proof 
that these choices require a simpler trans- 
lator than the other 18 possibilities. 

The logical expressions for translation to 
1, 2, 4, 8 output code are given for both con- 
ventional numbering (requiring logical de- 
tenting) and for reflected decimal numbering 
notation, 


June 


LoGIcaAL EXPRESSION FOR 
ENCODER TRANSLATORS 


a, b, c, d=Inputs from encoder brushes 
1, 2,4, 8=Outputs from translator 
e=1 output from translator for next 


most significant decimal digit 
Ko—LocGicaL DETENTING 
2=ad+b'c 

1 = ab’+(a+b)(c'd+cd’) 


Siesacide 
4 = d(c + d) 


Ko—REFLECTED NUMBERING 
8 = (ad'e’ + a'b’e)c’ 
4=d(c+ de’ + (a+ d')ce 
2=ad+0'c 
1 = [ab’ + (a + b’)(c'd + cd’) Je’ 
+ [a’b + (ab) (cd! + cd) |e 


K;—LocicaLt DETENTING 

yc vole)’ 2=(b6+c)d 

4=ab+cd' 1=ad'+[a'b+(a' +d)c' ld 
K;—REFLECTED NUMBERING 

8 = ab’e’ + a’cle 

4= (ab+cd')e’ + (a + d’)ce 

(b+ c)d 

= [ad + [a'b + (a + d)c'|de’ 

+ [a’d’ + b’c + ad(b’ +c) Je 


m bo 
ll 
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Translators for either of the above re- 
flected numbering schemes are more com- 
plex than translators for 4-bit reflected bi- 
nary-decimal codes.’ Two classes of codes 
permit conversion from reflected to conven- 
tional decimal number notation by comple- 
menting the output of one of four brushes in 
separate tracks where all pairs of words read 
from adjacent sectors are unit distance. The 
expression given in Tompkins’ paper may be 
used with logical detenting between decades, 
but not reflected numbering. Glixon’ has de- 
scribed a procedure using conventional num- 
bering notation, where each decimal digit is 
assigned one of a pair of words (which differ 
in one bit) depending on whether the next 
more significant digit of a complete number 
is even or odd. Using this procedure with the 
unit-distance two-track code K, Tompkins’ 
translator expressions apply. Indeed, Glix- 
on’s procedure extends the number of unit- 
distance codes that should be investigated 
for translator complexity because any unit- 
distance sequence of length ten may be used 
whether it closes or not. 
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Negative Base Number Systems* 


It is commonly stated that any positive 
integer 7 greater than unity can be used 
as the base of a number system. Thus, a 
number is expressed by a polarity sign 
(plus or minus) together with a shorthand 
notation for (a@n7"tdnar™ !+ - ++ +a +ao 
raf) ---), viz, EOn@n+4 °° G0'G1-*, 
the base being previously agreed upon, and 
the a; are integers 0 Sa;<r. However, it ap- 
pears that negative integers may also be 
used as bases. For example, 101101 would 
represent in a “negative two” system the 
decimal equivalent —35, derived from (—2)é 
egy — 2) 45-8 (2) P08. C2) l= 
—324+0—84+4+0+1 = —35. Similarly, 10011 
represents +15, and 1.11 represents +0.75. 
Notice that in a “negative base” system a 
separate sign indication is not required, the 
sign being self-contained. Table I of “nega- 
tive two” equivalents follows. 

It will be noted that any positive or nega- 
tive number can be represented. Whether a 
number in this notation is positive or nega- 
tive can be found by inspection, positive in- 
tegers having an odd number of figures, 
negative integers an even number of figures. 
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Correspondence 
TABLE I 
Negative Number Positive 
0 zero 0 
11 one 1 
10 two 110 
1101 three 111 
1100 four 100 
1111 five 101 
1110 six 11010 
1001 seven 11011 
1000 eight 11000 
1011 nine 11001 
1010 ten 11110 
110101 eleven 11111 
110100 twelve 11100 
110111 thirteen 11101 
110010 fourteen 10010 
110101 fifteen 10011 
110000 sixteen 10000 


The necessary addition and multiplica- 
tion tables are given below. 


Addition Multiplication 
0 1 0 1 

0 0 1 0 

1 1 110 1 


[Subtraction rules can be inferred 
from the above, but a straightforward 
procedure is simply to multiply by 
minus one (z.e., 11) and add. | 


Note the “double carry” that is required 
in adding 1 to 1. In many cases of practical 
addition or multiplication, this results in a 
nonending sequence of carries to the left; 
but once the carry passes out of field, the 
remaining representation is correct. For ex- 
ample, adding minus one to minus one: 


11 
+ 11 
sense i : 
A gh carries 
-++ 00010 sum 


No attempt is made to assess the practi- 
cal value of a “negative base” system. How- 
ever, its advantage of a set of rules for arith- 
metic free from dependence upon a prefixed 
polarity sign appears worth examining. 

Louts B. WADEL 
Chance Vought Aircraft, Inc. 
Dallas, Texas 


A Fast Circulating Memory* 


The recent development of diffused base 
transistors with alpha cutoff frequencies 
above 200 megacycles makes it possible to 
build very fast logical circuits. By making 
use of some of these circuits in the regenera- 
tion logic of a circulating memory, a very 
short access time can be realized. One ap- 
plication of such a fast small capacity mem- 
ory might be to replace many of the registers 
of a parallel computer. To do this, a number 
of delay lines equal to the number of digits 


* Received by the PGEC, March 14, 1957. 


123 


in the computer word would be used and 
would be circulated in synchronism by a 
common clock signal. The number of digits 
in any one delay line would be equal to the 
number of registers required. 

As an illustration of the speeds attain- 
able a small memory was built using an 
ordinary coaxial cable as the delay line and 
experimental Western Electric transistors in 
the regeneration logic. Eight digits were 
stored in the line and were circulated by 
means of the regeneration logic shown in 
Fig. 1. The regeneration is accomplished by 


GATES TO ‘I WHEN INPUT IS "J" 
“FOR CLOCK*| 


CLOCK IN fs) 


“O" FOR CLOCK=1 


GATES TO "0" WHEN INPUT IS “0” 
se 


COAX CONTAINS 


om 7 DIGITS 


Fig. 1—Logical diagram of delay line memory. 


using a flip-flop, at the input to the cable, 
which is set on each clock cycle to agree with 
the state of the output end of the cable. The 
timing of the various points of the circuit is 
shown in Fig. 2, which is drawn to represent 


| | | | 
CLOCK | | | | 
° | 1 


| | 
= 36.3 mys + | 
| | | 


OUTPUT Ip 
OF CABLE “ I™s a ~ 
° | 
Fig. 2—Timing diagram. 


the circuit operation near the maximum 
speed so that a clock period is equal to 
twice the time required to change the state 
of the flip-flop. Note that a nonreturn-to- 
zero system is used which permits the storage 
of a larger number of digits in the same 
length of cable. With the choice of clock 
period above, the output end of the cable 
remains in the same state throughout the 
length of the flip-flop gate signal (4 clock 
cycle). This is required for the maximum 
speed of gating. The length of the cable is 
chosen so that the delay is 7} clock cycles; 
the other 4 cycle is the time required to set 
the flip-flop. 

For the test circuit the cable was 212 
feet of RG62/U cable which has a delay of 
276 mys. For this delay to contain 74 clock 
cycles, the clock frequency was made equal 
to 27.2 megacycles. The transistors were 
Western Electric type GA-53233 diffused- 
base p-n-p. The memory so produced stored 
eight digits with a maximum access time of 
295 mus. Since the cable delay is only about 
1.3 mys per foot, the exact length of the 
cable for a given clock frequency is not criti- 
cal. An error of a foot or two in the above 
cable probably would not be noticeable in 
the circuit performance. 

Fig. 3 is a circuit diagram of the test 
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Fig. 3—Circuit diagram of test memory. 


memory. ‘wo transistors are required to 
generate a push-pull clock signal. In a large- 
scale memory this part of the circuit would 
be constructed only once for the whole 
memory. The rest of the circuit would be re- 
quired for each cable used in the memory. 
The write-in circuit which was used in this 
test memory was made as simple as possible 
to minimize the number of transistors re- 
quired. It therefore merely sets the flip-flop 
at some random time and is used only to 
produce typical storage patterns. 

A typical waveform at the input end of 
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the cable, point A of Fig. 3, is shown in Fig. 
4. The clock frequency is only 17 mec and 
there are only five digits stored in this 
figure, but the waveform is typical for the 
operation near the upper limit of speed. It 
has the same shape as did the same signal 
at 27.2 mc originally. The lower frequency of 
operation was necessary, because one of the 
original transistors drifted in its character- 
istics so as to be unusable before this wave- 
form was recorded and a replacement was 
not available because of a shortage of these 
transistors. When a transistor of lower 
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Fig. 5—Memory with READ and WRITE logic. 


alpha than the actual design limit of the 
circuits was substituted, the maximum speed 
of operation was reduced. 

In order to get information into and out 
of an actual memory, some additional logic 
is required. The logical diagram of one way 
to do this is shown in Fig. 5. The READ and 
WRITE signals coincide with the particular 
clock pulses to affect the position desired in 
the cable. The techniques for generating 
these signals could be the same as are al- 
ready in wide use in serial computers. 
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trical Engineers of Japan. 
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Amos Nathan (S’50-A’51-M’55) was 
born in Liegnitz, Germany, on October 5, 
1920. He entered Palestine in 1933 and re- 
ceived the degree of Dipl. El. Mech. Engi- 
neer from Technion, Israel Institute of Tech- 
nology, Haifa, Israel, in 1943. 

He served with the British Forces from 
1942 to 1946. From 1946 to 1948, he was an 
assistant at Technion. He then served in the 
Israel Defence Force until 1953, attaining 
the rank of major. After a year of study in 
the United States, he received the M.S. de- 
gree in electrical engineering from Columbia 
University in 1951. In 1956, he received the 
degree of doctor of science from Technion. 

Dr. Nathan rejoined the staff of Tech- 
nion in 1953, where he is now a senior lec- 
turer in electrical engineering, and he also 
has charge of an electronic differential ana- 
lyzer project. 

Dr. Nathan is an associate member of 
IEE, London. 
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Takeo Miura was born in Kyoto, Japan, 
on October 1, 1925. He graduated from 
Kyoto University, Japan, with the BS. 
degree (Kogakushi). He joined Hitachi, 
Ltd., in 1949 and was engaged in the design 
of eaecon motors at Kameido Works. 
In 1950, he was transferred to Central 
Reseich Laboratory, where he has since 
been engaged in the development of elec- 
tronic analog computers, regulators, and 
servomechanisms. 

Mr. Miura is a member of the Institute 
of Electrical Engineers of Japan. 
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PGEC News 


1957 WESCON 


The 1957 Western Electronic Show and 
Convention will be held in San Francisco 
for the four days of August 20-23. The dead- 
line for papers has passed. Inquiries about 
the meeting may be directed to the business 
office, 342 North La Brea, Los Angeles 36, 
Calif. 


DISTRIBUTION OF PGEC MEMBERS; 
Patp MEMBERSHIP BY SECTIONS 
AS OF DECEMBER 31, 1956 


: No. o 
Sections Mi pa 

Akron 29 
Alamogordo-Holloman ila 
Albuquerque-Los Alamos 26 
Atlanta 39 
Baltimore 114 
Bay of Quinte 4 
Beaumont-Port Arthur W 
Binghamton 61 
Boston 448 
Buenos Aires 2 
Buffalo-Niagara 26 
Cedar Rapids 19 
Central Florida 42 
Chicago 160 
China Lake 10 
Cincinnati 33 
Cleveland 35 
Columbus 13 
Connecticut Valley 118 
Dallas 31 
Dayton 40 
Denver 19 
Detroit 120 
Egypt 0 
Elmira-Corning 4 
El Paso 18 
Emporium 19 
Evansville-Owensboro 6 
Florida West Coast 1 
Ft. Huachuca i 
Ft. Wayne 16 
Ft. Worth Bil 
Hamilton 4 
Hawaii 1 
Houston 44 
Huntsville 30 
Indianapolis 22 
Israel 10 
Ithaca 22 
Kansas City 24 
Little Rock 9 
London 4 
Long Island 269 
Los Angeles TNS 
Louisville 5 
Lubbock 6 
Miami 12 
Milwaukee 49 
Montreal 38 
Newfoundland M6 
New Orleans 13 
New York 618 


eee No. 0 
Sections ete as : 

Northern Alberta 3 
Northern New Jersey 252 
Northwest Florida 7 
North Carolina-Virginia 33 
Oklahoma City 9 
Omaha-Lincoln 5 
Ottawa 14 
Philadelphia 447 
Phoenix 10 
Pittsburgh 52 
Portland 14 
Princeton 29 
Regina 1 
Rio de Janeiro i 
Rochester DS 
Rome-Utica 37 
Sacramento 7 
St. Louis 25 
Salt Lake City 5 
San Antonio 13 
San Diego 47 
San Francisco 284 
Schenectady 36 
Seattle 42 
South Bend-Mishawaka 10 
So. Alberta 2 
Syracuse 63 
Tokyo 22 
Toledo 8 
Toronto 39 
Tucson 10 
Tulsa Di 
Twin Cities 103 
Vancouver 6 
Washington 231 
Wichita 6 
Williamsport 1 
Winnipeg 5 
Total Paid Members 5625 


RESUME OF PGEC ADMINISTRATIVE 
COMMITTEE MEETING, FEBRUARY 
21 Sy. 


The PGEC Administrative Committee 
met during the Western Joint Computer 
Conference in Los Angeles on February 27, 
1957. In the absence of Chairman J. D. Noe, 
who was ill, Vice-Chairman W. Buchholz 
presided. Twelve voting members were 
present or represented by proxy. Several 
chapter and committee chairmen were wel- 
comed as guests. 

Secretary-Treasurer R. Y. Wing reported 
a financial balance of $24,920 as of Decem- 
ber 31, 1956. Of this, $6,000 has been ear- 
marked for the IRE Computer Fellowship 
for the next two years. Nearly $9,000 was 
spent on publications during the year. Not 
counting committed funds, the net balance 
increased by almost $1,000 during 1956. 

The PGEC membership stood at 5494 at 
the end of 1956, including 498 student mem- 
bers. 


Nominations Committee Chairman J. H. 
Felker presented a slate of new Administra- 
tive Committee members and officers for the 
coming year. There were no other nomina- 
tions. Voting takes place during the March 
meeting. 

The Committee approved a change in 
the PGEC by-laws to permit the PGEC 
Chairman to be elected during his third 
term on the Committee, so that he serves as 
Chairman during a fourth year and as Past- 
Chairman during a fifth year on the Com- 
mittee. This follows a recent change in the 
PGEC Constitution permitting the Admin- 
istrative Committee to consist of a maxi- 
mum of 18 members instead of 15, as before. 

Robert Roggenbuck from the Detroit 
Chapter was appointed to fill out the term 
of Norman Scott, who resigned. 

The new IRE Affiliate plan was adopted 
for the PGEC. An initial list of professional 
societies whose members will be eligible to 
join the PGEC as nonvoting Affiliates and 
receive the PGEC TRANSACTIONS (but not 
the IRE PROCEEDINGS) for a yearly fee of 
$6.50 includes: 


Association for Computing Machinery 

Society for Industrial and Applied 
Mathematics 

American Mathematical Society 

Operations Research Society of America 

National Association of Cost Account- 
ants 

National Machine Accountants Associa- 
tion 

American Management Association 

Institute of Aeronautical Sciences 

American Physical Society 

American Society of Mechanical Engi- 
neers 

Society of Automotive Engineers. 


The PGEC, for the time being, will not 
include advertising in the TRANSACTIONS 
since the additional revenue is not needed. 

WERNER BUCHHOLZ 


PGEC ADMINISTRATIVE 
COMMITTEE MEETING 


The following is a summary of the meet- 
ing of the PGEC Administrative Committee 
held in New York, N. Y., on March 18, 1957. 
A quorum of 13 members and proxies out of 
the 16 total members was present. 


Meeting Notices 


A letter from George Bailey, Executive 
Secretary, IRE, was read to the Administra- 
tive Committee, in which Mr. Bailey out- 
lined the change in IRE requirements for 
publicizing Chapter meetings. Heretofore, in | 
some Sections it has been considered essen- | 
tial to notify all members of the Section for 
every Chapter meeting, and the expense in- 
volved has prohibited the Section from pay- 
ing for occasional meeting notices sent to 
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members of a particular Chapter. It is the 
PGEC Administrative Committee’s present 
understanding that it is no longer an abso- 
lute requirement for the Section to notify all 
members for each meeting, therefore, it 
should be possible for the Section to handle 
the cost of an occasional special mailing to 
publicize a change in plan for a Chapter. 


IRE Predoctoral Computer Fellowship 


Dr. Willis H. Ware, Chairman of the 
Student Activities Committee, reported that 
the PGEC Fellowship is now officially estab- 
lished and the National Academy of Sciences 
is handling the fellowship for the PGEC. 
The award for the 1957-1958 fellowship is 
expected in April. 


Funds for PGEC Chapters 
The Committee approved a method of 
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providing financial assistance to the PGEC 
Chapters, since it is now possible to send 
funds to a Section earmarked for a particu- 
lar Chapter within the Section. The Nation- 
al PGEC will match IRE Headquarters’ 
funds ($10 per meeting attended by more 
than 25 members, up to five meetings per 
year) for each PGEC Chapter meeting. Pay- 
ment will be made upon submitting meeting 
reports, including date and number of people 
attending, to Stanley B. Disson, Chairman 
of the PGEC Sectional Activities Commit- 
tee, Burroughs Corp., Great Valley Labora- 
tories, Paoli, Pa. 


Elections 


New members and officers of the PGEC 
Administrative Committee were elected ef- 
fective April 1, 1957. The results are tabu- 
lated inside the front cover of this issue. 
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Publications 


As one of several moves to encourage the 
flow of worthwhile papers to the PGEC 
Transactions, the Committee reaffirmed the 
existing publication policy. Publication in a 
convention record or proceedings of a meet- 
ing shall not prejudice republication in the 
TRANSACTIONS, especially where the conven- 
tion records or proceedings of the meeting 
are sent only to a small percentage of PGEC 
membership. 


Next Meeting 


The next meeting of the PGEC Adminis- 
trative Committee will be held during 
WESCON in San Francisco during the pe- 
riod August 20-23, 1957. 

JERRE D. Nor 
Past Chairman, IRE-PGEC 
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Reviews of Current Literature 


Copies of books or of articles to be reviewed should be sent to H. D. Huskey, Depart- 
ment of Electrical Engineering, University of California, Berkeley, California. The 
editors wish to express their gratitude to the reviewers who, through their efforts, make 
this section possible-—Harry D. Huskey 


GENERAL 


57-31 
Moderne Rechenmaschinen—Hans 
Biickner. (Handbuch der Physik. Bd. II, 
Springer-Verlag, Berlin-Géttingen-Heidel- 
berg, pp. 471-498; 1955.) This review article 
is divided into two sections, covering analog 
and digital computers respectively. The em- 
phasis is, rightly, on electronic internally- 
programmed digital calculators. The history 
of the subject is briefly outlined, from Bab- 
bage and Lady Lovelace till about 1954. A 
discussion of programming with examples of 
coding is based on the system used for the 
EDSAC in Cambridge, England. The logical 
structure of arithmetic operations in com- 
puters, use of binary numbers, representa- 
tion of truth functions are among the theo- 
retical topics. A section is devoted to the 
various kinds of memory employed. 
W. Freiberger 
Courtesy of Mathematical Reviews 


57-32 
Gedanken-Experiments on Sequential 
Machines—Edward F. Moore. (Automata 
Studies, Annals of Math. Studies, no. 34, pp. 
129-153; 1956.) This paper studies distin- 
zuishability of sequential machines and dis- 
tinguishability of internal states within a 
machine by experiments performed on the 
machines. It is shown that, for general ma- 
chines, no finite experiment will suffice to 
identify one machine among the class of se- 
quential machines (Theorem 2). By restrict- 
ing the class of machines to those which are 
“strongly connected” and in “reduced form,” 
it is then shown that there is some finite ex- 
periment with copies of a machine which will 
suffice to distinguish the machine from oth- 
srs in the class (Theorem 9). Several other 
theorems on distinguishability and on the 
length of experiments are proved. 
Cuyenilee 
Courtesy of Mathematical Reviews 


57-33 

Some Uneconomical Robots—James T. 
Culbertson. (Automata Studies, Annals of 
Math. Studies, no. 34, pp. 99-116; 1956.) 
Using a number of each of three types of 
1eurons, receptor, central, and effector, and 
with proper interconnection of these neurons, 
the paper shows how a memoryless robot can 
9e constructed. Then, by chopping time into 
discrete levels and increasing the number of 
‘eceptor and central neurons, it is shown how 
1 robot with complete memory (a record, at 
iny time ¢, of all past receptor firings) can be 


constructed. A way of introducing proba- 
bilistic behavior without inserting unreliable 
elements is also indicated. 

C. Y.'Lee 


Courtesy of Mathematical Reviews 


681.142 57-34 
Design of Computer Circuits for Re- 
liability—W. Renwick. (Electronic Eng., vol. 
28, pp. 380-384; September, 1956.) Factors 
influencing the choice of components, and 
precautions taken during the initial circuit 
design and mechanical construction of the 
EDSAC II machine are described. 
Courtesy of Proc. IRE 
and Wireless Engineer 


57-35 
The Application of Matrix Methods to 
Multi-Variable Control Systems—R. J. 
Kavanagh. (J. Franklin Inst., vol. 262, pp. 
349-367; November, 1956.) This paper pre- 
sents the methods of matrix algebra which 
may be used for the description, analysis, 
and synthesis of the general ease of linear 
multivariable control systems. Some prac- 
tical examples of the analysis and synthesis 
of systems are given and criteria for the 
realizability of multivariable control sys- 
tems with desired characteristics are derived. 
The procedure adopted in this method is to 
describe the multivariable system by means 
of a matrix transfer function linking input 
to output. A computer matrix is then in- 
serted between true inputs and the actual in- 
puts to the multivariable system. The com- 
puter matrix forms linear combinations of 
the true inputs and the multivariable system 
outputs to form the desired actual inputs to 
the multivariable system. The synthesis 
problem reduces to the determination of the 
computer matrix elements. This paper is 
clearly written and provides information of 
considerable value to all who are engaged 
in the design of control systems, both analog 
and digital. 
Gene M. Amdahl 


57-36 
The Mathematics of Information Theory 
—Brockway McMillan. (1955 IRE Na- 


TIONAL CONVENTION RECORD, pt. 4, pp. 48- 
51.) It seems superfluous to review an article 
which is itself nothing but a review of the 
broad aspects of a complex subject. Need- 
less to remark, there is no mathematics in 
the article. There are no new approaches 
to any of the basic problems of information 
theory suggested; in fact, the limitations of 
present attempts at analysis are rather 


fatalistically accepted. The main point of 
the article is to observe that engineers and 
physicists have been working on the same 
problems as some mathematical statisticians, 
but with constraints which for some reason 
or other mathematical statisticians do not 
seem to permit themselves. Wald’s work is 
lauded as a unifying force, but is observed 
to be responsible for more philosophizing 
than discovery. In short, this article will 
pass the censor’s office, but will not be a 
big hit at the box office. 

R. H. DeLano 


57-37 
The Vernier Time-Measuring Technique 
—R. G. Baron. (Proc. IRE, vol. 45, pp. 21- 
30, January, 1957.) This paper describes a 
technique for precision time measurement of 
the interval between two pulses. The system, 
which is accurate to 10 millimicroseconds, 
employs an electrical equivalent of the me- 
chanical vernier scale to permit this ac- 
curacy while counting 2-mc pulses. This is 
accomplished by utilizing a coarse counter 
to count the 2-mc pulses between the start 
and stop pulses and a vernier counter to 
count the 2-mc pulses between the start or 
stop pulse and coincidence of a 2-mc pulse 
and a pulse from the vernier pulse generator. 
The vernier pulse generator is started by the 
start or stop pulse and produces pulses at a 
slightly greater rate than the 2-mc pulse 
generator until coincidence. The use of a 
single vernier frequency is made possible by 
proper control of the counter. In the interest 
of giving credit where credit is due, it should 
be noted that a paper, “Precision Automatic 
Time Measurement Equipment,” was pre- 
sented at the IRE National Convention, 
March 5, 1952, by D. W. Burbeck and W. E. 
Frady describing an almost identical system 
which was fabricated at Hughes. 
A. D. Scarbrough 


57-38 

High-Speed Digital Conversion—M. L. 
Klein, F. K. Williams, and H. C. Morgan. 
(Instruments and Automation, vol. 29, pp. 
1297-1302; July, 1956.) This is the third 
article of the series on analog-digital con- 
version by the authors. (See 56-158 and 
56-159, vol. EC-5, p. 261; December, 1956.) 
It is a survey of available equipment for 
high-speed data logging and digital conver- 
sion. The paper is divided into two sections, 
all-electronic data loggers and all-electronic 
high-speed converters. Four commercially 
available units are described in each cate- 
gory. Voltage conversion rates range from 
20-100,000 per second with an accuracy of 
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0.1 per cent. Block diagrams are given with 
a short outline of system operation. More 
detailed information in one or two cases 
would be helpful to the reader. However, the 
authors on the whole have presented a very 
excellent survey of existing equipment. 
Raymond Davis 


621.314.63 +621.314.7 57-39 
Maximum Power of Semiconductor 
Junction Elements—J. P. Vasseur. (Ann. 
Radio-élect., vol. 11, pp. 3-28; January, 
1956.) Universal curves are presented from 
which the maximum permissible power dis- 
sipation can be determined for Ge junctions 
cooled by 1) conduction or 2) convection. 
The influence of mounting, circuit arrange- 
ment, and signal waveform is discussed. The 
time taken to reach temperature equilib- 
rium is investigated. Both diodes and 
transistors are considered. 
Courtesy of Proc. IRE 
and Wireless Engineer 


621.314.7 57-40 
The Effect of Surface Treatments on 
Point-Contact Transistor Characteristics— 
J. H. Forster and L. E. Miller. (Bell Sys. 
Tech. J., vol. 35, pp. 767-811; July, 1956.) 
The properties of formed point contacts on 
Ge are discussed; technique for observing 
the equipotentials surrounding such con- 
tacts is described. Differences between 
donor-free and donor-doped contacts are 
emphasized. Differences are also observed 
between unformed point contacts subjected 
to different surface treatments. These con- 
siderations are applied to a study of transis- 
tor point-contact forming in contemporary 
practice. High yields from the forming 
process can be expected with oxidized sur- 
faces; the use of chemical washes to remove 
soluble Ge-oxide surface films greatly re- 
duces the forming yields. 
Courtesy of Proc. IRE 
and Wireless Engineer 


621.314.63:546:28 57-41 
Silicon Wide p-n Junction—J. Nishizawa. 
(Sci. Rep. Res. Inst. Tohoku Univ., Ser. B, 
vol. 6, pp. 183-215; 1955.) The effect of the 
depletion layer in Si p-” junctions is investi- 
gated by determining the V/J character- 
istics for various positions of the contact 
whisker, using a specially prepared junction 
having a very wide high-resistance layer. 
Courtesy of Proc. IRE 
and Wireless Engineer 


ANALOG COMPONENT 
RESEARCH 


681.142 57-42 
A Circuit for Analogue Formation of xy/Z 
—M. J. Somerville. (Electronic Eng., vol. 28, 
pp. 388-389; September, 1956.) “A quarter 
squares multiplier, using a triangle carrier 
waveform in the squaring circuits is ex- 
tended to give division simultaneously with 
multiplication. This is achieved by control- 
ling the slope of the triangle carrier waveform 
so as to be proportional to the divisor Z.” 
Courtesy of Proc. IRE 
and Wireless Engineer 
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ANALOG EQUIPMENT 


57-43 
Aircraft Performance Studied on an 
Electronic Analog Computer—L. B. Wadel 
and C. C. Wan. (Proc. Western Joint Com- 
puter Conf., Los Angeles, Calif., pp. 78-82; 
1955.) The relations between aircraft per- 
formance (maximum speed, rate of climb, 
fuel consumption) aerodynamic data, and 
engine data are discussed. The analog com- 
puter components for the solution of these 
equations are outlined. The equations are 
mostly nonlinear algebraic, so function 
generators and multipliers are the most fre- 
quently used components. The complete 
computer is a cross between a special-pur- 
pose and a general-purpose analog machine. 
Some iterative techniques for computing 
range are mentioned. The computations 
have two uses: for the preliminary design of 
aircraft, and for presenting the performance 
data of the finished aircraft to the pilot ina 
handbook form. No consideration is made 
of desirable input-output equipment for 
the computer. This equipment is probably 
one of the most important parts for practical 
use of the machine. In addition no indica- 
tion is given why an analog computer is 
superior to a general-purpose digital com- 
puter. 
Max V. Mathews 


57-44 

On the Continuous Solution of Integral 
Equations by an Electronic Analog. I— 
Michael E. Fisher. (Proc. Camb. Phil. Soc., 
vol. 53, pp. 162-174; January, 1957.) This 
paper describes the practical realization of a 
machine for solving a class of integral equa- 
tions by analog techniques. It is based on 
an idea due to Wallman and uses an inter- 
esting device developed by Bergman for 
the storage and reproduction of analog wave- 
forms. The machine can solve Fredholm 
integral equations of the second kind by a 
new form of iterative method which is 
claimed to be superior to that due to Neu- 
mann and the time taken for convergence 
is of the order of one second or less. Con- 
vergence and error analysis are discussed. 


Stanley Gill 


UTILIZATION OF ANALOG 
EQUIPMENT 


57-45 
The Representation of Constraints by 
Means of an Electronic Differential Ana- 
lyzer—D. T. Greenwood. (IRE TRANs., 
vol. EC-5, pp. 111-113; September, 1956.) 
The use of high-gain amplifiers is shown to 
be helpful in the representation of con- 
straints. This method enables one to repre- 
sent constrained systems in a manner such 
that all coordinates are available for the 
application of arbitrary forces or displace- 
ments. The procedure is explained by means 
of an example. 
Courtesy of Proc. IRE 


DIGITAL COMPONENT 
RESEARCH 
57-46 
The Physics of Solid-State Light 
Amplifiers—D. A. Cusano. (IRE TRAnNs., 
vol. NS-3, pp. 102-106; November, 1956.) 
An essential component of present light am- 
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plifying screens is a layer of inorganic, lumi- 
nescent material. The ability of this layer 
to derive energy from an electric field to 
which it is subjected and convert this energy 
into visible light (z.e., electroluminescence) 
is fundamental to, but not alone sufficient 
for the realization of solid-state light ampli- 
fication. What must be included is some 
means whereby the luminescent response to 
the field can be strongly controlled by inci- 
dent radiation. In one case this has been ob- 
tained through the discovery of a phosphor 
film with the necessary properties, and in 
another by the utilization of a photoconduct- 
ing material in contact with an electro- 
luminescent layer. In the latter case the 
photoconductor is the radiation-sensitive ele- 
ment, which serves to modify the potential 
across the luminescent component. 

The main characteristics of these light 
amplifying screens will be described, with 
particular emphasis on the interpretation 
of the phenomenon observed with the single 
phosphor film. The behavior, as it is rele- 
vant to their possible use in radiation detec- 
tion, will be discussed. 

Courtesy of Proc. IRE 


57-47 

Quarterly Report No. 12, Second Series 
—J. R. Bowman, C. H. T. Wilkins, e¢ al. 
(Quart. Rep. Computer Components Fellow- 
ship Mellon Inst., 45 pp.; July 1, 1956 to 
September 30, 1956.) In the first of these, 
the results of time and temperature studies 
on the preparation of printed enamel capaci- 
tors are discussed. Initial results obtained 
using a “flash evaporation” technique for 
the preparation of alloy resistor films by 
vacuum evaporation and deposition are re- 
ported in the second section. The third 
section contains descriptions of methods 
used for the preparation of stannic oxide re- 
sistor films. In Part II, changes in brightness 
wave characteristics that occur between 30 
and 120 volts at 60 cps in electroluminescent 
cells having eight different metal backings 
are discussed. Over this voltage range, the 
disturbance ratio and the phase angle are 
found to be independent of the choice of 
electrode metal. Also, as the voltage is in- 
creased 1) the disturbance ratio decreased, 
2) the modulation ratio decreased to zero 
then increased, 3) there is an over-all increase 
in the phase angle, but only in the nickel 
backed cell is this increase a smooth one. 
Cy Hats Wilkins 


57-48 

A Study of the Neon Bulb as a Non- 
linear Circuit Element—C. FE. Hendrix. 
(IRE TRANS., vol. CP-3, pp. 44-54; Septem- 
ber, 1956.) The ordinary NE-2 neon bulb 
possesses two-valued properties which make 
it valuable for use in logical gate circuits 
such as are used in digital computers. The 
cost of the NE-2 is of the order of $0.05 
compared to an average unit cost of $1.00 
for germanium diodes. In the logical gate 
application, circuit performance is relatively 
insensitive to parameter variations among 
bulbs. Both “or” and “and” gates, and com- 
binations thereof can be built, and operation 
up to 30 ke easily obtained. Statistical stud- 
ies of bulbs point out that the aging effects 
to be expected are a rise in voltage drop 
across the bulb and a reduction in the varia- 
bility among bulbs. The change appears 
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rather suddenly, after about 1000 milli- 
ampere-hours of operation. Equivalent cir- 
cuits are presented, and simple circuit de- 
sign techniques have been worked out. 
Courtesy of Proc. IRE 


57-49 
Transmission Secondary Electron Multi- 
plication for High-Speed Pulse Counting— 
E. J. Sternglass and M. M. Wachtel. (IRE 
Trans., vol. NS-3, pp. 29-32; November, 
1956.) A new type of high-speed electron 
multiplier employing transmission second- 
ary electron emission from thin insulating 
films is described. Electrons from a photo- 
surface are multiplied in a series of plane- 
parallel dynodes consisting of a thin scatter- 
ing layer of a heavy metal followed by a 
layer of a pure crystalline insulating ma- 
terial of high secondary emission yield. The 
characteristics of such dynodes are described 
and the de and pulse performance of an 
experimental seven-stage device employing 
this principle is presented. 
Courtesy of Proc. IRE 


57-50 
Supermendur, A New Rectangular- 
Loop Magnetic Material—H. L. B. Gould, 
D. H. Wenny. (Electrical Eng., vol. 76, pp. 
208-211; March, 1957.) This article dis- 
cusses the characteristics of a vanadium 
(2 per cent)—iron (49 per cent)—cobalt (49 
per cent) alloy named Supermendur. The 
electrical characteristics include: low coer- 
cive force, 0.26 oersted; low hysteresis loss; 
high maximum permeability, 66,000 at 
20,000 Gausses; high remanence, 21,150; 
saturation 24,000 Gausses; flux swing 45,150. 
The article describes the melting, processing, 
heat treating, electrical tests, and structure 
of this material. The material has been 
rolled to thicknesses of 0.3 mil. 
Tei warson 


57-51 
Modern Magnetic Ferrites and Their 
Engineering Applications—C. D. Owens. 
(IRE TrRans., vol. CP-3, pp. 54-62; Sep- 
tember, 1956.) Ferrites and their applica- 
tions constitute a rapidly developing field. 
This paper reviews the present status of the 
development, properties, and applications 
of ferrites and their significance to modern 
engineering. Magnetic ferrites were de- 
veloped from a “new look” at lodestone, the 
oldest magnetic material known to man. 
Proper control of composition and heat 
treatment yields ceramic type bodies with 
high permeabilities, which were formerly 
obtained only in metals. This material has 
brought the physicist new insight into the 
mechanisms of magnetism, and has given 
the design engineer a variety of new ma- 
terials with unprecedented combinations 
of properties. The unique characteristics 
of ferrites are bringing about revolutionary 
design techniques in inductors and trans- 
formers for use at carrier, radio, and video 
frequencies, and have made possible an 
important new family of nonreciprocal cir- 
cuit elements in microwave applications. 
Courtesy of Proc. IRE 
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57-52 
High-Speed Shift Registers Using One 
Core Per Bit—V. L. Newhouse and N. S. 
Prywes. (IRE TrANs., vol. EC-5, pp. 114- 
120; September, 1956.) A three, and a two 
winding per core, high-speed, current driven, 
one core per bit shift register is presented to- 
gether with an analysis of the basic circuit 
involved. An intermediate storage capacitor 
is used between successive logical elements. 
The charge and discharge of this capacitor 
are controlled in a positive manner by volt- 
age blocking pulses. The undesired feedback 
of energy from one stage to earlier stages is 
thereby prevented, giving high efficiency of 
operation. The three winding per core regis- 
ter described is reversible and capable of 
operating in the megacycle range. The appli- 
cation of the basic shift register element to 
computer logic applications is described. 
Courtesy of Proc. IRE 


621.314.63 57-53 
High-Frequency Silicon-Aluminum Al- 
loy Junction Diodes—M. B. Prince. (IRE 
TRANS., vol. ED-2, pp. 8-9; October, 1955.) 
A preliminary account of experiments on 
diodes made by alloying the point of a thin 
Al wire to a wafer of n-type Si; the diameter 
of the alloyed region is probably <0.0005 
inch. Rectification at frequencies up to 500 
mc has been obtained. The reverse satura- 
tion current is about 10~° A. 
Courtesy of Proc. IRE 
and Wireless Engineer 


621.314.63:546.28 57-54 
Reverse Current and Carrier Lifetime as 
a Fuction of Temperature in Silicon Junc- 
tion Diodes—E. M. Pell and G. M. Roe. 
(J. Appl. Phys., vol. 27, pp. 768-772; July, 
1956.) Previous work on Ge diodes (Ab- 
stract 3424 of Proc. IRE, 1955) is extended 
to Si grown-junction diodes; measurements 
over the temperature range —190° to 
+200°C are reported. The lifetime becomes 
constant over a low-temperature range; this 
is consistent with published recombination 
theory. The results indicate that charge 
generation from centers about 0.5 ev deep 
is responsible for most of the reverse current 
at temperatures up to well above room tem- 
perature. 
Courtesy of Proc. IRE 
and Wireless Engineer 


621.314.63 57-55 
On the Transient Behaviour of Semicon- 
ductor Rectifiers—B. R. Gossick. (J. Appi. 
Phys., vol. 27, pp. 905-911; August, 1956.) 
The investigation reported previously (Ab- 
stract 1236 of Proc. IRE, 1955) is extended 
to cover large-amplitude transients in sur- 
face-barrier and point-contact diodes. 
Courtesy of Proc. IRE 
and Wireless Engineer 


621.314.7:621.3.015.3 57-56 

Theory of Transient Process in Semi- 
conductor Triodes—E. [. Adirovich and 
V..G. Kolotilova. (G.-R. Acad. Sct. U.R.S-S., 
vol. 108, pp. 629-632; June, 1956. In Rus- 
sian.) Expressions are derived for the trans- 
fer characteristics of a transistor operating 
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in grounded-emitter and grounded-collector 
circuits. The corresponding expression for 
the grounded-base circuit was given pre- 
viously. (Abstract 1892 of Proc. IRE, 1955.) 
Courtesy of Proc. IRE 

and Wireless Engineer 


57-57 
Logic Circuits for a Transistor Digital 
Computer—G. W. Booth and T. P. Both- 
well. (RE Trans., vol. EC-5, pp. 132-138; 
September, 1956.) The reliability and per- 
formance which can be achieved in high- 
speed switching circuits using presently 
available high-frequency junction transis- 
tors suggests the use of transistor circuits 
in preference to more conventional circuits 
for many applications. The area of airborne 
digital computers is one in which the physi- 
cal characteristics of the transistor as well 
as of associated low dissipation components 
can be most favorably exploited. We have 
presented here a group of circuits which 
fulfill the requirements of speed, low dissipa- 
tion, size, and weight for most such applica- 
tions. The circuits shown will operate over 
the range —30°C to +60°C, and actually 
have been operated over the range —50°C 
to —90°C. Low dissipation of the circuits 
places minimum requirements on power 
supplies, and temperature control may be 
achieved with only a small amount of cool- 
ing. An estimate can be obtained from an 
example. For a large-scale computer, say 
2000 transistors, 800 of which are in bistable 
circuits, total dissipation will be less than 
140 watts. 
Courtesy of Proc. IRE 


57-58 
Semi-Conductor Diode Amplifier Con- 
siderations—Henry W. Kaufman. (1955 
IRE NationaL CONVENTION REcoRD, pt. 
4, pp. 146-149; 1955.) The author presents 
some experimental data on forward and 
reverse transients which was originally taken 
to aid in the design of conventional diode 
gating circuits. This data, which includes 
wafer thickness as a parameter, is then inter- 
preted from the novel standpoint of the 
design of diode amplifiers. In one example, 
the author calculates that a circuit using 
a 1N191 input diode, a 1N92B amplifier 
diode, and a resistor will give an energy 
gain of about 18. This is for a one-micro- 
second input and a two-microsecond output 
pulse. 
A. W. Holt 


57-59 

High-Speed Flip-Flops for the Mi£illi- 
microsecond Region—Z. Bay and N. T. 
Grisamore. (IRE TRANs., vol. EC-5, pp. 
121-125; September, 1956.) The problem 
of designing high-speed flip-flops has been 
approached by dividing the circuit operation 
into steady-state and switching functions, 
the steady-state function being assigned to 
a slave flip-flop and the switching function 
to a driving circuit. Circuits, using conven- 
tional components and having a resolving 
time of 10 mysec are described. Resolving 
times as low as 2 musec have been attained 
by using a special beam-deflection tube as 
the slave flip-flop. These circuits dissipate 
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considerably less power in the intervals be- 
tween switching than conventional circuits. 
Courtesy of Proc. IRE 


57-60 

VHF Pulse Techniques and Logical Cir- 
cuitry—D. E. Rosenheim and A. G. Ander- 
son. (Proc. IRE, vol. 45, pp. 212-219; 
February, 1957.) This interesting paper 
presents techniques and components for 
handling 10 millimicrosecond pulses at 50 
mc repetition rates. The active circuits make 
use of the Philips EFP-60 secondary emis- 
sion pentode. The EFP-60 has an additional 
active electrode, the dynode which allows 
the tube to be used as a noninverting ampli- 
fier with both current and voltage gain and 
a grid to dynode transconductance of 20,000 
pmhos The Western Electric 417A/5842 
triode with g,,=25,000 wmhos is used as a 
cathode follower for the fast pulses. Logical 
gating is performed with selected gold- 
bonded germanium diodes having high for- 
ward conductance and good transient re- 
sponse. Storage is accomplished by electro- 


magnetic delay lines and EFP-60 reshaping 
circuits. A number of schematics of other 


computer circuits are given together with 
helpful hints on generating, observing, and 
transmitting these extremely fast pulses. 

A. D. Scarbrough 


537.226/.227:546.431.824-31 57-61 

Polarization Reversal in the Barium 
Titanate Hysteresis Loop—R. Landauer, 
D. R. Young, and M. E. Drougard. (J. 
Appl. Phys., vol. 27, pp. 752-758; July, 
1956.) The af hysteresis loop of tetragonal 
BaTiO; is examined in relation to Merz’s 
finding (Abstract 445 of Proc. IRE, 1955) 
that the rate at which the polarization re- 


verses in an applied field of intensity E is 


proportional to e“”, where a depends on 
temperature. The switching rate adjusts to 
variations in the applied field with a time 
lag whose order of magnitude does not ex- 
ceed 107-8 sec. The results imply that do- 
mains cross the crystal with a speed of the 
order of that of sound. Increases in dielec- 
tric constant measured during switching re- 
main unexplained. 
Courtesy of Proc. IRE 
and Wireless Engineer 


537.226/.228.1 57-62 
Effect of Hydrostatic Pressure on the 
Hysteresis Loop of Guanidine Aluminium 
Sulphate Hexahydrate—W. J. Merz. (Phys. 
Rev., vol. 103, pp. 565-566; August, 1956.) 
Measurements at a frequency of 60 cps indi- 
cate that both the spontaneous polarization 
and the coercive field strength of this ma- 
terial increase with pressure up to about 500 
atm. 
Courtesy of Proc. IRE 
and Wireless Engineer 


537.226/.227 57-63 
Ferroelectricity in Ammonium Sulphate 
—B. T. Matthias and J. P. Remeika. (Phys. 
Rev., vol. 103, p. 262; July 1, 1956.) A brief 
note reporting that (NH4)2S5O4 becomes 
ferroelectric below its transition point at 
—49.5°C. 
Courtesy of Proc. IRE 
and Wireless Engineer 
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537.226/227:546.431.824-31:539.234 57-64 
Time Changes in Thin Films of BaTiO;— 
C. Feldman. (J. Appl. Phys., vol. 27, pp. 
870-873; August, 1956.) “The decrease of 
dielectric constant with time under an 
applied alternating field, less than the coer- 
cive field, has been studied in films of BaTiO; 
between 1 and 3 yp thick. The phenomenon 
may be interpreted as being associated with 
the process of switching the domains to a 
position more nearly parallel to the applied 
field.” 
Courtesy of Proc. IRE 
and Wireless Engineer 


537.226/.227:546.431.824-31 57-65 
Switching Time in Ferroelectric BaTiO; 
and its Dependence on Crystal Thickness— 
W. J. Merz. (J. Appl Phys., vol. 27, pp. 
938-943; August, 1956.)“The switching 
time ¢, and the switching current tmax have 
been measured as a function of applied 
field E and of the size of the sample. It has 
been observed that the “activation field” a 
for the nucleation of new domains is inversely 
proportional to the thickness of the sample. 
This behavior can be explained by assuming 
a surface layer. The thickness of this layer 
has been calculated to be of the order of 1074 
cm. The same way we can explain the thick- 
ness dependence of the 60-cps coercive field 
strength. Furthermore, it has been found 
that the switching time depends to a first 
approximation linearly on the thickness of 
the sample if the field E is kept constant. 
This can be explained by assuming a domain 
wall motion primarily in the forward direc- 
tion or by assuming a nucleation mechanism. 
The maximum velocity of the domain growth 
was found to be of the order of the velocity 
of sound. The switching time does not de- 
pend on electrode area.” 
Courtesy of Proc. IRE 
and Wireless Engineer 


681.142 +-621.397.5]:621.395.625 57-66 
Recording and Reproduction of Fre- 
quencies above 100 kc/s on Magnetic 
Media. Application to Storage Elements and 
to the Recording of Television Images— 
R. Charlet. (Rev. Gén. Elect., vol. 65, pp. 
359-365; June, 1956.) A combined summary 
of papers presented at a conference on the 
recording of sound and information held in 
Paris in April, 1955; current techniques are 
described and necessary improvements are 
indicated. 
Courtesy of Proc. IRE 
and Wireless Engineer 


57-67 
500,000,000—Bit Random—Access Mem- 
ory—George E. Comstock. (Instruments 
and Automation, vol. 29, pp. 2208-2211, 
November, 1956.) This paper describes a 
large random access memory. The storage 
medium is standard 3-inch instrument— 
quality Mylar-base magnetic recording 
tape. Strips of tape approximately two feet 
long are strung vertically on stainless steel 
frames. About 40,000 feet of tape are acces- 
sible on a random basis within less than a 
second, with average access time less than 
% second. 
Raymond Davis 
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DIGITAL SYSTEMS 


RESEARCH 
57-68 
A Topological Method for the Determi- 
nation of the Minimal Forms of a Boolean 
Function—R. H. Urbano and R. K. Mueller. 
(IRE Trans., vol. EC-5, pp. 126-132; Sep- 
tember, 1956.) The topology of the -dimen- 
sional cube is used to reduce the problem of 
determining the minimal forms of a Boolean 
function of » variables to that of finding the 
minimal coverings of the essential vertices 
of the basic cell system associated with the 
given function. The proof of this statement 
is contained in the central Theorem 4. A nu- 
merical easily programmed procedure is 
given with which it is possible to treat prob- 
lems with a greater number of variables 
than has heretofore been practical. The 
procedure bypasses the determination of the 
basic cells (the prime implicants of W. V. 
Quine) and locates the essential vertices, 
from which in turn the irredundant and 

minimal forms are obtained. 
Courtesy of Proc. IRE 


57-69 

Circuit Minimization: Minimal and 
Irredundant Boolean Sums by Alternative 
Set Method—E. W. Samson and R. Mueller. 
(Communications Laboratory, Electronics Re- 
search Directorate, AF Cambridge Res. Cen- 
ter, Cambridge, Mass., Tech. Rep. 55-109, 
v+11 pp.; 1955.) 

Circuit Minimization: Sum to One 
Process for Irredundant Sums—E. W. 
Samson and R. K. Mueller. (Communica- 
tions Laboratory, Electronics Research Duirec- 
torate, AF Cambridge Res. Center, Cam- 
bridge, Mass., Tech. Rep. 55-118, vii+16 
pp.; 1955.) The problem of simplifying 
switching circuitry is closely allied with 
the problem of simplifying truth functions, 
and has given rise to investigations of 
Boolean functions and Boolean matrices. 
The authors, beginning with the analysis 
by Quine (Amer. Math Monthly, vol. 59, 
pp., 521-531; 1952, MR 14, 440) have 
developed an algorithm which tests whether 
a Boolean function is equal to one. This 
“sum-to-one” process is based upon the 
following theorem: if a;, be, c: and R are 
Boolean functions not involving the Boo- 
lean variable s explicitly, then 


soo a,+s >) bet > cope R 
j k l 


if and only if Soj% ajbe+ Dor ZR. 

This sum-to-one algorithm is then ap- 
plied to find all minimal and irredundant 
Boolean sums equivalent to a given Boolean 
Function. It seems to the reviewer that the 
method of application of the basic algorithm 
to the main problem could be considerably 
clarified by, for example, attempting to 
program it for some automatic computer. 

S. Gorn 
Courtesy of Mathematical Reviews 


57-70 

A Note on Microprogramming— Herbert 

T. Glantz. (J. Assoc. Comp. Mach., vol. 3, 
pp. 77-84; April, 1956.) Since the term 
“microprogramming” was introduced early 
in this decade, the word has become widely 
used, but with many different meanings. 
Unfortunately, in most cases, the use of the 
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term has not occurred in the technical litera- 
ture—a situation which has retarded the 
development of a common definition 
throughout the field. Hence, it is consider- 
able service in itself that the author has 
chosen to publish his interpretation of this 

| technique. In discussing this concept. the 
author suggests that a possible criticism of 
| current commercial machines is that there is 
a certain lack of versatility which arises 
from the customary fixed set of commands 
which provide for the use of only a fraction 
of the total possible different sequences of 
| combinations of the internal operations of 
| the machine. On the other hand, the essen- 
tial characteristic of a microprogrammed 
| machine is that it permits the programmer, 
in one manner or another, to manufacture 
his own command which could produce 
| whatever sequence of combinations of inter- 
| nal operations he might desire. These arbi- 
trary sequences are produced in the example 
which he discusses by placing a particular 
serial-parallel pattern of binary characters 
in a magnetic core memory. A control com- 
mand which initiates the sequential read-out 
of this pattern into appropriate control 
circuits initiates the desired operation. Al- 
though it lacks any rigorous analysis of 
the possible advantages or disadvantages of 
such a system as compared to more conven- 
tional ones, and the definitions may be 
considered to be incomplete, this article 
represents a considerable contribution to the 
common understanding of this technique. 
Gordon Morrison 


| 
| 
; 


| 57-71 

Reliable Circuits Using Less Reliable 
Relays, Part I—E. F. Moore and C. E. 
Shannon. (J. Franklin Inst., vol. 262, pp. 
191-208; September, 1956.) The fuil paper 
describes an investigation made of relays 
whose reliability can be described in simple 
terms by means of probabilities. It is shown 
that by using a sufficiently large number of 
' these relays in the proper manner, circuits 
can be built which are arbitrarily reliable, 
regardless of how unreliable the original 
relays are. In Part I, it is shown that if 
relay contacts may be characterized by 
the probability a of some closing when the 
coil is energized, and with probability c 
of others closing when the relay is not ener- 
gized, then a highly reliable computer may 
be constructed from these relays if and only 
if ac. A set of theorems are proven which 
make possible the design of a network in 
which a series of contacts may be replaced 
by the network and have the probability 
h(a) of the contacts being open when ener- 
gized be less than $ and the probability h(c) 
of the contacts being open when not ener- 


gized be greater than 3. 


Gene M. Amdahl 


57-72 

Reliable Circuits Using Less Reliable 
Relays, Part II—E. F. Moore and C. E. 
Shannon (J. Franklin Inst., vol. 262, pp. 
281-297; October, 1956.) Part II develops 
a three-step procedure for designing a relay 
network of arbitrary reliability. The first 
step consists of finding a network which for 
a given a and c makes h(a) <3 and h(c) >}. 
The second step consists of finding networks 
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for which h[h(a) <3]<4 and h[h(c) >}] =3. 
The third step consists of finding networks 
which take these probabilities of 1 and 3 
or better and move them arbitrarily close 
to 0 or 1 respectively. The results of these 
procedures are compared with von Neu- 
mann’s results using the Sheffer stroke or- 
gans. For an individual contact (or Sheffer 
stroke organ) error probability of 0.005, 
to achieve a circuit error probability of 
about 10-9, a redundancy of about 100 to 1 
is required with contacts but about 60,000 
to 1 is required with Sheffer stroke organs. 

Gene M. Amdahl 


57-73 
Digital Process Control—M. L. Klein, 
F. K. Williams, and H. C. Morgan. (Instru- 
ments and Automation, vol. 29, pp. 1979- 
1984; October, 1956.) This paper discusses 
briefly the general problem of closed-loop 
automatic control systems. The problem of 
constructing a digital control system is 
presented. As an example a digital control 
system is designed to control a simple chemi- 
cal process in which it is required to control 
the temperature at one point and the pres- 
sure at two points. The solution is a fixed 
program machine. The material is well pre- 
sented; however, some mention might have 
been given of stored program machnies 
for ditigal control applications. 
Raymond Davis 


DIGITAL EQUIPMENT 


57-74 
Rechenautomaten im Dienste der Tech- 
nik. Erfahrungen mit dem Zuse-Rechen- 
automaten Z4—Eduard Stiefel. (Arbeitsge- 
meinschaft fiir Forschung des Landes Nord- 
rhein-Westfalen. No. 45, pp. 29-45, discussion 
pp. 47-65; 1955.) This is mainly a report on 
the work done and the experience gathered 
with the low-speed, small-memory relay 
computer Z4 when tabulating functions, 
solving polynomial equations, linear and 
nonlinear differential equations, partial dif- 
ferential equations, eigenvalue problems, 
etc. The author emphasizes that a general- 
purpose computer, used for problems of 
various length and complexity, will be useful 
only if its speed is matched by adequate in- 
and output facilities and if programming is 
as easy as possible. This has led to the fol- 
lowing requirements for the new machine 
(ERMETH) replacing Z4: 1) ERMETH 
is a decimal machine throughout with 
floating point arithmetics. 2) The “coding 
box” on which codes will be punched will 
show the conventional mathematical sym- 
bols (+, —, X, + etc.). 3) A few simple 
orders together with a “Small Manual” will 
suffice for the simple computations of “noy- 
ices,” but additional facilities will be avail- 
able for “experts.” After the lecture, H. 
Behnke, L. Brandt, G. Hoheisel, H. Kaiser, 
J. von Neumann, E. Sperner, E. Peschl, 
H. Peterson, R. Sauer, E. Stiefel, A. Wal- 
ther and W. Weizel discussed various aspects 
ranging from storage devices and analogy 
machines to economics and training of 
mathematical personnel. 
H. H. Goldstine 
Courtesy of Mathematical Reviews 
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57-75 
The G-15 Digital Computer—S. Herbert 
Lewis. (Instruments and Automation, vol. 
29, pp. 1773-1779; September, 1956.) This 
paper was presented at the Second Com- 
puter Clinic, Second International Automa- 
tion Exposition, Chicago, Illinois, 1955. It 
is an introduction to the principles, opera- 
tion, and application of a general-purpose 
digital computer. The author does a good 
job of explaining to the novice what a 
general purpose computer is supposed to do. 
Elementary binary arithmetic is presented 
in a clear manner. After this introduction 
the major elements of a general purpose 
computer are described in terms of the G-15 
machine. The command word structure is 
given; however, the list of the commands 
available in the G-15 are omitted. Program- 
ming techniques are well covered and a 
typical flow chart is presented. Subroutines, 
scaling, and optimizing access time are 
briefly explained. A highway construction 
problem is taken as an example of a problem 
which illustrates the capabilities of the G-15 
computer. The author has succeeded in 
writing a clear and understandable paper 
for the reader who has no previous knowledge 
of the digital computing field. 
Raymond Davis 


57-76 
The UNIVAC and UNIVAC Scientific— 
W. Allen and G. E. Smith. (The Computer 
Handbook, ed. by Milton H. Aronson, The 
Instruments Publishing Co., Pittsburgh, 
Pa., pp. 46-55; 1955.) Brief descriptions are 
provided of the functions and interrelation- 
ships of the major components of the 
UNIVAC and UNIVAC Scientific (former- 
ly ERA 1103) computing systems, including 
input-output and auxiliary equipment. The 
instruction repertories are described in gener- 
al terms but are not listed in detail. Some 
fields of application of the two computer sys- 
tems are indicated, emphasizing the com- 
mercial uses of the UNIVAC and the scien- 
tific uses of the UNIVAC Scientific. A few 
simple examples are presented of the use of 
some of the special coding features of the 
UNIVAC Scientific. 
E. Bromberg 


57-77 
Digital Time and Sequence Control— 
M. L. Klein, F. K. Williams and H. C. 
Morgan. (Justruments and Automation, vol. 
29, pp. 2403-2407; December 1956.) The 
paper deals with the application of digital 
techniques to time and sequence control 
problems. The circuits of a variable-interval 
automation control are clearly discussed. 
Raymond Davis 


57-78 

Digital Telemetering and Print-Out of 
Pipeline Information—M. T. Nigh. (Elec- 
trical Eng., vol. 76, pp. 216-219; March, 
1957.) In the past, common methods of 
telemetering pipeline information have been 
the time-impulse method and variation of 
the frequency of a signal proportional to the 
measured variable. These are analog meth- 
ods. Noise or signal interruption cause dis- 
turbing readings on dispatcher recorders. 
This article describes a system which con- 
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verts readings to digital form at the source 
and transmits binary-coded decimal infor- 
mation to a central station. The system 
includes provisions for rejecting nonsense 
information. A system now in operation sig- 
nals twelve different conditions at a remote 
pumping station, including seven malfunc- 
tion conditions, with a few of these mal- 
function indications really representing a 
set of several malfunctions. The station also 
signals when deviations of certain variables 
exceed preset amounts, and when the rate 
of change of certain variables exceeds pres- 
ent amounts, either case causing automatic 
print-out in an electric typewriter at the 
central station. The transmitters and en- 
coders are described in general terms—no 
circuits. Telephone transmission lines are 
used. Bits are transmitted as long or short 
pulses. These channels can be used for trans- 
mitting control signals back to the remote 
station. 

H. T. Larson 


UTILIZATION OF DIGITAL 
EQUIPMENT 


57-79 

Programming the Computer—S. M. 
Rock and W. W. Klammer. (Control Engi- 
neering, vol. 4, pp. 119-123; March, 1957.) 
This paper describes the various steps in- 
volved in preparing a program for scientific 
and, particularly, business problems. It men- 
tions the formulation and numerical analysis 
of scientific problems very briefly, but con- 
centrates mainly on the formulation and 
programming of business data processing 
problems. The article describes how the 
general analysis of a business problem might 
be performed using check reconciliation, as 
an example, specifying how the memory 
space must be assigned for the program and 
each type of input and how the detailed flow 
chart is prepared from which the individual 
program steps are written. The authors 
present an example of a program for the 
bank (check) reconciliation problem. They 
then describe debugging operations and the 
various alternatives available here; such 
as a program comparing computer results 
with manually obtained results in a checking 


‘routine, monitoring programs, and “post 
_mortems.” They describe ‘how the “B regis- 


ter” or “index register” is used, and briefly 
discuss subroutines and automatic program- 
ming. The article is an excellent, but brief, 
description of the operations involved in 
preparing programs for business applica- 
tions. 


Neal J. Dean 


57-80 

Office Automation—W. F. McCelland, 
D. R. Jarema, and R. M. James. (/nstru- 
ments and Automation, vol. 29, pp. 2426— 
2430; December, 1956.) The authors have 
divided this paper into three sections, a gen- 
eral discussion of linear programming and 
simulation, and two specific applications of 
the IBM 650 data processing system. The 
linear programming discussion is centered 
around the savings that might be achieved 
by applying this mathematical tool to busi- 
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ness problems. No attempt is made to ex- 
plain the techniques involved. Both of the 
installations, one at the A. C. Nielson Com- 
pany, and the other at Montgomery Ward, 
have replaced punched card equipment. No 
major procedural changes were made. The 
A. C. Neilson Company is processing mar- 
ket research data and Montgomery Ward 
is using the 650 to prepare operating budg- 
et and profit and loss statements. 

Raymond Davis 


BOOK REVIEWS 


57-81 

Analog Computer Techniques—C. L. 
Johnson. (McGraw-Hill Book Co., New 
York, N. Y., 246 pp.+6 index+11 appen- 
dix+xi; 1956.) The material covered in 
this book should be a valuable aid to users 
and operators of de electronic analog com- 
puting equipment. As the author points 
out, the analog computer and the digital 
computer become much more effective tools 
for the engineer, the scientist, and the 
mathematician when he understands the 
principles of operation of such equipment 
and has a working knowledge of the tech- 
niques involved in its utilization. Brief, but 
in general, adequate descriptions are given 
of the principles of operation and design of 
the computing components which comprise 
present day analog computers. In view of 
the rapid changes occurring in the state of 
the art, references made by the author to 
the characteristics of components in com- 
mercially available computers must be ac- 
cepted by the reader largely as examples for 
the point being made by the author, rather 
than as information on the capability of the 
latest available equipment. Good discus- 
sions are presented on how to program prob- 
lems including the considerations of time 
scale, amplitude scale, linear and nonlinear 
function representation, problem stability, 
and computer stability. One chapter presents 
the techniques involved in using an elec- 
tronic differential analyzer for special appli- 
cations such as solving algebraic equations, 
partial differential equations, and eigenvalue 
problems. Control and checking features 
included in present-day computers are sum- 
marized. The use of numerical analysis tech- 
niques and digital computers as auxiliary 
tools for analog computer checking is also 
covered. Two chapters at the end of the 
book are devoted to related-type computing 
equipment; there is a discussion of the prin- 
ciples of operation of the repetitive type 
analog computer and of the digital differen- 
tial analyzer. It should be noted that two 
or three areas of computer operation are not 
included in the book. Although the drift and 
bandwidth problems are touched on in the 
discussion of time scale, certain mathemati- 
cal and graphical information should be 
available to the operator so that he may 
more readily adapt his problem to the com- 
puter and also chose a computer for his 
problem. The latter consideration requires 
a discussion of using analog computers for 
real time simulation. More information 
should have been provided regarding co- 
ordinate systems and vector transforma- 
tions. It would have been helpful also for the 
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book to include a discussion of ac electro- 
mechanical differential analyzers and pas- 
sive network analyzers. The author is to 
be congratulated for his accomplishment of 
assembling in one place the techniques which 
are being used by most of the large analog 
computer operating groups, in some cases, 
as written instructions, but most often as 
operator trade secrets. Although the book: 
is written as a text for an introductory 
course on using analog computers, it should 
prove valuable to computer operators and 
users out of the class room. 
C. M. Edwards 
Courtesy of Proc. IRE 


57-82 
Automatic Digital Computers—M. V. 
Wilkes. (John Wiley and Sons, New York, 
N. Y.,-305 pp.-++-x; 1956.) This book is a 
good survey of the digital computing field. 
Being a survey, it will probably not satisfy 
each specialist in his own line, but it will 
give him a good understanding of the prob- 
lems arising in other lines. As the author 
states, the book will not enable a competent 
computer-programmer to design automatic 
digital computers, nor will it enable a com- 
puter designer to become a competent com- 
puter programmer. The book contains excel- 
lent sections on the history of automaticcom- 
puters, the design of many of the present 
computers, the philosophy of coding, and 
the philosophy of computer design. The 
author goes quite deeply into the electronic 
circuitry of computers and into coding tech- 
niques in order to delineate the still unre- 
solved controversies, such as serial against 
parallel design and automatic machine com- 
piling of programs against direct coding in 
machine language. The typography is clear 
and the illustrations both profuse and of 
good quality. One of the most valuable 
contributions is the comprehensive anno- 
tated bibliography. This book should cer- 
tainly be in the library of every large-scale 
computing facility and will be of real inter- 
est to all who are associated in any way with 
the digital computing field. 
John E. Maxfield 
Courtesy of American Mathematical Monthly 


57-83 

Handbook of Semiconductor Electronics 
—Edited by L. P. Hunter. (McGraw-Hill 
Book Co., New York, N. Y., 132 pp.+13 
index +68 bibliography +xxvili; 1956.) This 
book, which includes sections by fourteen 
contributors, is the result of an effort to 
collect together the major principles in the 
semiconductor electronics field. An idea 
of its scope can be obtained from the table 
of contents. There are twenty sections dis- 
tributed among four main divisions: physics 
of transistors, diodes, and photocells; tech- 
nology of transistors, diodes, and photo- 
cells; circuit design and application for 
transistors, diodes, and photocells; and refer- 
ence material. Section headings are: 1) 
Transistor Characteristics; 2) Electronic 
Conduction in Solids; 3) Rectification in 
Solids; 4) Transistor Action; 5) Photocon- 
ductivity and Photovoltaic Cells; 6) Prep- 
aration of Semiconductor Materials; 7) 
Methods of Preparing PN Junctions; 8) 
Metal-Semiconductor Contacts; 9) En- 
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capsulation; 10) Device Design Considera- 
tions; 11) Low-Frequency Amplifiers; 12) 
High-Frequency and Video Amplifiers; 13) 
Directly Coupled Amplifiers; 14) Transistor 
Oscillators; 15) Transistor Switching Cir- 
cuits; 16) Circuits Using Special Semicon- 
ductor Devices; 17) Graphical Analysis of 
Nonlinear Circuits; 18) Matrix Methods of 
Circuit Analysis; 19) Measurement of Semi- 
conductor-Device Parameters; 20) Meas- 
urement of Material Parameters. An exten- 
sive bibliography of the transistor, semi- 
conductor, and rectifier fields covering the 
period 1940-1955 inclusive closes the book; 
this bibliography is the most comprehensive 
listing in print for these fields. As its title 
indicates, the book is organized as a hand- 
book, simplifying the location of desired 
information but carrying with it the dis- 
advantage of small type, the same size as 
that of this review. There are very few 
typographical errors. The excellent selection 
and balance of topics, and the lack of ap- 
appreciable duplication of material in sec- 
tions written by different contributors are a 
tribute to the editor’s skill and patience. 
All the contributions are authoritative and 
reflect the wide experience of their authors. 
Especially noteworthy, however, are the 
initial sections on the physics of semicon- 
ducting devices, the section on low-frequency 
amplification, and that on transistor switch- 
ing circuits. Many of the sections de- 
part from being strictly a collection and dis- 
cussion of past work and in addition include 
valuable material unavailable elsewhere. It 
will be a rare worker in the semiconductor 
field who will not find much to learn from 
the book. Nevertheless, its topics are so well 
developed from basic principles that the 
newcomer to semiconductors may use it to 
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familiarize himself with almost any subject 

in the area. The book fulfills its aims very 

well, may be highly recommended, and is 

likely to remain for some time the standard 
reference in this rapidly growing field. 

J. R. Macdonald 

Courtesy of Proc. IRE 


57-84 

Handbook of Industrial Electronic Con- 
trol Circuits—J. Markus and V. Zeluff. 
(McGraw-Hill Book Co., New York, N. Y., 
352 pp. +8 index +xiv; 1956.) Engineers who 
have made use of the 1948 edition of Markus 
and Zeluff’s Handbook of Industrial Elec- 
tronic Circuits will be delighted to know 
that a sequel, Handbook of Industrial Elec- 
tronic Control Circuits, is now available. This 
handbook, the fifth by these electronics 
editors, fully lives up to and probably will 
somewhat enhance the reputation estab- 
lished by its predecessors. As electronics 
engineers well know, much that is worth- 
while in the way of industrial electronics 
circuits has been published since 1948 but 
the problem of finding just what one wants 
when it is needed is a burdensome one and 
this reviewer for one finds it utterly im- 
possible to remember all the circuits that 
are of frequent use, much less those that in- 
volve only occasional use. This book then 
brings together in conveniently indexed 
form most (perhaps all) of the circuits 
published in Electronics during the years 
1948 to 1956 which are of primary interest 
to the industrial electronics engineer. As in 
its predecessor volume, each circuit is ac- 
companied by a concise description giving 
the general nature of the circuit, an explana- 
tion of how it works, data on critical com- 
ponents, and suggested applications. New 
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circuits are included under most of the main 
subject titles used in the 1948 book and four 
new subjects—radiation circuits, remote con- 
trol circuits, sorting circuits and transistor 
circuits—are introduced. Transistor circuits 
in particular, were just not in existence 
prior to 1948 but forty circuits are included 
in this issue and the list seems to encom- 
pass just about everything suitable for use 
in the industrial field, plus some that are 
probably of use solely in communications 
work. Even so, the transistor application 
field is so new that further expansion of this 
section would be a very fine thing. It is in- 
teresting to note that no limiter circuits were 
published between 1948 and 1956, but in 
such a well-developed field as motor control 
twenty-three new circuits are shown. This 
handbook, like its predecessor, is intended 
primarily for experienced engineers working 
in the industrial electronics field. Instructors 
of industrial electronics courses will find it 
valuable as a reference book and recent 
graduates will find it of great help in bridg- 
ing the gap between college theory and 
practical circuit design. This reviewer was 
disappointed not to find at least one good 
memory circuit for use in automatic test 
equipment, also a good thermocouple tem- 
perature controller circuit, but apparently 
these are yet to appear anywhere. In a few 
cases, notably the neon lamp flip-flop circuit 
on page 77, the description is oversimplified 
to the point of not warning about change of 
characteristics of the neon tube with time; 
hence, it may lead to disappointment if fol- 
lowed too literally. All in all, the book is 
extremely well done and fills a very sub- 
stantial and rapidly growing need in its field. 
C. A. Priest 

Courtesy of Proc. IRE 
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